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Abstract 



We study the effect of global topology of tlie spatial geometry on the cosmic mi- 
crowave background (CMB) for closed flat and closed hyperbohc models in which 
the spatial hypersurface is multiply connected. If the CMB temperature fluctua- 
tions were entirely produced at the last scattering, then the large-angle fluctuations 
would be much suppressed in comparison with the simply connected counterparts 
which is at variance with the observational data. However, as we shall show in this 
thesis, for low matter density models the observational constraints are less stringent 
since a large amount of large-angle fluctuations could be produced at late times. 
On the other hand, a slight suppression in large-angle temperature correlations in 
such models explains rather naturally the observed anomalously low quadrupole 
which is incompatible with the prediction of the "standard'' Friedmann-Robertson- 
Walker-Lemaitre models. Interestingly, moreover, the development in the astro- 
nomical observation technology has made it possible to directly explore the imprint 
of the non-trivial topology by looking for identical objects so called "ghosts" in 
wide separated directions. For the CMB temperature fluctuations identical patterns 
would appear on a pair of circles in the sky. Another interesting feature is the 
non-Gaussianity in the temperature fluctuations. Inhomogeneous and anisotropic 
Gaussian fluctuations for a particular choice of position and orientation are regarded 
as non-Gaussian fluctuations for a homogeneous and isotropic ensemble. If adiabatic 
non-Gaussian fluctuations with vanishing skewness but non-vanishing kurtosis are 
found only on large angular scales in the CMB then it will be a strong sign of the 
non-trivial topology. Even if we failed to detect the identical patterns or objects, 
the imprint of the "finiteness" could be still observable by measuring such statistical 
property. 



Contents 



I introduction 





1.1 Finite or infinite?] 




1.2 Imprint ot topology 






2 


3-dimensional Topology 




2.1 Flat geometry 




2.2 Spherical geometry 




2.3 Hyperbolic geometry] 



3 Mode Functions 

3.1 Boundary element method 

p. 1.1 Formulatioii 

3.1.2 Computation ot low-lying eigenmodes . . 

3.1.3 Statistical property of eigenmodes .... 

3.2 Periodic orbit sum method 

3.2.1 Length spectra 

3.2.2 'I'race formula 

p. 2. 3 First eigenvalue and geometrical quantity 

3.3 Spectral measurements of global anisotropy . . . 



4 CMB Anisotropy 

4.1 Linear perturbation theory . . . 

4.2 Sachs -Wolfe cfi'ect| 

4.3 Angular power spectrum .... 

4.3.1 Closed flat model^ . . . 

4.3.2 Closed hyperbolic models 
O Bayesian analysis 

[4.4.1 Closed flat moclcl^ . . . 
4.4.2 Closed hyperbolic models 



5 Search for Topology 

5.1 Circles in the sky . . . 
Mon-Caussian signature 



ii 



5.2.1 Cosmic variance 



5.2.2 'ibpologicai measurements 



|6 Summary 



A Boundary integral equation 

B Algorithm of computing length spectra 

|C Estimate of volume in terms of diameter 



83 
85 

90 

100 

102 

103 



iii 



List of Figures 



2.1 Volume Spectra of CH manifolds 


2.2 Cusped Hyperbolic Manifold 




2.3 Smallest CH Manifolds 





3.1 1168 Boundary Elements 

3.2 Lowest Eigenmodes 

3.3 Error Estimation 

3.4 Random Nature of Expansion Coefficients! 

3.5 Z'^ Symmetry 

3TB Kolomogorov-Smirnov 'lest 111 

3.7 Variance of Expansion Coefficients 

O Classical Staircases 

3.9 Locally Averaged Multiplicities 

3.1(J Spectral Staircases 

3.11 Diameter vs 1^'irst Eigenvalue 

3.12 Diameter vs Volume 

3.13 Volume vs I'lrst Eigenvalue 

3.14 Length of Shortest Periodic Orbit vs fci or Diameter 

3.15 Deviation from Weyl's I'brmula 

3.16 Spectrum Measurements 

3.17 Anisotropy in Ceometry 



4.1 COBE-DMR Anisotropy Sky Mi^ 

O Suppression in Large-angle Power for 'ibroidal Models 

4.3 Ceometric Variance 

4.4 Comparison with COBE data 

4.5 Suppression in Large-angle Power for CH models . . 

4.6 Relative Likelihoods of Toroidal Models I 

4.7 Relative Likelihoods of 'Ibroidal Models 11 

4.8 Relative Likelihoods of Smallest CH Models 

O Relative Likelihoods of Open and 'i'hurston Model . . 

4.10 Relative Likelihoods of 27 CH Models! 



Matched Circles 



iv 



5.2 Detectability in COBE Data 

5.3 Detectability in IVlAJr' Data) 

5.4 Non-Gaussian Distribution 

5.5 Deviation from Distribution 

5.6 Cumulative Distribution 

5.7 Cosmic Variance .... 

5.8 Simulated CIVIB Sky Map 

5.9 Minkowski Functional . 

5.10 JNon-Gaussian Signature in Variance 



A.l Boundary Integral 



C.l Estimate of Volume of "thin" Part 



V 



List of Tables 



3.1 Eigenvalue k and Multiplicity mlk) 

3.2 Kolmogorov-Smirnov 'lest 1 

3.3 Kolmogorov-Smirnov Test 11 

3.4 Run Test 

3.5 Coefficients of l:''itting Curves for Averaged Multiplicities 

3.6 First Eigenvalues 



vi 



vii 



1 



Chapter 1 
Introduction 



And so I'll follow on, and whereso 'er thou set the extreme coasts, I'll query, 
"what becomes thereafter of thy spear?" 'Twill come to pass that nowhere can 
a world's- end he, and that the chance for further flight prolongs forever the 
flight itself. 

-De Rerum Natura(Lucretius, 98?-55? BC) 

1.1 Finite or infinite? 

To ancient people the sky had been appeared as the immense dome with stars placed 
just inside or on it. The daily movements of the stars were attributed to the rotation 
of the dome called the celestial sphere around the Earth. Aristotle concluded that 
the size of the celestial sphere must be finite since the movement of an object of 
infinite size is not allowed by his philosophy. He denied the existence of any forms 
of matter, space and time beyond the boundary. On the other hand Epicurus has 
thought that the universe has no beginning in time and the space is unlimited in 
size. If the universe were limited in size, he said, one could go to the end, throw a 
spear and where the spear was located would be the new 'limit' of the universe. 

In the 16th century Copernicus proposed a new cosmology where the Sun is 
placed at the center of the universe instead of the Earth. Historically, the helio- 
centric cosmology had already been proposed by Aristarchus in the third B.C but 
soon the idea was rejected under the strong influence of Aristotle's philosophy. The 
failure to detect the parallax at those days must have provided a ground for the old 
geocentric (Earth-centered) cosmology. Copernicus' new idea remained obscure for 
about 100 years after his death. However, in the 17th century, the works of Kepler, 
Galileo, and Newton built upon the heliocentric cosmology have completely swept 
away the old geocentric cosmology. The parallax effect owing to the movement of 
the Earth has been finally confirmed in the 19th century. Thus the argument of 
Aristotle about the finiteness of the universe lost completely its ground. It seems 
that most of people in modern era are comfortable with the concept of an infinite 
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space which extends forever as far as we could see. 

However, one should be aware that the no-boundary condition does not require 
that the space is unlimited in size. Let us consider the surface of the Earth. If the 
surface were flat, one might worry about the "end" of it. Actually there is no end 
since the surface of the Earth is closed. In the 19th century, Riemann proposed a 
new cosmological model whose spatial geometry is described by a 3-sphere, a closed 
3- space without boundary. If we were lived in such a space, an arrow which had 
been shot in the air would be able to move back and hit the archer. This could have 
surprised Epicurus who did not know the non-Euchdean geometry. 

In the standard Friedmann-Robertson-Walker-Lemaitre (FRWL) models, the 
spatial geometry is described by homogeneous and isotropic spaces, namely, a 3- 
sphere 5*^, an Euclidean space and a hyperbolic space which have positive, 
zero and negative constant curvature, respectively. The former one is spatially 
closed (finite) but the latter twos are spatially open (infinite). The standard models 
have succeeded in explaining the following important observational facts:the ex- 
pansion of the universe, the content of light elements and the cosmic microwave 
background (CMB). However, we are still not sure whether the spatial hypersurface 
is closed (finite) or not (infinite). It seems that measuring the curvature of the spatial 
geometry is sufficient to answer the question. 

In fact, the above three spaces might not be the unique candidates for the cos- 
mological model which describes our universe. We can consider a various kinds of 
closed multiply connected^ spaces with non-trivial global topology other than 0. 
For instance, a flat 3-torus which is locally isometric to E^ can be obtained by 
identifying the opposite faces of a cube. Identifying antipodal points in yields a 
projective space RP^. Furthermore, a plenty of examples of closed hyperbolic (CH) 
spaces have been known. It would be inappropriate to represent hyperbolic (=con- 
stantly negatively curved) spaces as "open" spaces which have been widely used in 
the literature of astrophysics. The local dynamics of a closed multiply connected 
space are equivalent to those of the simply connected counterpart since Einstein's 
equations can only specify the local structure of spacetime and matter. Therefore 
these multiply connected models do not confront with the above-mentioned three 
observational facts. One can also consider non-compact finite-volume spaces which 
extends unlimitedly in particular directions. However, in what follows we shall only 
consider spatially "closed" (compact) models just for simplicity. Regardless of each 
topology, the effect of the non-trivial topology appears only on scales of the order 
of the length of the closed curves that are not continuously contractable to a point. 



^If there exists a closed curve which cannot be continuously contracted to a point then the 
space is multiply connected otherwise the space is simply connected. 

^It has been conjectured that is the only example of a closed connected and simply connected 
3-space (Poincare's conjecture). 
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1.2 Imprint of topology 

Suppose that the spatial hypersurface is multiply connected on scales of the order 
of the horizon or less. In this case, the universe is called the "small universe" |I| for 
which we would be able to observe surprising periodical patterns in far distant place 
as if we were looking into a kaleidoscope. Apparently it would look like the space 
extending to infinitely remote points, but the observed images represent a series of 
"snapshots" of one astronomical object at simultaneous or different epoch as the 
photons go around the space. Note that for a model without the cosmological 
constant, we would not be able to see such images since the universe would be col- 
lapsed by the time the photons go around the space. 

Surprisingly the recent advance in the observation technology has made it pos- 
sible to determine the spatial global topology of the universe. The observation 
methods are divided broadly into two categories which are complementary: the de- 
tection of the periodical structure in the astronomical objects and the detection of 
the peculiar property in the statistics of density or CMB temperature fluctuations. 
The former method can directly prove the presence of the non-trivial topology but it 
needs precise measurements of relevant physical quantities. On the other hand, the 
latter methods do not require high-precision measurements while the direct proof 
might be difficult since it relies on a certain assumption of the initial conditions. 

In multiply connected spaces a number of geodesic segments that connects ar- 
bitrary two points exist. We call the nearest image "real" for which the geodesic 
segment that connects the light source and the observer is the shortest, and oth- 
ers are called "ghosts". If the comoving radius of the last scattering surface were 
much larger than the injectivity radius (i.e. half of the shortest length of a loop 
which cannot be continuously contracted to a point) in the space theoretically we 
would be able to find such ghost images. In general, the lengths of these geodesic 
segments are different. The longer segments correspond to older images. Therefore, 
the life time of the light source must be long and stable for recognizing the images 
as "ghosts". To date a number of authors have tried to search for "ghosts" using 
cluster catalogs(see |0 and references therein) but have failed to obtain any positive 
signals. 

The appearance of the "ghosts" is relevant to the periodical structure in the 
distribution of observed objects since each geodesic segment corresponds to a copy 
of one single domain called the fundamental domain which tessellates the developed 
space (the apparent space) called the universal covering space. One way for detect- 
ing the periodical structure is to look for the spike in the distribution of the physical 
distance between two arbitrary objects[^ Note that the method is applicable 
to only flat and some spherical geometry. Using a 3-dimensional catalog of galaxy 
clusters by Bury a constraint L>600/i~^Mpc for closed flat 3-torus models in which 
the fundamental domain is a cube with side L without cosmological constant has 
been obtained 0. If the geodesies that connect the light source and the observer 
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have equivalent lengths, then we can see perfectly identical objects without time 
delay in wide separated directions in the sky. For instance, one way is to look for 
any identical sets of neighboring QSOs[Q. If the relative position and the orienta- 
tion of such QSOs were equivalent to those of another set of QSOs that had been 
observed in a different direction then it would be a sign that we were looking at one 
object (neighboring QSOs) from different directions. In general such a special con- 
figuration in which one observes an identical object at a simultaneous time appears 
on a pair of circles. Suppose that the observer sits at the center of a sphere. If the 
sphere is large enough then it wraps around the space and intersects with itself on 
a circle. Then one would easily notice that the two geodesic segments that connect 
the center with a point on the circle have the same length. All we have to do is 
to search for a pair of identical circles in the sky. One might use QSOs or another 
high-z objects such as 7-ray bursts for detecting the perfectly identical "ghosts" but 
ideal one is to use the CMB temperature fluctuations]^. The so-called "circle in the 
sky" method is a powerful tool for detecting the non-trivial topology of the spatial 
geometry if the space is so small that the last scattering surface can wraps around 
the space. 

Another signature of the non-trivial topology appears in the statistical prop- 
erty of the CMB temperature fluctuations. The break of the global homogeneity 
and isotropy in the spatial geometry naturally leads to a non-Gaussian feature in 
the fluctuations regardless of the type of the primordial perturbation (adiabatic or 
isocurvature). For a particular choice of orientation and position of the observer, 
the fluctuations form an inhomogeneous and anisotropic Gaussian random field. 
Marginalizing over the orientation and the position, then the ensemble of fluctua- 
tions can be regarded as a homogeneous and isotropic non-Gaussian random field. 
The non-Gaussian signals might appear in the topological quantities (total length 
and genus) of isotemperature contours in the sky map although the detection of 
the signal using the COBE data|^ is unlikely owing to the significant instrumental 
noises. However, the observed non-Gaussian signals in the bispectrum might be 
the positive sign of the non-trivial topology if they are not related with any system- 
atic errors. 

Using the COBE data, a lower bound L>4800 h^^Mpc for compact flat 3-torus 
toroidal models (for which the fundamental domain is a cube with equal sides L) 
without the cosmological constant (assuming that the initial perturbation is adia- 



batic and scale-invariant (^ = 1) ) has been obtained ||T0|, |TT], |T^, |T^. The suppression 



of the fluctuations on scales beyond the size of the fundamental domain leads to a 
decrease in the large-angle power spectra Ci which is at variance with the observed 
temperature correlations. 

In contrast, for low matter density models, the constraint could be considerably 
milder than the locally isotropic and homogeneous flat (Einstein-de-Sitter) models 
since a bulk of large-angle CMB fluctuations can be produced at late epoch by the 
so-called (late) integrated Sachs- Wolfe (ISW) effect fl^, |l5l which is the gravita- 
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tional blueshift effect of tfie free streaming pliotons by tlie decay of tlie gravitational 
potential. Because the angular sizes of the fluctuations produced at late time are 
large, the suppression of the fluctuations on scale larger than the topological iden- 
tification scale may not lead to a significant suppression of the large-angle power if 
the ISW effect is dominant. 

Thus the investigation of the CMB anisotropy in low matter density models 
is quite important. If there is no cosmological constant, then the geometry must 
be hyperbolic, the observational aspects of a limited number of closed(compact) 
hyperbolic (CH) models have been studied by Gott[|TB] and Fagundes [0, [T^, |TU| 
but the property of the CMB in these models has not been investigated until re- 
cently although they have a number of interesting features. For instance, we expect 
that the initial perturbations are smoothed out since the geodesic flows are strongly 
chaotic This may provide a solution to the pre- inflationary initial value 

problem[0]. Another interesting property is the existence of the lower bound for the 
volume. It is known that the volume of CH manifolds must be larger than 0.16668 
times cube of the curvature radius although no concrete examples of manifolds with 
such small volumes are known [^]. If the deep fundamental theory predicts higher 
probability of the creation of the universe with smaller volume then we may answer 
the question why we see the periodical structure at present. 

Unfortunately, the simulation of the CMB anisotropy in CH models is not an 
easy task. Unlike closed flat spaces, the detailed property of the eigenmodes of the 
Laplacian has not yet been unveiled. One of the object of this thesis is to establish 
a concrete method for analyzing the eigenmodes of CH 3- spaces necessary for sim- 
ulating the CMB anisotropy which enable us to extract a generic property of the 
eigenmodes. For this purpose we consider two numerical methods which had origi- 
nally been developed for the analysis of semiclassical behavior in classically chaotic 
systems. 

In chapter 2, we briefly review the necessary ingredients of mathematics of 
the three types of closed locally homogeneous and isotropic spaces. In chapter 3 
we describe two numerical methods, namely, the direct boundary element method 
(DBEM) and the periodic orbit sum method (POSM) for calculating eigenmodes of 
the Laplacian of CH manifolds. Then we apply the methods to a number of known 
examples of CH manifolds with small volume and study the statistical property of 
eigenmodes. We also characterize the geometric property of manifolds in terms of 
the eigenvalue spectrum. Chapter 4 is devoted to the study of CMB anisotropy 
in closed low mater density models with flat and hyperbolic geometry. First we 
estimate the degree of the suppression in the large- angle power. Next, we carry out 
the Bayesian analysis using the COBE-DMR 4-year data and obtain observational 
constraints on the models. In chapter 5, we discuss methods for detecting non-trivial 
topology in the CMB. First we describe the direct method for searching identical 
fluctuation patterns then we explore the expected non-Gaussian signatures. In the 
last chapter, we summarize our results and discuss future prospects. 
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Chapter 2 

3-dimensional Topology 



. . . / think it is fair to say that until recently there was little reason to expect 
any analogous theory for manifolds of dimension 3 (or more)-except perhaps 
for the fact that so many 3-manifolds are beautiful. 

(W.P. Thurston) 



The classification of tlie 2-dimensional closed surfaces had already been completed 
in the 19th century. A closed surface is homeomorphic to either a sphere or a 2-torus 
with or without handles if orientable. A 2-torus without handles and that with han- 
dles can be endowed with flat geometry and hyperbolic geometry, respectively. In 
other words, 2-dimensional closed (orientable) surfaces can be classified into three 
types of geometry of homogeneous and isotropic spaces, namely a 2-sphere 5^, a 
Euclidean plane E"^ and a hyperbolic plane H^. 

Similarly, it has been conjectured that the topology of closed 3-manifolds (com- 
plete and without any singularities) can be classified into eight types of geometry 



of homogeneous spaces]^. Suppose a closed 3-manifold M. First we repeatedly 
cut M along 2-spheres embedded in M and glue 3-balls to the resulting boundary 
components. Next we cut the pieces along tori embedded nontrivially in M. The 
above-procedure is called a canonical decomposition. Thurston conjectured that ev- 
ery closed complete 3-manifold M has a canonical decomposition into pieces which 
can be endowed with a complete, locally homogeneous Riemannian metric. Let us 
define a geometry to be a pair (X, G) where X is a 3-manifold and G acts transi- 
tively on X such that the stabilizer of any point x G X is a compact subgroup of 
G. A geometry is called maximal if G is not contained in a larger group G' for a 
given X. By considering G-invariant metrics on X, one can recover the ordinary 
viewpoint of differential geometry (the spaces are homogeneous Riemannian mani- 
folds). (X, G) and (X', G') are equivalent if there is a diffeomorphism : X — X' 
with an isomorphism from G to G' (f) : g i — > (po go g^^. Then we have the following 
theorem proved by Thurston 
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Any maximal simply connected geometry which admits a compact quotient is equiva- 
lent to one of the following eight geometries {X, IsomX) where X is either E^., H^, S^, 
^2 X R,H^ X R, SL{2, R), Nil, or Sol (see [|6| for full details). Thus in addition to 
three geometries of constantly curved spaces E^, S^, (flat, spherical, hyperbolic), 
there exist five other geometries of homogeneous spaces. 

Although it is not known whether the canonical decomposition of a closed 3- 
manifold M into pieces of locally homogeneous spaces always exist, we do know 
how to construct every closed 3-manifold. Suppose a link L to be a 1-dimensional 
compact submanifold. By removing a regular neighborhood of L and gluing it back 
by some new identification determined by a diffeomorphism of the torus for each 
component of L, one can obtain a new manifold A^. The procedure is called the 
Dehn surgery. It has been proved that every closed 3-manifold can be obtained by 
a Dehn surgery along some link in a 3-sphere whose complement is hyperbolic. 
Suppose L C M is a link whose complement M — L can be endowed with hyper- 
bolic geometry. Then except for a finite number of choices, all the remaining Dehn 
surgeries yield closed hyperbolic manifolds [^]. This fact suggests that most closed 
3-manifolds are hyperbolic in some sense. 

In this chapter we describe 3-types of geometry of constant curvature E^, and 

which are most relevant to the cosmology. 



2.1 Flat geometry 

First we consider the fiat geometry (M, G) = {E\IS0{3) = x 50(3)). The ac- 
tions of the isometry group on E^ consist of an identity, a translation, a glide reflec- 
tion, i.e. a reflection in a plane through the origin followed by translation parallel to 
the plane and a helicoidal motion, i.e. a translation accompanied by a rotation, and 
their combination, which generate 18 distinct types of locally Euclidean spaces [P^]. 
Eight types of these spaces are open and other ten types are closed. The closed 
orientable ones are the following six types T^, T^/Z2, T^/Z4, T^/Z2 x Z2, T^/Z^ and 
T^/Zq where denotes a 3-torus and Zm is a cyclic group of order m. The funda- 
mental domains of the first four spaces can be a cube or a parallelepipe. Identifying 
the opposite faces by three translations one obtains a 3-torus T^. Similarly, identi- 
fying opposite faces by translations but one pair being rotated by vr or 7r/2 yields 
T^/Z2 and T^/Z^, respectively. Identifying the opposite faces by translations with 
a IT rotation, then we obtain T^/Z2 x Z2. The fundamental domain of the last two 
can be a hexagonal prism. Identifying the opposite faces by translations but the 
hexagonal ones being rotated by 27r/3 or n/3 yield T^/Z^ and T^/Zq, respectively. 
It should be noted that there is no upper or lower bound for the volume of the closed 
manifolds since there is no particular scale in flat geometry. 

geometry is said to admit a compact quotient if G has a subgroup T which acts in X as a 
covering group so that X/T becomes compact 
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T^ is globally homogeneous but the others are not globally homogeneous. If a 
Riemannian manifold M is globally and locally homogeneous then every element 7 
of the discrete isometry group F C G is equivalent to a Clifford translation of M, 
i.e. a transformation 7 : x — > 7(x) such that the distance between x and its image 
7(0;) is the same for all x G M. Thus the simplest model is belonging to a rather 
special class of multiply connected spaces. 



2.2 Spherical geometry 

Next we describe the topology of the spherical geometry (S''^, 5*0(4)). The most 
direct way to define is as the unit sphere {x^ + x^ + xl + xl = 1} in H'^ but here 
we consider the unit sphere + \z2\'^ = 1 in C^. Then S'^/Zm, (rn > 1) is called 
the lens space Lm which is defined as an orbit space of the 3-sphere under the 
action of Zm, given by p8| , 

(zi, Z2) I — > (zie^"^/"^, 226^"^/"^). (2.1) 

In the case of m = 2 the space is called the real projective space denoted as RP^. 
Each piece of a cell decomposition of of Lm consists of points {zi, Z2) G such that 
Z2 = pe"^ , p > and 

?^ < ^ < Mi±i), (2.2) 

m m 

where q = 0, 1, ... ,m — 1 specifies the cell. In stereographic projection, each cell 
looks like a convex lens. 

Identifying as the space of quaternions Q we can obtain another representa- 
tion of S^/T. The basis of Q is denoted as k} which satisfy 

z2 = f = e = -i (2.3) 

and 

ij = k = —ji, jk = i = —kj, ki = j = —ik. (2-4) 

The subspace spanned by the identity 1 is identified with R, and its elements are 
called rea/;the subspace Ri + Rj + Rk is called the space of pure quaternions. The 
conjugate of a quaternion q = a + bi + cj + dk is given hj q = a — bi — cj — dk and 
its square root is the absolute value of q denoted by The set 

Q' = {qeQ:\q\ = i} (2.5) 

forms a multiplicative group, called the group of unit quaternions which can be 
identified as 5*^. The elements of Q' are called the unit quaternions. It is known 
that the right or left multiplication by a unit quaternion gives a self-action of by 
orientation-preserving isometries while conjugation gives a self-action which takes 
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any two-sphere onto itself (which fixes 1)|2^. Conjugation can be described as 
follows. Let us consider q, q' G Q'. Then the action by conjugation is given by 
p{q)W) = qq'q~^ = p{—(i){q') which fixes 1. In other words, it leaves invariant a 
2-sphere. Thus p{q) can be regarded as an element of 50(3) but p : i — > 5*0(3) is 
a two to one {p{q) = p{—q)) map (homomorphism) . Therefore, can be regarded 
as the universal covering of 50(3), i.e. 5^/{±l} is isomorphic to 50(3). 

Thus the discrete isometry group of 50(4) (without fixed points)is described by 
the discrete subgroups of the rotations 50(3) which have fixed points, namely the 
cyclic group Zm(order=m), the dihedral group Dm(order=2m) (i.e. the symmetry 
group of a regular m-gon lying in a 2-plane), three polyhedral groups, i.e., the 
symmetry groups T(order=12) of a regular tetrahedron, 0(order=24) of a regular 
octahedron and /(order=60) of a regular icosahedron in R^. The binary dihedral 
and binary polyhedral groups are defined by 

D:,^P'\DJ, T*^p-\T), 0*^p-\0), r^p-\l). (2.6) 

Globally homogeneous spherical spaces are classified as follows: (i) 5^, (ii) RP^, (iii) 
5VZ^, (m > 2), (iv) S^/D*^, (m > 2), (v) 5VT*,(vi) S^/0\ (vii) S^/P Note 
that all the groups are Clifford groups. On the other hand, globally inhomogeneous 
spaces are obtained by the quotient of 5^ by a group of the following: (i) y<H{m > 
2) where H is either Z)*, T*, O* or /* and m is relatively prime to the order of H, (ii) 
a subgroup of index(i.e. a number of left(right) cosets with respect to the subgroup)3 
in Z^rn X T*, where m is odd (iii) a subgroup of index 3 in Z2n x -D2m; where n is even 
and m and n are relatively prime. Note that any isometry of that has no fixed 
point is orientation-preserving. Therefore even if we consider the maximal group 
0(3) instead of 50(3) the result does not change ||29[ . 
The volume of M = S^/T is simply given by 

vol(M) = 2tt^R^/\T\ (2.7) 

where R is the curvature radius and |r| is the order of F. Thus the largest manifold 
is itself. On the other hand, there is no lower bound for the volume since one 
can consider a group with arbitrary large number of elements (for instance, consider 
S^/Zfn with arbitrary large m). 



2.3 Hyperbolic geometry 

Finally we consider the hyperbolic geometry [H^, PSL[2,C)). The discrete sub- 
group F of PSL{2, C) which is the orientation-preserving isometry group of the 
simply-connected hyperbolic 3-space is called the Kleinian group. Any CH 3- 
spaces (either manifold or orbifold) can be described as compact quotients M = 
H^/T. If we represent as an upper half space {xi,X2,X3), the metric is written 
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as 

where i? is the curvature radius. In what follows, R is set to unity without loss of 
generality. If we represent a point p on the upper-half space, as a quaternion whose 
fourth component equals zero, then the actions of PSL{2, C) on U C U {00} take 
the form 

f ap + b . / \ 

: p p = ad — be = 1, p = z + X3J, z = Xi + X2I, (2.9) 

cp + d 

where a, b, c and d are complex numbers and 1, i and j are represented by matrices, 

-(o?)- -(o-O' ^K--)- 

The action 7 is explicitly written as 

7:i7^UCU{oo} ^ if^UCU{oo}, 

, , / (az + b) (cz + d) + acxl x^ \ , ^ 

where a bar denotes a complex conjugate. Elements of V for orientable CH manifolds 
are SL{2, C) conjugate to 

/exp(//2 + i0/2) 
\ exp(-//2 - i(p/2) 

which are called loxodromic if </> 7^ and hyperbolic if = 0. An orientable CH 
manifold is obtained as a compact quotient of by a discrete subgroup which 
consist of loxodromic or hyperbolic elements. 

Topological construction of CH manifolds starts with a cusped manifold with 
finite volume Mc obtained by gluing ideal tetrahedra. Mc is topologically equivalent 
to the complement of a knot K or link L (which consists of knots) in 3-sphere S'^ 
or some other closed 3-spaces. A surgery in which one removes the tubular neigh- 
borhood N oi K whose boundary is homeomorphic to a torus, and replace by 
a solid torus so that a meridian^ in the solid torus goes to {p, q) curv^ on A^ is 
called (p, q) Dehn surgery. Except for a finite number of cases, Dehn surgeries on 

^Given a set of generators a and b for the fundamental group of a torus, a closed curve which 
connects a point x in the torus with ax is called a meridian and another curve which connects a 
point X with bx is called a longitude. 

^If C connects a point x with another point (pa + qb)x where p and q are co-prime integer, C 
is called a (p, q) curve. 
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Figure 2.1: Volume Spectra of CH manifolds 

A prefix "m" in the labeling number represents a cusped manifold obtained 
by gluing five or fewer ideal tetrahedra. The numbers in the right side 
denote the volumes of the corresponding cusped manifold. The volumes 
have been computed by using the SnapPea kernel. 



K always yield CH 3-manifolds which implies that most compact 3-manifolds are 
hyperbolic ||25||. However, we are not still sure in what Dehn surgery yields 

a hyperbolic manifold for any kinds of cusped manifolds. Furthermore there are a 
variety of ways to construct isometric CH manifolds by performing various Dehn 
surgeries on various cusped hyperbolic manifolds. Thus the classification of CH 
manifolds has not yet been completed. 

However, with the help of the computer, we now know a number of samples 
of CH manifolds(and orbifolds). For instance a computer program "SnapPea" by 
Weeks can numerically perform Dehn surgeries on cusped manifolds which have 
made it possible to construct a large number of samples of CH 3-manifolds. Snap- 
Pea can also analyze various property of CH manifolds such as volume, fundamental 
group, homology, symmetry group, length spectra and so on. 

In order to describe CH 3-manifolds, the volume plays a crucial role since only a 
finite number of CH 3-manifolds with the same volume exist [^. The key facts are: 
the volumes of CH 3-manifolds obtained by Dehn surgeries on a cusped manifold 
Mc are always less than the volume of M^, CH 3-manifolds converge to Mc in the 
limit \p\, \q\ oo. As shown in figure p2.1| , the volume spectra are discrete but there 
are many accumulation points which correspond to the volumes of various cusped 
manifolds. The smallest cusped manifolds in the known manifolds have volume 
2.0299 which are labeled as "m003" and "m004"in SnapPea. m003 and m004 are 
topologically equivalent to the complement of a certain knot in the lens space L^^i 
and the complement of a "figure eight knot" (figure [2?^) , respectively ||3^. (3,-1) and 
(-2, 3) Dehn surgeries on m003 yield the smallest and the second smallest known 
manifolds, which are called the Weeks manifold (volume=0.9427) and the Thurston 



Figure 2.2: Cusped Hyperbolic Manifold 

The figure in the left side represents a "figure eight knot" and the figure 
in the right shows the Dirichlet domain of a cusped manifold m004 viewed 
from two opposite directions in the Klein(projective) coordinates where 
geodesies and planes are mapped into their Euclidean counterparts. The 
two vertices on the left and right edges of the polyhedron which are iden- 
tified by an element of the discrete isometry group correspond to a cusped 
point. Colors on the faces correspond to the identification maps. One ob- 
tains the Dirichlet domain of m003 by interchanging colors on quadrilateral 
faces in the lower(or upper) right figure. 
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manifold (volume=0.9814), respectively. It should be noted that the volume of CH 
manifolds must be larger than 0.16668... times cube of the curvature radius although 
no concrete examples of manifolds with such small volumes are known [p^] . As \p\ 
and |g| becomes large, the volumes converge to that of M^. Similarly, one can do 
Dehn surgeries on m004 or other cusped manifolds to obtain a different series of CH 
manifolds. 



Figure 2.3: Smallest CH Manifolds 

Dirichlet domains of the Weeks manifold (left) and the Thurston mani- 
fold(right) viewed from opposite directions in the Klein(projective) coor- 
dinates. Colors on the faces correspond to the identification maps. The 
manifolds can be obtained by performing (3,-1) and (-2,3) Dehn surgeries 
on m003, respectively. 
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Chapter 3 
Mode Functions 



Science grows slowly and gently; reaching the truth by a variety of errors. One 
must prepare the introduction of a new idea through long and diligent labour; 
then, at a given moment, it emerges as if compelled by a divine necessity ... 

(Karl Gustav Jacob Jacobi, 1804-1851) 

In locally isotropic and homogeneous background spaces each type (scalar, vector 
and tensor) of first-order perturbations can be decomposed into a decoupled set of 
equations. In order to solve the decomposed linearly perturbed Einstein equations, 
it is useful to expand the perturbations in terms of eigenmodes (=eigenfunctions that 
are square-integrable) of the Laplace-Beltrami operator which satisfies the Helmholtz 
equation with certain boundary conditions 

{V^ + k^)ukix) = 0, (3.1) 

since each eigenmode evolves independently. Then one can easily see that the time 
evolution of the perturbations in multiply-connected locally isotropic and homo- 
geneous spaces coincide with those in the locally and globally homogeneous and 
isotropic FRWL spaces whereas the global structure of the background space is de- 
termined by these eigenmodes. In the simplest closed flat toroidal spaces obtained 
by gluing the opposite faces of a cube by three translations, the mode functions are 
given by plane waves with discrete wave numbers. In general, the mode functions 
in locally isotropic and homogeneous spaces can be written as a linear combination 
of eigenf unctions on the universal covering space. 

The property of eigenmodes of CH spaces have not been unveiled until recently 
although a number of theorems concerning with lower bounds and upper bounds 
of eigenvalues in terms of diffeomorphism-invariant quantitities have been known in 
mathematical literature. However, despite the lack of knowledge about the prop- 
erty, eigenmodes of CH spaces have been numerically investigated by a number of 
authors in the field of number theory and harmonic analysis and of quantum chaol- 
ogy which are connected each other at a deep level. For the last 15 years various 
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numerical techniques have been apphed to the eigenvalue problems for solving the 
Helmholtz equation with periodic boundary conditions (manifold case), Neu- 
mann and Dirichlet boundary conditions (orbifold case). Eigenvalues of CH 2-spaces 



has been numerically obtained by many authors[^ ^ 0, |35|, ^ Eigenvalues 



of cusped arithmetic and cusped non-arithmetic 3-manifolds with finite volume have 



been obtained by Grunewald and Huntebrinker using a finite element methodjSS 



Aurich and Marklof have computed eigenvalues of a non-arithmetic 3-orbifold us- 
ing the direct boundary element method(DBEM) |p9|| . The author has succeeded in 
computing low-lying eigenmodes of the Thurston manifold, the second smallest one 
using the DBEM^^], and later the Weeks manifold, the smallest one in the known 
CH manifolds using the same method PT|. Cornish and Spergel have also succeeded 
in calculating eigenmodes of these manifolds and 10 other CH manifolds based on 
the Trefftz method p2|. 

In quantum mechanics, an eigenmode can be interpreted as a wave function of a 
free particle at a stationary state with energy E = k'^. The statistical property of the 
energy eigenvalues E and the eigenmodes m^; of the classically chaotic systems has 
been investigated for exploring the imprints of classical chaos in the corresponding 
quantum system (e.g. see|Q and other articles therein). Because any classical dy- 
namical systems of a free particle in CH spaces (known as the Hadamard-Gutzwiller 
models) are strongly chaotic (or more precisely they are K-systems with ergodic- 
ity, mixing and Bernoulli properties |31|]), it is natural to expect a high degree of 
complexity for each eigenstate. In fact, in many classically chaotic systems, it has 
been found that the short-range correlations observed in the energy states agree 
with the universal prediction of random matrix theory (RMT) for three universal- 
ity classes:the Gaussian orthogonal ensemble(GOE), the Gaussian unitary ensem- 



ble(GUE) and the Gaussian symplectic ensemble (GSE)[^ For the Hadamard- 
Gutzwiller models the statistical property of the eigenvalues and eigenmodes is de- 
scribed by GOE (which consist of real symmetric NxN matrices H which obey the 
Gaussian distribution oc exp (— Tr_f/'^/(4a^)) (where a is a constant) as the systems 
possess a time-reversal symmetry. RMT also predicts that the squared expansion 
coefficients of an eigenstate with respect to a generic basis are distributed as Gaus- 
sian random numbers |^ which had been numerically confirmed for some classically 
chaotic systems ( [p6i 46]). Although these results are for highly excited modes in 
semiclassical region, the low-lying modes may also retain the random property. For 
the analysis of CMB temperature fiuctuations, the statistical property of the expan- 
sion coefficients of low-lying eigenmodes is of crucial importance since it is relevant 
to the ensemble averaged temperature correlations on the present horizon scales. 

For classically chaotic systems, the semiclassical correspondence is given by the 
Gutzwiller trace formula]^ which relates a set of periodic orbits in the phase space 
to a set of energy eigenstates. Note that the correspondence is no longer one-to-one 
as in classically integrable systems. Interestingly, for CH spaces, the Gutzwiller trace 
formula gives an exact relation which had been known as the Selberg trace formula in 
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mathematical literature 1 48]. The trace formula gives an alternative method to com- 
pute the eigenvalues and eigenmodes in terms of periodic orbits. The poles of energy 
Green's function are generated as a result of interference of waves each one of which 
corresponds to a periodic orbit. Roughly speaking, periodic orbits with shorter 
length contribute to the deviation from the asymptotic eigenvalue distribution on 
larger energy scales. In fact, a set of zero-length orbits produce Weyl's asymptotic 
formula. Because periodic orbits can be obtained algebraically, the periodic orbit 
sum method enables one to compute low-lying eigenvalues for a large sample of man- 
ifolds or orbifolds systematically if each fundamental group is known beforehand. 
The method has been used to obtain eigenvalues of the Laplace-Beltrami operator 
on CH 2-spaces and a non-arithmetic 3-orbifold|]52|, |3^, However, it has not 
been applied to any CH 3-manifolds so far. 

In section 1 we first formulate the DBEM for computing eigenmodes of the 
Laplace-Beltrami operator of CH 3-manifolds. In section 2. we study the statistical 
property of low-lying eigenmodes on the Weeks and the Thurston maniolds, which 
are the smallest examples in the known CH 3-manifolds. In section 3 we describe the 
method for computing the length spectra and derive an explicit form for computing 
the spectral staircase in terms of length spectrum which is applied for CH manifolds 
with volume less than 3. In section 4 we analyze the relation between the low- ly- 
ing eigenvalues and several diffeomorphism-invariant geometric quantities, namely, 
volume, diameter and length of the shortest periodic geodesies. In the last section 
the deviation of low-lying eigenvalue spectrum from the asymptotic distribution is 
measured by C— function and the spectral distance. 



3.1 Boundary element method 
3.1.1 Formulation 

The boundary element methods (BEM) use free Green's function as the weighted 
function, and the Helmholtz equation is rewritten as an integral equation defined 
on the boundary using Green's theorem. Discretization of the boundary integral 
equation yields a system of linear equations. Since one needs the Discretization on 
only the boundary, BEM reduces the dimensionality of the problem by one which 
leads to economy in the numerical task. To locate an eig envalue, the DBEM [] 
requires one to compute many determinants of the corresponding boundary matrices 
which are dependent on the wavenumber k. 

Firstly, let us consider the Helmholtz equation with certain boundary conditions, 

(V' + fc2)M(x) = 0, (3.2) 



^The DBEM uses only boundary points in evaluating the integrand in Eq.(3.6). The indirect 
methods use internal points in evaluating the integrand in Eq.(p]^) as well as the boundary points. 
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which is defined on a bounded M-dimensional connected and simply-connected do- 
main Q which is a subspace of a M-dimensional Riemannian manifold Ai and the 



boundary is piecewise smooth. V = V*Vi 



1,2, ■ ■ ■,M), and Vj is the 



covariant derivative operator defined on A^. A function u in Sobolev space H^{Q) 
is the solution of the Helmholtz equation if and only if 



7^[M(x),^;(x)] = ((V' + P) m(x), t;(x) ) = 



(3.3) 



where v is an arbitrary function in Sobolev space H^{Q) called weighted function 
and ( ) is defined as 

{a,b) = I ab^dV. (3.4) 
Jn 

Next, we put m(x) into the form 



M 

U = ^Uj^j, 



(3.5) 



where (f)j^s are linearly independent square- integrable functions. Numerical methods 
such as the finite element methods try to minimize the residue function TZ for a fixed 
weighted function f (x) by changing the coefficients Uj. In these methods, one must 
resort to the variational principle to find the u/s which minimize TZ. 

Now we formulate the DBEM which is a version of BEMs. Here we search m(x)'s 
for the space C^{Q) n C^iQ) n L'^{Q). First, we slightly modify Eq.(|H3) using the 
Green's theorem 

{V^u)v^dV - [ iV^v)uy/gdV = I {V,iu)v^dS' - [ {V^v)u^dS\ (3.6) 
a Jn Jan Jan 

where g = det {gij} and dV = dxi . . . dxM'i the surface element dS^ is given by 

dSi : 



(3.7) 



where ^jj^-.-jM+i denotes the Mfl- dimensional Levi-Civita tensor. Then Eq. ( p.3|) 
becomes 

[ {V^v + k\)uy/^ dV + [ {Viu)v^dS'-f iViv)uy/^dS' = 0. (3.8) 
Jn Jan Jan 
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As the weighted function v, we choose the fundamental solution G£;(x, y) which 
satisfies 

iV' + E)GE{^,y) = Sni^-y), (3.9) 

where E = k'^, and ^^(x — y) is Dirac's delta function, ^^(x, y) is also known as 
the free Green's function whose boundary condition is given by 

lim GB(x,y) = 0, (3.10) 

d(x,y)^oo 

where (i(x, y) is the geodesic distance between x and y. Let y be an internal point 
of Q. Then we obtain from Eg. (^.81) and Eq. 



w(y)+/ GEi^,y)V,u^dS'~ [ i\/,GE{^,y))u^dS' = 0. (3.11) 
Jan Jan 

Thus the values of eigenfunctions at internal points can be computed using only 
the boundary integral. If y G dQ, we have to evaluate the limit of the boundary 
integral terms as G'£;(x, y) becomes divergent at x = y (see appendix A). The 
boundary integral equation is finally written as 

lu{y)+ [ GE{^,y)ViU^dS'- [ {V,GE{^,y))u^dS' = 0, (3.12) 
^ Jan Jan 

or in another form, 

\u{y)+! GE{^,y)^^n^V9dS- [ ^^^j^n^ u ^ dS = 0, (3.13) 
2 Jan ox' Jq^ dx' 

where n* = dS'/dS and dS = ^fdSUiSi. Note that we assumed that the boundary 
surface at y is sufficiently smooth. If the boundary is not smooth, one must calcu- 
late the internal solid angle at y (see appendix A). Another approach is to rewrite 
Eq. ( p.llj ) in a regularized form . We see from Eq. ( p.l2| ) or Eq. ( p.l3|) that the ap- 
proximated solutions can be obtained without resorting to the variational principle. 
Since it is virtually impossible to solve Eq.( |3.13| ) analytically, we discretize it using 
boundary elements. Let the number of the elements be N. We approximate u by 
some low-order polynomials (shape function) on each element as m = ci -|- C2 ?7 + C3 
where rj and ^ denote the coordinates on the corresponding standard element ^ 
Then we have the following equation: 

[H]{u} = {G]{q}, q^—, (3.14) 

where {u} and {q} are A^-dimensional vectors which consist of the boundary val- 
ues of an eigenfunction and its normal derivatives, respectively. [H] and [G] are 

^It can be proved that the approximated polynomial solutions converge to w(x) as the number 
of boundary elements becomes large [50, 
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N X A^-dimensional coefficient matrices which are obtained from integration of the 
fundamental solution G'£;(x, y) and its normal derivatives multiplied by Ui and g,, 
respectively. Note that the elements in [H] and [G] include k implicitly. Because 
Eq. (|3.14|) includes both u and q, the boundary element method can naturally incor- 
porate the periodic boundary conditions: 

m(x) = M(5(i(x)), g(x) = -g(5(i(x)), ondil, (3.15) 

where (^j's are the face-to-face identification maps defined on the boundary. The 
boundary conditions constrain the number of unknown constants to N. Application 
of the boundary condition ( p.l5| ) to Eg. ( p.l4|) and permutation of the columns of 



the components yields 

[A]{x} = 0, (3.16) 

where N x iV-dimensional matrix A is constructed from Gij and Hy and A^-dimensional 
vector X is constructed from Wj's and g^'s. For the presence of the non-trivial solution, 
the following relation must hold, 

det[A]=0. (3.17) 

Thus the eigenvalues of the Laplace-Beltrami operator acting on the space (Q) (1 
C'^{Q) n L'^{Q) are obtained by searching for /c's which satisfy Eq. (|3.17|) . 



3.1.2 Computation of low-lying eigenmodes 



In this section, we apply the DEEM for computing the low-lying eigenmodes of CH 
3-manifolds. The Helmholtz equation in the Poincare coordinates is written as 



1 - Ixl 



1 - Ixl 



X- V? 



u + k'^u = 0, 



(3.18) 



where A^; and V^; are the Laplacian and the gradient on the corresponding three- 
dimensional Euclidean space, respectively. Note that we have set the curvature 
radius R = 1 without loss of generality. By using the DEEM, the Helmholtz equation 
( |3.18 ) is converted to an integral representation on the boundary. Here Eq.( |3.13 ) 
can be written in terms of Euclidean quantities as 



Qu 

G'fc(x,y) — dS - 
an ox' 



9Gfc(x, y) 



wn^ dS 



0, 



idQ dx' 

where dS = 2(1 — |xp)^^ dSs- The fundamental solution is given as |3 



(3.19) 



Gk (x,y) 



4^ 7^^^ 



TT TT 

- - < arg s < -, 



(3.20) 
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Figure 3.1: 1168 Boundary Elements 



where s = — k"^ and cr = cosh(i(x, y). Then Eq.( |3.19| ) is discretized on the 
boundary elements Fj as 



1 ^ 
-n(x,) + J2 



Myj) f 9Gk{^i^yj'^..r.. _o, (3.21) 



j=i 



/r 



Gk{xi,yj)—^;^dS - j u{yj)dS 



where denotes the number of the boundary elements. An example of = 1168 
elements on the boundary of the fundamental domain in the Poincare coordinates 
is shown in figure p.l| . These elements are firstly generated in Klein coordinates in 
which the mesh-generation is convenient. The maximum length of the edge Al in 
these elements is 0.14. The condition that the corresponding de Broglie wavelength 
271 /k is longer than the four times of the interval of the boundary elements yields a 
rough estimate of the validity condition of the calculation as fc< 11. On each Fj, u 
and q = du/dn are approximated by low order polynomials. For simplicity, we use 
constant elements: 

m(xj) = = Const. , 5'(xj) = = Const. , on Fj. (3.22) 

Substituting Eq. (1X221) into Eq.(g), we obtain 



N N 

j=i j=i 

Hij I^J 



H 



ij 



Hij — 2 I = J, 



where 



Hij = ^{xi,yj)dS{yj), Gjj = Gk{^i,yj) dSiyj). (3.23) 

The singular integration must be carried out for I-I components as the fundamental 
solution diverges at (x/ = y/). This is not an intractable problem. Several numerical 
techniques have already been proposed by some authors [Q, ^]. We have applied 
Hayami's method to the evaluation of the singular integrals |5^. Introducing coor- 
dinates similar to spherical coordinates centered at x/, the singularity is canceled 
out by the Jacobian which makes the integral regular. 

Let gi {i = l,2, . . . ,8) be the generators of the discrete group F which identify a 
boundary face Fj with another boundary face gi{Fi): 

gi{xi) = Xi, Xi G Fi. (3.24) 
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The boundary of the fundamental domain can be divided into two regions OQa and 
OQb and each of them consists of N/2 boundary elements, 



on A = UFi, dVtB = Ugi{Fi), i 
The periodic boundary conditions 

u{gi{xi)) = M(xi), g(fi'i(xi)) = -g(xi), 



1,2, 



1,2 8 



(3.25) 



(3.26) 



reduce the number of the independent variables to N, i.e. for all G dfls, there 
exist Qi ET and x^ G OQa such that 



u{xb) = u{gi{xA)) = u{xa), g(xB) = -q{gi{xA)) 
Substituting the above relation into Eq. (|3.23|) , we obtain 



Haa + Hab 
Hba + Hbb 



-Gaa + Gab 
-Gba + Gbb 



Ua 



(3.27) 



(3.28) 



where ua = {u^, m^, . . . m^^^) and qa = {q^, (f-,--- q^^"^) and matrices H = {Hjj} and 
G = {Gij} are written as 



H 



Haa Hab 
Hba Hbb 



G 



Gaa Gab 
Gba Gbb 



(3.29) 



Eq. d^) takes the form 

[A{k)]{x} = 0, (3.30) 

where x A^-dimensional matrix A is constructed from G and H and A^-dimensional 
vector X is constructed from ua and qA- For the presence of the non-trivial solution, 
the following relation must hold. 



det[A{k)] = 0. 



(3.31) 



Thus the eigenvalues of the Laplace-Beltrami operator in a CH space are obtained 
by searching for fc's which satisfy Eq. (|3.31|) . In practice, Eq.( p.31 ) cannot be exactly 
satisfied as the test function which has a locally polynomial behavior is slightly de- 
viated from the exact eigenfunction. Instead, one must search for the local minima 
of det[A(k)]. This process needs long computation time as A{k) depends on k im- 
plicitly. Our numerical result (A;<13) is shown in table p.l| . 

The first "excited state" which corresponds to k = ki is important for the under- 
standing of CMB anisotropy. Our numerical result ki = 5.41 is consistent with the 
value 5.04 obtained from Weyl's asymptotic formula 



N[u] 



Vol{M)u^ 



= a/F - 1, Z/»l, 



(3.32) 
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Weeks 


Thurston 


k 


m{k) 


k 


m{k) 


k 


m{k) 


k 


m{k) 


5.268 


1 


11.283 


1 


5.404 


1 


10.686 


2 


5.737 


2 


11.515 


1 


5.783 


1 


10.737 


1 


6.563 


1 


11.726 


4 


6.807 


2 


10.830 


1 


7.717 


1 


12.031 


2 


6.880 


1 


11.103 


2 


8.162 


1 


12.222 


2 


7.118 


1 


11.402 


1 


8.207 


2 


12.648 


1 


7.686 


2 


11.710 


1 


8.335 


2 


12.789 


1 


8.294 


1 


11.728 


1 


9.187 


1 






8.591 


1 


11.824 


1 


9.514 


1 






8.726 


1 


12.012 


2 


9.687 


1 






9.246 


1 


12.230 


1 


9.881 


2 






9.262 


1 


12.500 


1 


10.335 


2 






9.754 


1 


12.654 


1 


10.452 


2 






9.904 


1 


12.795 


1 


10.804 


1 






9.984 


1 


12.806 


1 


10.857 


1 






10.358 


1 


12.897 


2 



Table 3.1: Eigenvalue k and Multiplicity m{k) 



assuming that no degeneracy occurs. One can interpret the first excited state as 
the mode that has the maximum de Broglie wavelength 27T/ki. Because of the 
periodic boundary conditions, the de Broglie wavelength can be approximated by the 
"average diameter" of the fundamental domain defined as a sum of the inradius r_ 
and the outradius which yields ki = 4.9 just KWo less than the numerical value. 
From these estimates, supercurvature modes in small CH spaces {Vol{/iA) ~ 1) are 
unlikely to be observed. 

To compute the value of eigenfunctions inside the fundamental domain, one 
needs to solve Eq.( p.30|) . The singular decomposition method is the most suitable 



numerical method for solving any linear equation with a singular matrix A, which 
can be decomposed as 

A = U^DV, (3.33) 

where U and V are unitary matrices and D is a diagonal matrix. If Da in D is almost 
zero then the complex conjugate of the i-th row in V is an approximated solution of 
Eq.( p.30| ). The number of the "almost zero" diagonal elements in D is equal to the 
multiplicity number. Substituting the values of the eigenfunctions and their normal 
derivatives on the boundary into Eq. (|3.11|) , the values of the eigenfunctions inside 



■^The inradius r_ is the radius of the largest simply-connected sphere in the fundamental domain, 
and the outradius r+ is the radius of the smallest sphere that can enclose the fundamental domain. 
r_ =0.535, r+ =0.7485 for the Thurston manifold. 
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Figure 3.2: Lowest Eigenmodes 

The lowest eigenmodes of the Weeks manifold (A; = 5.268) and the Thurston 
manifold {k = 5.404) on a slice z = in the Poincare coordinates. The 
boundaries of the copied Dirichlet domains {solid curves) are plotted in 
full curve. 



the fundamental domain can be computed. Adjusting the normalization factors, one 
would obtain real-valued eigenfunctions. Note that non-degenerated eigenfunctions 
must be always real-valued. 

The numerical accuracy of the obtained eigenvalues is roughly estimated as fol- 
lows. First, let us write the obtained numerical solution in terms of the exact solution 
as k = ko + dk and tifc(x) = Ukoi^x) + Suki^x.) , where ko and Mfco(x) are the exact eigen- 
value and eigenfunction, respectively. The singular decomposition method enables 
us to find the best approximated solution which satisfies 

[A]{x} = e, |e| « 1, (3.34) 

where e is a N-dimensional vector and | | denotes the Euclidean norm. It is expected 
that the better approximation gives the smaller |e|. Then Eq. ( p.34 ) can be written 



/ Gk,+sk{^,yj)i^ + ih + Skf){uko{x) + 6uk{x))y^dV^ = e{yj). (3.35) 
Jn 

Ignoring the terms in second order, Eq. ( p.35| ) is reduced to 

/ Gk{x,yj){{A + kl)Sukix) + 2k5kukix)),/^dV^ = e{yj). (3.36) 
Jn 

Since it is not unlikely that (A + /cQ)(5Mfc(x) is anticorrelated to 2k6kuk{x), we obtain 
the following relation by averaging over yj, 

2k\6k\(jj^Gk{x,yj)uk{x)^dV^\^ ~ < |e| >, (3.37) 

where <> denotes the averaging over yj Thus one can estimate the expected devi- 
ation of the calculated eigenvalue \6k\ from Mfe(x) and e(yj). We numerically find 
that \6k\ = 0.005 for k = 5.41 and \6k\ = 0.01 for k = 9.91(Thurston manifolds). 
The other deviation values lie in between 0.005 and 0.01. 

By computing the second derivatives, one can also estimate the accuracy of the 
computed eigenfunctions. The accuracy parameter err is defined as 

err{k, x) = (A + k^)ukix), (3.38) 



where Ufc(x) is normalized {0{uk{x)) ~ 1). We see from figure 



that the accuracy 
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Figure 3.3: Error Estimation 

Errors (err) in terms of hyperbolic distance d to the boundary from 291 
points in the fundamental domain of the Thurston manifold. 



becomes worse as the evaluation point approaches the boundary. However, for points 
with hyperbolic distance d > 0.1 between the evaluating point and the nearest 
boundary, the errors are very small indeed: err = 10~^ ^. This result is considered 
to be natural because the characteristic scale L of the boundary elements is ~ 0.07 
for our 1168 elements. If d < L, the integrands in Eq. (|3.11| ) become appreciable on 
the neighborhood of the nearest boundary point because the free Green's function 
approximately diverges on the point. In this case, the effect of the deviation from 
the exact eigenfunction is significant. If d » L, the integrand on all the boundary 
points contributes almost equally to the integration so that the local deviations are 
cancelled out. 

As we shall see in the next section, expansion coefficients are calculated using the 
values of eigenf unctions on a sphere. Since the number of evaluating points which 
are very close to the boundary is negligible on the sphere, expansion coefficients can 
be computed with relatively high accuracy. 



3.1.3 Statistical property of eigenmodes 

In pseudospherical coordinates {R,x,0,(f)), the eigenmodes are written in terms of 
complex expansion coefficients C,uim and eigenmodes on the universal covering space, 

u. = Y.UmXM)yiUOA), (3.39) 

Im 

where u = ^/k"^ — 1 and Yim is a (complex) spherical harmonic. The radial eigen- 
function Xi^i is written in terms of the associated Legendre function P, and the F 
function as 



Then the real expansion coefficients a^^im are given by 

O'uoo = ~^^{^uoo), Q-uio = -\/c^Re(^y/o), 
a^ijn = V2Re{^^im), m > 0, 

ttyim = -V2lm{^^i_rn), m < 0, (3.41) 

where 

2 



[l + Re{F{uJ))y 



F(z/«) F(/ — m + 1) 
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Figure 3.4: Random Nature of Expansion Coefficients 

Expansion coefficients a^im are plotted in ascending order as /(/ + 1) + m + 
1,0 < / < 20 for eigenmodes A: = 5.268(left) and /c = 12.789(right) on the 
Weeks manifold at a point which is randomly chosen. 



o-uim 's can be promptly obtained after the normalization and orthogonalization of 
these eigenmodes. The orthogonalization is achieved at the level of 10~^ to 10"'^ (for 
the inner product of the normalized eigenmodes) which implies that each eigenmode 
is computed with relatively high accuracy. 

In figure ^]4| one can see that the distribution of a^im'^ which are ordered as 
/(/ + 1) + m + 1 are qualitatively random. In order to estimate the randomness 
quantitatively, we consider a cumulative distribution of 



-'ulm 



\0'vlm 0'v\ 



(3.43) 



where is the mean of a^im^ and a1 is the variance. If Oi^^m's are Gaussian then h^im 
's obey a distribution P{x) = (l/2)^/^r(l/2)x~^/^e~^/^ with 1 degree of freedom. 
To test the goodness of fit between the the theoretical cumulative distribution I{x) 
and the empirical cumulative distribution function In{x), we use the Kolmogorov- 
Smirnov statistic Dj^j which is the least upper bound of all pointwise differences 
|J^(x)-/(x)| [H, 

Dn = sup \lNix) - I {x)\. (3.44) 

In{x) is defined as 

0, X < yi, 

In{x) = {j/N, y,<x<y,+i, j = 1, 2, . . . , iV-1, (3.45) 

1, Vn < X, 

where yi < y2 < ■ ■ ■ < yN are the computed values of a random sample which 
consists of N elements. For random variables Z^at for any 2; > 0, it can be shown 
that the probability of D]\f<d is given by 



lim P{Dn < d = zN-^/^) = L{z), (3.46) 

where 

00 

L{z) = 1 - 2^(-l)^-ie-2^'"'. (3.47) 
i=i 

From the observed maximum difference Dj^ = d, we obtain the significance level 
= 1 — P which is equal to the probability oi Dn> d. If is found to be large 
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Thurston 


k 


ao 


k 


ao 




5.404 


0.98 


10.686b 


7.9 X 10- 


-4 


5.783 


0.68 


10.737 


0.96 




6.807a 


0.52 


10.830 


0.67 




6.807b 


7.1 X 10-* 


11.103a 


0.041 




6.880 


1.00 


11.103b 


8.8 X 10- 


15 


7.118 


0.79 


11.402 


0.98 




7.686a 


0.26 


11.710 


0.92 




7.686b 


2.3 X 10-^ 


11.728 


0.93 




8.294 


0.45 


11.824 


0.31 




8.591 


0.91 


12.012a 


0.52 




8.726 


1.00 


12.012b 


0.73 




9.246 


0.28 


12.230 


0.032 




9.262 


0.85 


12.500 


0.27 




9.754 


0.39 


12.654 


0.88 




9.904 


0.99 


12.795 


0.76 




9.984 


0.20 


12.806 


0.42 




10.358 


0.40 


12.897a 


0.87 




10.686a 


0.76 


12.897b 


6.9 X 10- 


-4 



Table 3.2: Kolmogorov-Smirnov Test I 

Eigenvalues k and the corresponding significance levels an for the 
test of the hypothesis /Ar(x) ^ I{x) for the Thurston manifold are 
listed. The injectivity radius is maximal at the base point. 
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Figure 3.5: Symmetry 

A Dirichlet domain of the Thurston manifold in the Klein coordinates 
viewed from opposite directions at Q where the injectivity radius is lo- 
cally maximal. The Dirichlet domain has a Z2 symmetry (invariant by 
7r-rotation)at Q. 



enough, the hypothesis /Ar(x) ^I{x) is not verified. The significance levels for 
< / < 20 for eigenmodes < 13 on the Thurston manifold are shown in table |3.2| . 
The agreement with the RMT prediction is fairly good for most of eigenmodes which 
is consistent with the previous computation in However, for five degenerated 
modes, the non-Gaussian signatures are prominent (in PD[, two modes in (A; < 10) 



have been missed). Where does this non-Gaussianity come from? 

First of all, we must pay attention to the fact that the expansion coefficients 
o-uim depend on the observing point. In mathematical literature the point is called 
the base point. For a given base point, it is possible to construct a particular class of 
fundamental domain called the Dirichlet [fundamental) domain which is a convex 
polyhedron. A Dirichlet domain VL{x) centered at a base point x is defined as 

n{x) = Pi H{g, x) , H{g, x) = {z\d{z, x) < d{g{z),x)}, (3.48) 



where g is an element of a Kleinian group r(a discrete isometry group of PSL{2, C)) 
and d{z, x) is the proper distance between z and x. 

The shape of the Dirichlet domain depends on the base point but the volume 
is invariant. Although the base point can be chosen arbitrarily, it is a standard to 
choose a point Q where the injectivity radiusQ is locally maximal. More intuitively, 
Q is a center where one can put a largest connected ball on the manifold. If one 
chooses other point as the base point, the nearest copy of the base point can be 
much nearer. The reason to choose Q as a base point is that one can expect the 
corresponding Dirichlet domain to have many symmetries at Q ||58|| . 

As shown in figure |3.5| the Dirichlet domain at Q has a Z2 symmetry (invariant 



by TT-rotation) if all the congruent faces are identified. Generally, congruent faces 
are distinguished but it is found that these five modes have exactly the same values 
of eigenmodes on these congruent faces. Then one can no longer consider ayim^ as 
"independent" random numbers. Choosing the invariant axis by the vr- rotation as 
the 2;-axis, a^imS are zero for odd m's, leading to a non-Gaussian behavior. It should 
be noted that the observed Z2 symmetry is not the subgroup of the isometry group 
(or symmetry group in mathematical literature) D2 (dihedral group with order 2) 
of the Thurston manifold since the congruent faces must be actually distinguished 



^The injectivity radius Vinj at a point p is equal to the radius of the largest connected and 
simply connected ball centered at p which does not cross itself. 
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Weeks 


Thurston 


k 


< Od > 


k 


< ao > 


k 


< aD > 


k 


< an > 


5.268 


0.58 


10.452b 


0.62 


5.404 


0.63 


10.686b 


0.62 


5.737a 


0.61 


10.804 


0.63 


5.783 


0.61 


10.737 


0.62 


5.737b 


0.61 


10.857 


0.62 


6.807a 


0.62 


10.830 


0.63 


6.563 


0.62 


11.283 


0.57 


6.807b 


0.62 


11.103a 


0.59 


7.717 


0.59 


11.515 


0.61 


6.880 


0.63 


11.103b 


0.60 


8.162 


0.61 


11.726a 


0.63 


7.118 


0.61 


11.402 


0.61 


8.207a 


0.65 


11.726b 


0.59 


7.686a 


0.61 


11.710 


0.62 


8.207b 


0.61 


11.726c 


0.61 


7.686b 


0.63 


11.728 


0.64 


8.335a 


0.59 


11.726d 


0.61 


8.294 


0.60 


11.824 


0.62 


8,335b 


0.62 


12.031a 


0.60 


8.591 


0.60 


12.012a 


0.63 


9.187 


0.59 


12.031b 


0.60 


8.726 


0.60 


12.012b 


0.61 


9.514 


0.56 


12.222a 


0.61 


9.246 


0.60 


12.230 


0.60 


9.687 


0.61 


12.222b 


0.62 


9.262 


0.63 


12.500 


0.63 


9.881a 


0.61 


12.648 


0.59 


9.754 


0.62 


12.654 


0.62 


9,881b 


0.62 


12.789 


0.59 


9.904 


0.60 


12.795 


0.62 


10.335a 


0.63 






9.984 


0.60 


12.806 


0.62 


10.335b 


0.60 






10.358 


0.62 


12.897a 


0.62 


10.452a 


0.63 






10.686a 


0.60 


12.897b 


0.56 






Table 3.3 


: Kolmofforov 


-Smirnov 


Test II 







Eigenvalues k and corresponding averaged significance levels < a_D > 
based on 300 realizations of the base points for the test of the hy- 
pothesis /jv(ic) / I{x) for the Weeks and the Thurston manifolds. 



in the manifoIdQ. 

Thus the observed non-Gaussianity is caused by a particular choice of the base 
point. However, in general, the chance that we actually observe any symmetries 
(elements of the isometry group of the manifold or the finite sheeted cover of the 
manifold) is expected to be very low. Because a fixed point by an element of the 
isometric group is either a part of 1-dimensional line (for instance, an axis of a 
rotation) or an isolated point (for instance, a center of an antipodal map). 

In order to confirm that the chance is actuahy low, the KS statistics an of a^imS 
are computed at 300 base points which are randomly chosen. As shown in table 
( |3.3| ), the averaged significance levels < an > are remarkably consistent with the 
Gaussian prediction, la oi are found to be 0.26 to 0.30. 

Next, we apply the run test for testing the randomness of a,y/m's where each set of 

^The observed Z2 symmetry might be a "hidden symmetry" which is a symmetry of the finite 
sheeted cover of the manifold (which can be tessellated by the manifold) . 
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Weeks 


Thurston 


k 


< Or > 


k 


< ar > 


k 


< ar > 


k 


< Or > 


5.268 


0.51 


10.452b 


0.52 


5.404 


0.48 


10.686b 


0.51 


5.737a 


0.48 


10.804 


0.52 


5.783 


0.45 


10.737 


0.49 


5.737b 


0.45 


10.857 


0.53 


6.807a 


0.53 


10.830 


0.53 


6.563 


0.54 


11.283 


0.49 


6.807b 


0.50 


11.103a 


0.52 


7.717 


0.50 


11.515 


0.51 


6.880 


0.47 


11.103b 


0.53 


8.162 


0.54 


11.726a 


0.51 


7.118 


0.50 


11.402 


0.51 


8.207a 


0.52 


11.726b 


0.48 


7.686a 


0.49 


11.710 


0.51 


8.207b 


0.49 


11.726c 


0.49 


7.686b 


0.52 


11.728 


0.49 


8.335a 


0.53 


11.726d 


0.48 


8.294 


0.50 


11.824 


0.54 


8,335b 


0.50 


12.031a 


0.54 


8.591 


0.50 


12.012a 


0.51 


9.187 


0.53 


12.031b 


0.51 


8.726 


0.51 


12.012b 


0.49 


9.514 


0.55 


12.222a 


0.54 


9.246 


0.43 


12.230 


0.51 


9.687 


0.53 


12.222b 


0.50 


9.262 


0.50 


12.500 


0.48 


9.881a 


0.51 


12.648 


0.54 


9.754 


0.54 


12.654 


0.48 


9,881b 


0.51 


12.789 


0.48 


9.904 


0.52 


12.795 


0.50 


10.335a 


0.54 






9.984 


0.49 


12.806 


0.51 


10.335b 


0.51 






10.358 


0.53 


12.897a 


0.57 


10.452a 


0.53 






10.686a 


0.51 


12.897b 


0.55 



Table 3.4: Run Test 

Eigenvalues k and corresponding averaged significance levels < > 
for the test of the hypothesis that the auirnS are not random 
numbers for the Weeks and Thurston manifolds. a^'s at 300 
points which are randomly chosen are used for the computation. 



o-uim 's are ordered as Z(Z + l) + m + l (see Suppose that we have n observations 



of the random variable U which falls above the median and n observations of the 
random variable L which falls below the median. The combination of those variables 
into 2n observations placed in ascending order of magnitude yields 

UUU LL UU LLL U L UU LL, 

Each underlined group which consists of successive values of f/ or L is called run. 
The total number of run is called the run number. The run test is useful because 
the run number always obeys the Gaussian statistics in the limit n — > oo regardless 
of the type of the distribution function of the random variables. As shown in table 
37^ , averaged significance levels <ar> are very high (la is 0.25 to 0.31). Thus each 



set of a^/m's ordered as /(/ + 1) + m + 1 can be interpreted as a set of Gaussian 
pseudo-random numbers except for limited choices of the base point where one can 
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Figure 3.6: Kolomogorov-Smirnov Test III 

Plots of m-averaged significance levels aD{i^,l) based on 300 realizations 
for the Weeks and the Thurston manifolds (0 < / < 20 and k < 13). n 
denotes the index number which corresponds to an eigenmode Uk where 
the number of eigenmodes less than k is equal to n {k(n = l) is the lowest 
non-zero eigenvalue). The accompanying palettes show the correspondence 
between the level of the grey and the value. 



observe symmetries of eigenmodes. 

Up to now, we have considered / and m as the index numbers of a,^im at a fixed 
base point. However, for a fixed (/, m), the statistical property of a set of cty/m's at a 
number of different base points is also important since the temperature fluctuations 
must be averaged all over the places for spatially inhomogeneous models. From 



figure |3lB one can see the behavior of m-averaged significance levels 



aD[y,l)= 2^ ^ (3.49) 

•m=—l 

which are calculated based on 300 realizations of the base points. It should be noted 
that each Qyim at a particular base point is now considered to be "one realization" 
whereas a choice of / and m is considered to be "one realization" in the previous 
analysis. The good agreement with the RMT prediction for components / > 1 has 
been found. For components / = 1, the disagreement occurs for only several modes. 
However, the non-Gaussian behavior is distinct in / = components. What is the 
reason of the non-Gaussian behavior for I = 0? Let us estimate the values of the 
expansion coefficients for / = 0. In general, the complex expansion coefficients C,uim 
can be written as, 

UmiXo) = I «.(Xo, e, 0) YlUe, <f))dQ. (3.50) 

For / = 0, the equation becomes 

^^oo(xo) = -— 7=^ / u^ixo,0,^)dn. (3.51) 

2V2smi/xo J 

Taking the limit Xo ^ 0, one obtains, 

e„„ = (3.52) 



Thus 0^,00 can be written in terms of the value of the eigenmode at the base point. 
As shown in figure p.2| the lowest eigenmodes have only one "wave" on scale of 
the topological identification scale L (which will be defined later on) inside a single 
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Dirichlet domain which imphes that the random behavior within the domain may 
be not present. Therefore, for low-lying eigenmodes, one would generally expect 
non-Gaussianity in a set of 0,^00 's- However, for high-lying eigenmodes, this may 
not be the case since these modes have a number of "waves" on scale of L and they 
may change their values locally in an almost random fashion. 

The above argument cannot be applicable to a,yim 's for / 7^ where X,yi ap- 
proaches zero in the limit Xo ~^ while the integral term 

u,{Xo,0,(P)YiU9,(P)dn (3.53) 

also goes to zero because of the symmetric property of the spherical harmonics. 
Therefore ayim 's cannot be written in terms of the local value of the eigenmode for 
/ 7^ 0. For these modes, it is better to consider the opposite limit Xo 00. It is 
numerically found that the sphere with very large radius Xo intersects each copy of 
the Dirichlet domain almost randomly (the pulled back surface into a single Dirichlet 
domain chaotically fills up the domain). Then the values of the eigenmodes on the 
sphere with very large radius vary in an almost random fashion. For large Xoj we 
have 

XH(xo)oce-^>^«+'^(^'')\ (3.54) 

where 0(i/, /) describes the phase factor. Therefore, the order of the integrand in 
Eq. ( |3.50| ) is approximately e~^^° since Eq. ( p.50| ) does not depend on the choice 



of Xo- As the spherical harmonics do not have correlation with the eigenmode 
Uv{xoidi'P)i the integrand varies almost randomly for different choices of (/,m) or 
base points. Thus we conjecture that Gaussianity of a,y/m's have their origins in the 
chaotic property of the sphere with large radius in CH spaces. The property may 
be related to the classical chaos in geodesic flowsQ . 

Let us now consider the average and the variance of the expansion coefficients. 
As the eigenmodes have oscillatory features, it is natural to expect that the averages 
are equal to zero. In fact, the averages of < a^im > 's over < / < 20 and —l<m<—l 
and 300 realizations of base points for each z/-mode are numerically found to be 
0.006 ± 0.04 - 0.02(lcr) for the Weeks manifold, and 0.003 ± 0.04 - 0.02(la) for the 
Thurston manifold. Let us next consider the i/-dependence (/c-dependence) of the 
variances Var{auim)- In order to crudely estimate the i/-dependence, we need the 



angular size 69 of the characteristic length of the eigenmode Ui^ at Xo[0 



/c2(smh(2(xo + rave)) - smh(2(xo - rave)) - ^rave) ' 



®If one considers a great circle on a sphere with large radius, the length of the circle is very 
long except for rare cases in which the circle "comes back" before it wraps around in the universal 
covering space. Because the long geodesies in CH spaces chaotically (with no particular direction 
and position) wrap through the manifold, it is natural to assume that the great circles also have 
this chaotic property. 



3.2. PERIODIC ORBIT SUM METHOD 



31 



Figure 3.7: Variance of Expansion Coefficients 

Averaged squared a^^s {k< 13) based on 300 realizations of the base points 
for the Weeks and the Thurston manifold with ztlcr run-to-run variations. 
Uu is defined to be Var{auim) averaged over 0<Z<20 and —l<m<l. The 
best-fit curves for the Weeks and the Thurston manifolds are 21.0z^~2 and 
20.3i/~^, respectively. 



where Vol{M) denotes the volume of a manifold M and Vave is the averaged radius 
of the Dirichlet domain. There is an arbitrariness in the definition of Vave- Here we 
define Vave as the radius of a sphere with volume equivalent to the volume of the 
manifold, 

VoliM) = 7r(sinh(2r„„e) - 2r,„e), (3.56) 

which does not depend on the choice of the base point. Here we define the topolog- 
ical identification length L as L = 2rave- For the Weeks and the Thurston manifold, 
L = 1.19 and L = 1.20 respectively. From Eq. ( p.55|) , for large xo, one can ap- 



proximate Uiy{xo) ~ '^'uix'o) by choosing an appropriate radius x'o which satisfies 
z/~^ exp(— 2xo) = z/'~^ exp(— 2x'q). Averaging Eq. ( |3.50|) over Us and m's or the base 
points, for large Xoy one obtains. 



(IWH ~ f'm .lrnl'), (3.57) 

exp(-2x'o) 

which gives {\^uim\'^) ~ Thus the variance of aj^^^'s is proportional to z/~^. The 
numerical results for the two CH manifolds shown in figure p.7| clearly support the 
z/~^ dependence of the variance. 

As we have seen, the property of eigenmodes on general CH manifolds is sum- 
marized in the following conjecture: 

Conjecture: Except for the base points which are too close to any fixed points by 
symmetries, for a fixed v, a set of the expansion coefficients a^im over (/, m) 's can 
be regarded as Gaussian pseudo-random numbers. For a fixed {vim) {I > 0), the 
expansion coefficients at different base points that are randomly chosen can be also 
regarded as Gaussian pseudo-random numbers. In either case, the variance is pro- 
portional to z/~^ and the average is zero. 



3.2 Periodic orbit sum method 
3.2.1 Length spectra 

Computation of periodic orbits (geodesies) is of crucial importance for the semiclas- 
sical quantization of classically chaotic systems. However, in general, solving a large 
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number of periodic orbits often becomes an intractable problem since the number 
of periodic orbits grows exponentially with an increase in length. For CH manifolds 
periodic orbits can be calculated algebraically since each periodic orbit corresponds 
to a conjugacy class of hyperbolic or loxodromic elements of the discrete isometry 
group r. The conjugacy classes can be directly computed from generators which 
define the Dirichlet domain of the CH manifold 

Let gi, {i = 1, ■■■,N) be the generators and X be the identity. In general these 
generators are not independent. They obey a set of relations 



n 



9ii9i2---9in=^, (3-58) 



which describe the fundamental group of M. Since all the elements of F can be 
represented by certain products of generators, an element (? G F can be written 

9 = 9ii9i2 ■ ■ ■ 9i^, (3.59) 

which may be called a ^^word" . Using relations, each word can be shorten to a 
word with minimum length. Furthermore, all cyclic permutations of a product of 
generators belonging to the same conjugacy class can be eliminated. Thus conju- 
gacy classes of F can be computed by generating words with lowest possible length 
to some threshold length which are reduced by using either relations ( |3.58|) among 
the generators or cyclic permutations of the product. 

In practice, we introduce a cutoff length l^ut depending on the CPU power be- 
cause the number of periodic orbits grows exponentially in / which is a direct conse- 
quence of the exponential proliferation of tiles (copies of the fundamental domain) 
in tessellation. Although it is natural to expect a long length for a conjugacy class 
described in a word with many letters, there is no guarantee that all the the periodic 
orbits with length less than l^ut is actually computed or not for a certain threshold 
of length of the word. 

Suppose that each word as a transformation that acts on the Dirichlet funda- 
mental domain D. For example, D is transformed to D' = gD by an element g. 
We can consider gD's as tessellating tiles in the universal covering space. If the 
geodesic distance d between the center (basepoint) of D and that of gD is large, 
we can expect a long periodic orbit that corresponds to the conjugacy class of g. 
Tessellating tiles to sufficiently long distance d > Icut makes it possible to compute 
the complex primitive length spectra {Lj = Ij exp{i(f)j), \lj < Icut} where Ij is the 
real length of the periodic orbit of a conjugacy class with one winding number and 
(pj is the phase of the corresponding transformation. We also compute multiplicity 
number m{lj) which counts the number of orbits having the same lj and (pj. 

In general, the lower limit of the distance d for computing a complete set of 
length spectrum for a fixed Icut is not known but the following fact has been proved 



by Hodgson and Weeksj^Oj. In order to compute a length spectra of a CH 3-manifold 
(or 3-orbifold) with length less than /, it suffices to compute elements {g} satisfying 

d{x,gx) < 2 cosh" ^( cosh _R cosh/ /2), (3.60) 
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Figure 3.8: Classical Staircases 

The classical staircases A'^(/) for the Weeks manifold and the Thurston 
manifold are well consistent with the asymptotic distribution ( 3.61| ) for 
large /. 

where x is a basepoint and R is the spine radiusQ. Note that there is a unique geodesic 
which lies on an invariant axis for each hyperbolic or loxodromic element. Snap Pea 
can compute length spectra of CH spaces either by the "rigorous method" based 
on the inequality ( p.60| ) or the "quick and dirty" method by setting the tessellating 
radius d by hand. The former method has been used for manifolds with small volume 
(< 1.42), but the latter method {d = lcut + 0.5) has been also used for some manifolds 
with large volume (> 1.42) since the tessellating radius given by the former method 
is sometimes so large that the computation time becomes too long. The detailed 
algorithm is summarized in appendix B. 

The asymptotic behavior of the classical staircase A^(/) which counts the number 
of primitive periodic orbits with length equal to or less than / for CH 3-spaces can 
be written in terms of / and the topological entropy r ||6^ 



Nil) ~ Ei(r/) ~ E^^EM, / ^ oo. (3.61) 
tI 

The topological entropy for D-dimensional CH spaces is given hj t = D — 1. A 
larger topological entropy implies that the efficiency in computation of periodic orbit 
is much less for higher dimensional cases p9|. 



In figure ^]^, the computed classical staircases {Icut = 7.0) are compared with 
the asymptotic formula for the smallest (Weeks) manifold and the second smallest 
(Thurston) manifold. For both cases, an asymptotic behavior is already observed 
at / ~ 3.5. 

Although the asymptotic behavior of the classical staircase N{1) does not depend 
on the topology of the manifold, the multiplicity number m{l) does. In fact, it was 
Aurich and Steiner who firstly noticed that the locally averaged multiplicity number 

<m(/)>=— g{k), (A^ = total number of terms) (3.62) 

l-Al/2<li<l+Al/2 

grows exponentially < m{l) >~ e'/^// as / ^ oo for arithmetic 2-spaces (manifolds 
and orbifolds)[^ |6l|, since the length I of the periodic orbits are determined 
by algebraic integers in the form 2cosh(//2) = algebraic integeT^. The failure of 

'^Spine radius R is equal to the maximum over all the Dirichlet fundamental domain's edges 
of the minimum distance from the edge to the basepoint. Note that R is finite even for a cusped 
manifold with finite volume. 

^For a two-dimensional space, the classical staircase has an asymptotic form N{1) ~ exp(Z)/L 
On the other hand, the classical staircase for distinct periodic orbits has a form N{1) ~ exp(Z/2) 
for arithmetic systems. Because m{l)dN = dN we have m{l) ~ exp(^/2)/Z as Z ^ oo. 
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Figure 3.9: Locally Averaged Multiplicities 

Plots of locally averaged multiplicities (A/ = 0.2) with one-sigma errors and 
the fitting curves for the Weeks manifold m003(-3,l) which is arithmetic 
(star) and a non-arithmetic manifold m004(l,2) (diamond). The fitting 
curves (a exp(/)// for the former and exp(6/) / (cl) for the latter) are obtained 
by the least square method using data 3.0 < / < 7.0. 
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Table 3.5: Coefficients of Fitting Curves for Averaged Multiplicities 

Coefficients of the fitting curves which describe the average behavior of lo- 
cally averaged multiplicities <m{l) > for arithmetic (A) and non-arithmetic 
(N) 3-manifolds. G denotes the isometry group (symmetry group). Fitting 
parameters a, b and c are obtained by the least square method using data 
3.0 < I < 7.0. t For m007(3,l), the isometry group may be larger than D2. 



application of the random matrix theory to some CH spaces may be attributed to the 
arithmetic property. For non-arithmetic spaces, one expects that the multiplicities 
are determined by the symmetries (elements of the isometry group) of the space. 
However, in the case of a non-arithmetic 3-orbifold, it has been found that < m{l) > 
grows exponentially in the form e^^ /{cl) where b and c are ffiting parameters ||39||. 
This fact might implies that the symmetries of long periodic orbits are much larger 
than that of the space even in the case of non-arithmetic systems. 

For 3-manifolds, one can expect that the property of < m{l) > for 3-orbifold 
also holds. The locally averaged multiplicities (A/ = 0.2) for the smallest twelve 
examples which include seven arithmetic and five non-arithmetic 3-manifolds |6^ 



have been numerically computed using SnapPea. From figure one can see that 
the difference in the behavior of <m{l) > between the non-arithmetic manifold and 
the arithmetic one is manifest. As observed in 3-orbifolds, averaged multiplicities 
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behave as 

a — J—, (arithmetic), (3.63) 

exp bl . . /o ^ .\ 

— J — , (non-arithmetic) , (3.64) 

where a depends on the discrete isometry group while b and c are fitting parameters. 
From table I, one observes that arithmetic manifolds having a larger symmetry group 
have a larger value of a. The growth rates for non-arithmetic manifolds (6 ~ 0.56) are 
always less than that for arithmetic manifolds {b = 1) but nevertheless exponential. 



<m{l)> = 



3.2.2 Trace formula 

Gutzwiller's periodic orbit theory provides a semiclassical quantization rule for clas- 
sically chaotic systems. The theory is expressed in form of a semiclassical approx- 
imation {h — >■ 0) of the trace of the energy Green's operator (resolvent operator) 
Ge = (A + E)~^ in terms of the length of periodic orbits (geodesies) {Li} which is 
known as the Gutzwiller trace formula^^. For the dynamical system of a free mas- 



sive particle on a CH space known as the Hadamard-Gutzwiller model, the periodic 
orbits give the exact eigenvalues, and the relation is no longer semi-classical approx- 
imation. In mathematical literature, the trace formula is known as the Selberg trace 
formula [Q. In what follows we consider only orientable CH 3-manifolds(denoted 
as CH manifolds) (for general cases including orbifolds, see [§9|). The Selberg trace 
formula for a CH manifold M = H^/T (F is a discrete isometry group containing 
only hyperbolic or loxodromic elements) can be written as 

TriGE-GE') = -^(P-P) 

K9to) f expi-ipligr)) exp{-ip'l{g. 



E 



4(cosh/(5(^) - cos0(5f^)) V ip ¥ 

where = E — 1, v{M) denotes the volume of M, l{gr) is the (real) length of 
the periodic orbit of transformation gr E T. g^g is a transformation that gives the 
shortest length of the periodic orbit ({gro) which commutes with gr- (pigr) is the 
phase of the transformation g^. The sum in ( |3.65| ) extends over F conjugacy classes 

{9r}:={9r\9r = hgrh-\heT} (3.66) 

of hyperbolic (0 = 0) or loxodromic elements (0 7^ 0). However the periodic orbit 
sum in ( p.65| ) which is known as MaaB-Selberg series converges at only complex 
energy such that Im p < —1 and Im p' < — 1. In order to obtain real eigenvalues, 
one needs to multiply the trace by some suitable analytic "smoothing" function h{q) 
that satisfies: 



36 



CHAPTER 3. MODE FUNCTIONS 



{i):h{q) = h{-q); 

{u):h{q) = for 5 > as |g| oo; 

(iii):/i(g) is analytic in the strip |Im g| < 1 + e for e > 0. 

Multiplying (|3.65| ) by q h{q) /{tt i) and integrating it over q from — oo to oo, one 
obtains the general Selberg trace formula, 

^ — ' ZTT 

n=0 

+ Y-^ , „ ^1^""^ —--hil{gr)), (3.67) 

2 cosh /Q?, - cos 007. ^ ^ 

{9t} 

where 

1 r°° 

h{l) = — dq h{q) exp(-zg/), (3.68) 

and /i"(0) is the second derivative of h and p„ denotes a wavenumber of the cor- 
responding eigenmode. The sum in (|3.67|) is absolutely convergent for any real 
eigenvalues = + 1. One can obtain various functions of eigenvalues such as 
heat kernels and energy level densities from periodic orbits by choosing an appro- 
priate "smoothing" function h{q). 

In order to obtain eigenvalues, a simple approach is to compute the spectral 
staircase 

oo 

iV(E) = ^^(E-EO, (3.69) 

n=0 

where i?o = 0, Ei, E2, ■ ■ ■ are the eigenvalues of the Laplacian and 6 is the Heaviside 
function ||32||. To explore supercurvature modes ue, < E < 1, we choose the 
"smoothing" function of N{E) as 



V(p') = i(l-Erf(^))^ 



which is real for E > 0. Note that hp{p') satisfies all the conditions (i) to (iii). By 
taking the limit e ^ 0, one obtains the spectral staircase 

00 

N{E)=hmY,hUPn). (3.71) 

n=0 

Let us first estimate the behavior of the trace in ( |3.67| ). From a straightforward 
calculation, the zero-length contribution can be written as 

6 vra 



^ (r(3/4))~iF(5/4, 1/2, p^e') + |^r(3/4)F(7/4, 3/2, p^i^.j)} 
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where r(a;) is the Gamma function and F{a,b,z) is the confluent hypergeometric 
function. For x=p'^/e'^^oo, it is asymptotically expanded as 

+ g^e-x-^l V - - (x)-" + (3.73) 
r(a) 1^ n! I 

■ n=0 



where (a)^ = r(a + m)/r(a). From (|3.72| ) and (|3.73|) , in the lowest order, we have 
the average part, 



v{M)~„ v{M). 



, 3 



that gives the dominant term in the Weyl asymptotic formula for p>>l. 
Next, we estimate the oscillating term in ( p. 67 ) 

hp^^{l) = a j dq exp{-iql) q exp(^- ^^ ^ a = ^^2^3/2 • {'i-^^) 

In the long orbit-length limit I » 1 with p> 0, the integrand in Eg. ( p. 75]) oscillates so 
rapidly that the dominant contribution comes from or —p where (g^ —p^Y/e^ 
can be approximately given as — 4(g — p^p^ /e^. Then ( p.75|) can be written as 

V(0 ~ ^sin(p/)exp(^^^j. (3.76) 

Thus each periodic orbit corresponds to a wave with wavelength 2tt/1 and an ampli- 
tude which is exponentially suppressed with an increase in Z or a decrease in p. For 
a finite subset of length spectra Ij < Icut, the appropriate choice for the smoothing 
scale is given by e = ap^^"^ for > where a depends on l^ut since for a reasonable 
value of the proportional factor a all the contributions from periodic orbits with 
large length I > Icut can be negligible. For < 0, an optimal choice can be ob- 
tained from a numerical computation of Eq.( p.75| ) directly. Comparing the obtained 
smoothed spectral staircase with the one based on the computed "true" eigenvalues 
using the direct boundary element method (DBEM), it is numerically found that 
for /cut = 7.0, an appropriate smoothing scale is given by 

/0.116A;2 + 0.184 fc + 1.2 (A; < 1) , . 

'^^^ ~ \ 0.832A;i/2 + 0.668 (A; > 1) ' ^"^'^^^ 

where E = k'^ = p"^ + 1. The eigenvalues can be computed by searching E at which 
N{E) — 0.5 becomes positive integer. It should be emphasized that the precision of 
computation depends on the value of Icut which determines the resolution scale in 
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Figure 3.10: Spectral Staircases 
The spectral staircases N{k) for the Weeks manifold and the Thurston man- 
ifold obtained by the DBEM are compared with the average parts N(k) 
(solid curves), namely, Weyl's asymptotic formula ( |3.74| ) and smoothed 
spectral staircase (dotted curves) obtained by the periodic orbit sum 
method(POSM) using all periodic orbits / < 7.0. 
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Table 3.6: First Eigenvalues 

1/2 

The first (non-zero) wavenumbers ki = are calculated using the 

DBEM and using the POSM for seven smallest manifolds, kis agree with 
relative accuracy AA:i//ci = 0.01 — 0.08. The multiplicity number of the 
first eigenmode is calculated using the DBEM. 



the eigenvalue spectra. 

In order to get eigenvalues from the smoothed spectral staircase, one must take 
into account the effect of the multiplicity number(degeneracy number) for each eigen- 
value since the spectral staircase is smoothed on larger scales for degenerated modes. 
Therefore, the numerical accuracy becomes worse if the eigenmode has a large mul- 
tiplicity number. Fortunately, the order of the symmetry group is not so large for 
a small manifold (volume< 3). For instance, of the twelve smallest examples, nine 
manifolds have a symmetry group with order 4. If one assumes that the multiplicity 
number is either 1 or 2 then the deviation A/c from a precise value is approximately 
given by 

i(i_Erf(^^^(^))^i±M. ,,r8) 

For instance, if one uses a length spectrum / < 7.0 then ( |3.78|) gives Afc = 0.30, 0.33 
for fc = 5.0,3.0, respectively. If one permits the multiplicity number as much as 6, 
then the expected precision becomes A/c = 0.49, 0.54 for /c = 5.0, 3.0, respectively. 

We can see from table II that the first eigenvalues calculated by using length 
spectra / < 7.0 for some smallest known CH manifolds lie within several per-cent 
of those obtained by the DBEM. Note that the eigenvalues are also consistent with 
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those obtained by the Trefftz method P2 
non-zero mode is degenerated 
owing to the smoothing effect. From figure |3.10 



For two examples in which the first 
the eigenvalues are much shifted to lower values 
one can see that the curves of 



the obtained smoothed spectral stairs cross the "true" stairs at almost half height. 
A slight deviation from the average part of the spectral staircase is caused by the 
interference of waves each one of which corresponds to a periodic orbit. 



3.2.3 First eigenvalue and geometrical quantity 

The estimate of the first (non-zero) eigenvalue Ei = kf of the Laplace-Beltrami oper- 
ator plays a critical role in describing the global topology and geometry of manifolds. 
A number of estimates of Ei for n-dimensional compact Riemannian manifolds M 
using diffeomorphism-invariant quantities have been proved in mathematical litera- 
ture. 

First of all, we consider the relation between the first eigenvalue Ei and the 
diameter d which is defined as the maximum of the minimum geodesic distance 
between two arbitrary points on M. Various analytic upper and lower bounds of Ei 
in terms of d have been known. Suppose M with Ricci curvature bounded below by 
—L{L > 0). Cheng and Zhou proved that Ei satisfies 



El > max 



Ivr^ 1 TT^ L2 3 



2 



12 <P - 4^- + 16 - Ji. ^exp(-C„v/I^/2)J . (3.79) 



where C„ = ma,x\y/n — 1, v^]||64||. Another lower bound has been obtained by 
Lu[^. Suppose that the Ricci curvature of M is bounded below as Rat > —Kgab, {K > 
0) for a some real number K where Qab is the metric tensors of M. Then Ei satisfies 

rvr^ ^ 8 8 / d'^K\ 8 /d r— -\ / d r— -\- 

i!.>max[--A;---,-exp(-— ),-(l + -Vii-(n-l))exp(--V/f(„-l)) 

(3.80) 

As for upper bounds, the following theorem has been proved by Cheng [|6^] . Suppose 
M with Ricci curvature larger than (n — l)c, then we have 

Ei<Ei{Vnic,d/2)) (3.81) 

where Vn{c,r) denotes a geodesic ball with radius r in the n-dimensional simply- 
connected space with sectional curvature c and Ei is the first Dirichlet eigenvalue. 
Setting L = K = 2,n = 3, and c = — 1, we obtain the upper and lower bounds of Ei 
for CH 3-manifolds. For the upper bound, (|3.81|) gives a simple relation, d < 2tx/vi, 



where uf = kf — 1. The physical interpretation is clear: the wavelength Ai = 27r/z/i 
of the lowest non-zero mode must be larger than the diameter. 

Now, we compare the first eigenvalues of 263 examples of CH 3-manifolds with 
volume less than 3 which have the length of the shortest periodic orbit Imm > 0.3 (the 
Hodgson- Weeks census[^]) and of 45 other examples obtained by Dehn surgeries 
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Figure 3.11: Diameter vs First Eigenvalue 

Diameter d versus ki for 263 examples of CH 3-manifolds with Imin > 0.3 
and v<3 (group A, circle) and 52 examples that are obtained by performing 
Dehn surgeries (51 examples for \p\, \q\ < 10 plus 1 example {p, q) = (16, 13) 
group B, diamond) on a cusped manifold m003 with the best-fit curves for 
group B corresponding to ( 3.82| ) (dashed curve) assuming /ci(cusp) = 1 and 



( 3.83| ) (solid curve) assuming A;i(cusp) = 0.1. The fitting curves also agree 



with the computed values for group A. The upper and lower thick curves 
denote the analytic bounds (^1^), (|3^ ) and ( |3^ ). 

(IpI < 17, |g| < 14) on a cusped manifold mOO30 with the analytic bounds. The 
diameter of a CH 3-manifold is given by the supremum of the outradius(3 over all 
the basepoints, which has been numerically computed using the SnapPea kernelQ. 
The numerical accuracy is typically Ad = 0.03 — 0.09 depending on the topology of 
the manifold. 



As shown in figure p.ll| , the eigenvalues are well described by an empirical fitting 
formula Ai = jSd, or 

Applying the least square method for the 263 manifolds in the Hodgson-Weeks 
census (group A), and 53 manifolds obtained by Dehn surgeries on m003 (group B), 
the best-fit values P = 1.70, 1.73 have been obtained for each group, respectively. 
Note that /5 ~ 1.7 is slightly larger than the values 1.3-1.6 for 12 examples in the 
previous result by Cornish and Spergel]!^. The deviation from the fitting formula 
( |3.82[ ) is found to be remarkably small (with one sigma error Aki = 0.18,0.19 for 



each group, respectively), which implies the existence of much sharper bounds. 
The empirical formula ( |3.82| ) asserts that no supercurvature modes {i.e. ki<l) 



exist in the limit d ^ oo where the manifold converges to the original cusped 
manifold. However, cusped manifolds may have some supercurvature modes even 
for those with small volume. Instead of (p.82| ), we consider a generalized empirical 
formula 

h = J{h{cnsp)f + ^, (3.83) 

where A;i(cusp) < 1 is the smallest wavenumber for the original cusped manifold. 
Although no supercurvature modes were observed in this analysis, the non-existence 



^The Hodgson- Weeks census with volume less than 3 also includes 8 manifolds obtained by 
Dehn surgeries on m003. 

^"The outradius at a basepoint x is equal to the minimum radius of the simply-connected ball 
which encloses the Dirichlet domain at x. 

^^I would like to thank J. Weeks for providing me a code to compute the diameter using the 
SnapPea kernel. 
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Figure 3.12: Diameter vs Volume 
Diameter d versus volume v for 79 manifolds (/mm. > 0.5) (star), 184 man- 
ifolds (0.5 > Imin > 0.3) (circle) and 53 manifolds obtained by performing 
Dehn surgeries on m003 (group B, diamond) with the best-fit curves (3^ 



(dashed-dotted curve) and ( |3^ ) (dashed curve). The solid curve denotes 
the volume of m003. 



of such modes in the limit d ^ oo was not confirmed since the numerical accuracy 
becomes worse for manifolds with small ki and large d. 

Next, we consider the relation between diameter d and volume v of the manifold. 
Since diameter is given by the supreme of the outradius (minimum radius of a sphere 
which circumscribes the Dirichlet domain) all over the basepoints, one expect that 
V is estimated as the volume of a sphere (in a hyperbolic space) with radius r 
somewhat smaller than d if there is no region which resembles the neighborhood of 
a cusp ( "thin part" 0) or equivalently Imin is sufficiently large. Suppose that r = ad 
with a < 1 then the volume of a sphere 



V = 7r(sinh(2arf) - 2ad) (3.84) 

gives the approximate value of a CH manifold with diameter d. The best fit value 
for a sample of 79 manifolds with Zmj„>0.5 is a = 0.69. If a manifold has a "thin" 
part then the relation (|3.84|) is no longer valid since the diameter becomes too long. 
Let Vc be the volume of a cusped manifold (with only one cusp) and v{d) be 
the volume of a CH manifold M obtained by a Dehn surgery on Mc. In the limit 
d ^ CO, one can show that the following approximation holds (see appendix C): 

(a exp(-2(rf- rfo))\ ocrN 

v[d) = v^[l ^— ^ j , (3.85) 

where 6 denotes a ratio of the volume of the complementary part to that of the "thin" 
part and do is the diameter of the complementary part. For a sample of 41 manifolds 
with diameter longer than 1.65 in group B, the best fits are 5 = 0.0 and do = 0.25. 
As shown in figure 3.12| , the volume-diameter relation for CH manifolds with large 



Imin is well described by the fitting formula ( |3.84| ). As Imin becomes smaller, or 
equivalently, d becomes larger, a CH manifold M converges to the original cusped 
manifold Mc in which ( |3.85| ) holds. 

Finally, we look into the relation between the first eigenvalue and the volume. 
For CH manifolds with sufficiently large Imin, ( ^-831 ) and ( ^.84[ ) give 



47ra 

v{h)=7r{smh{g{h))-g{h)), g{h) = (3.86) 

PV(^i) - (A;i(cusp))^ 



12 



A "thin" part is defined as a region where the injectivity radius is short. 
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Figure 3.13: Volume vs First Eigenvalue 

Volume V versus ki for 79 manifolds with Imin > 0.5 (star), 184 mani- 
folds with 0.5 > Imin > 0.3 (circle) and 53 manifolds obtained by per- 
forming Dehn surgeries on m003 (group B, diamond) with fitting curves 
( |3.86| ) (solid curve), ( |3.87D (dashed-dotted curve), and ( |3. 88]) (dashed curve) 



where A;i(cusp) = 1 is assumed. 

Figure 3.14: Length of Shortest Periodic Orbit vs ki or Diameter 

The length of the shortest periodic orbit Imin versus the deviation Aki 
from the fitting formula ( 3.86| ) where /i;i(cusp) = 1 is assumed(left) and 



^min versus diameter d (right) for 263 CH manifolds (group A). 
To be consistent with the Weyl's asymptotic formula which is valid for ki » 1 



Mv) = \l{ — y + h (3.87) 



the fitting parameters should satisfy a/P = 3 ■ 4~^/^7r"^/^ ^ 0.44 which well agree 
with the numerically computed values a/ [3 = 0.69/(1.6 — 1.7) = 0.41 — 0.43 provided 
that /ci(cusp) = 1. One can see from figure |3.13| that both ( p.86| ) and ( |3.87|) give a 



good estimate of the first eigenvalue k\ for globally "slightly anisotropic" manifolds 
with t> < 3 and Imin > 0.5. For manifolds with large d, v and small Imin, ( |3.86| ) 
and ( p.87| ) give incorrect estimates for ki. In that case, the effect of the curvature 



cannot be negligible. In contrast to compact flat spaces, the volume of a sphere in 
hyperbolic space increases exponentially as the radius increases. Assuming that the 
relation ( p. 86] ) holds, for sufficiently globally "isotropic" CH manifolds (large Imin), 
ki is significantly larger than that for compact flat spaces with the same volume 
even if one assumes that ki{cusp) ~ 0. However, for CH manifolds converging to 
the original cusped manifold M^, the formula ( |3.86| ) has to be modified. If Imin is 
sufficiently small and d is large while keeping the volume finite then M has a "thin" 
part similar to the neighborhood of a cusp. Then one can use an asymptotic formula 
i ^M) instead of (jpl). 



where Vc is the volume of M^. As shown in figure |3.14| , ki shifts to a smaller value for 
manifolds with smaller Imin which have larger d. In the limit M —>■ M^, the length 
of periodic orbits of M also converges to that of except for the shortest orbit 
whose length Imin goes to zero. From the general Selberg trace formula ( |3.67| ), one 
can see that the wave which corresponds to the shortest orbit has a large amplitude 
with a long wavelength ~ 2t\: /Imin- Therefore, the presence of a very 
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Figure 3.15: Deviation from Weyl's Formula 

Spectral staircases for m003(9,7) (dashed-dotted staircase, left) , 
m003(16,13) (solid staircase, right) and m003(16,13) without the shortest 
periodic orbit (dotted staircase, right) are shown in comparison with the 
corresponding Weyl's asymptotic formula (solid curves) which have been 
numerically computed by the POSM using all periodic orbits with length 
/ < 7. N{k) = k perfectly fits the staircase of m003(16,13) for k < 5 (thick 
line, right). 



short periodic orbit results in the deviation of the energy spectrum at low-energy 
region (small k) from the asymptotic distribution. To confirm this, the spectral 
staircase using the length spectra of m003(16,13) but the shortest periodic orbit is 
removed has been computed. Note that m003(16,13) is very similar to the original 
cusped manifold m003 having d = 2.75 and Imin = 0.0086. As shown in figure 
p.l5| , the computed spectrum staircase agrees well with Weyl's asymptotic formula. 
The computed spectrum may coincide with the one for a manifold in which the 
"thin" part is cut off. On the other hand, the spectral staircases for m003(9,7) 
and m003(16,13) which have small deviate from the average part (=Weyl's 
asymptotic formula( p.74 )) (figure |3.15 ). These manifolds have a "thin" part which 



is virtually one-dimensional object. Let us remind that Weyl's asymptotic formula 
for a n-dimensional compact manifold M is given by 

m-^^^, (3.89) 

where a;„ = 2tc"'^'^ / (nT{n/2)) is the volume of the unit disk in the Euclidean n- 
space and v{M) is the volume of M. If one assumes that the Weyl formula still 
holds for small k, then the spectrum in the low-energy region which corresponds to 
fluctuations on large scales for these manifolds can be approximately described by 
the Weyl formula for n = 1, N{k) oc k/n. For m003(16,13), it is numerically found 
that N{k) = k provides a good fit for k < 5. In the next section, we will measure the 
deviation from the asymptotic distribution using the low-lying eigenvalue spectra. 



3.3 Spectral measurements of global anisotropy 

Among many possibilities, we should choose physically well-motivated quantities for 
measuring the global "anisotropy" in geometry in terms of eigenvalue spectra. First, 
we consider (^-function which is relevant to the CMB anisotropy. The angular power 
spectra for CH universes are approximately written as 



(3.90) 
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where Fi{k) can be approximated as a polynomial function of k for a given /. In 
order to measure the "anisotropy" , we define the following parameter in terms of 
the ^-function, 

A(s) = as) /Us) = E A;r^y ^ A:-f , k'^^ ~ ^ (3.91) 

i=l i=l * 

where kf are the eigenvalues of the Laplace-Beltrami operator on a CH 3-manifold M 
with volume v and fc^j are the eigenvalues obeying Weyl's asymptotic formula with 
volume V ( |3.87|) . Note that the zero- mode fco = is not included in the summation. 
Here we only consider the case s > 1 which ensures the convergence of the sum. 
The numerical result shows the clear difference between the 'slightly anisotropic" 
manifold m003(-3,l) {d = 0.82, Z^m = 0.58) , "somewhat anisotropic" manifold 
m003(9,7) {d = 2.11, /™„ = 0.028) and "very anisotropic" manifold m003(16,13) 
{d = 2.75, /mm = 0.0086) (figure |3.16|) . The presence of the fluctuations on large 
scales in these "anisotropic" manifolds shifts the corresponding (^-function to a larger 
value. For the case in which the shortest periodic orbit is removed from the length 
spectra of m003(16,13), one can see that the spectrum coincides with that obeying 
Weyl's asymptotic formula. 

Next, we consider the spectral distance ds proposed by Seriu which measures 
the degree of semi-classical quantum decoherence between two universes having one 
massless scalar field ||6^, 



..|M,Afl4f 1„1(| + |), (3.92) 

i=l * 

where kf and kf are the eigenvalues of the Laplace-Beltrami operator on a com- 
pact n- manifold M and M, respectively. Here we choose eigenvalues fc^j as kf. In 
practice we introduce a cutoff in the summation. It is numerically found that the 
summation converges rapidly for the 3 examples, namely, m003(-3,l), m003(9,7) and 
m003(16,13). The contribution of the first several terms dominates the summation 
(figure [3?T6| ) . The result implies that a universe having a spatial geometry m003(9,7) 
or m003(16,13) semiclassically decoheres with a universe having a spatial geometry 
m003(-3,l) (figure [3T7|). 
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Figure 3.16: Spectrum Measurements 

The figure in the left shows A(s)=(^(s)/C«)(s) for 3 examples of CH 3- 
manifolds ,m003(9,7) (dashed-dotted curve), m003(16,13) (solid curve), 
m003(-3,l) (dashed curve) and m003(16,13) where the shortest periodic 
orbit is removed (dotted curve). The figure in the right shows ds{n) = 
h X]r=i \ {ki/kwi + kwi/ki) for the same examples where n is the cutoff 
number in the summation. 



Figure 3.17: Anisotropy in Geometry 

Plots of a Dirichlet domain of the Weeks manifold m003(-3,l) (left) 
and that of m003(16,13) (right) viewed from two opposite directions in 
the Klein coordinates. The Dirichlet domain of m003(16,13) is quite 
similar to that of the original cusped manifold m003 (see figure 2.2). 
m003(16,13) has a "thin" part which is similar to the neighborhood of 
a cusp. The colors on the faces correspond to the identification maps. 
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Chapter 4 

CMB Anisotropy 



The most incomprehensible thing about the universe is that it is comprehensi- 
ble. 

(Albert Einstein, 1879-1955) 

The surprising discovery of the cosmic microwave background (CMB) by Penzias 
and Wilson in 1964 provided firm evidence that our universe started from a hot 
big-bang. The CMB is considered as the fossil of the photons as long ago as redshift 
2; ~ 10^ when the scale factor was 10^^ of the present value and the temperature 
was hot as T = 3 X lO^K. The almost perfect blackbody spectrum in the CMB 
implies that the radiation and matter were once in good thermal contact because 
of the interaction between photons and electrons through Compton scattering when 
the universe was very hot. As the temperature decreased gradually owing to the 
cosmic expansion, the electrons have started to combine with protons. Eventually, 
the density of free electrons became too low to maintain the thermal contact and 
photons decoupled with matter. During the subsequent cosmic expansion the wave- 
lengths of the decoupled photons are stretched which leads to the low temperature 
2.7 K in the background today. 

The COBE satellite launched in 1989 discovered a slight imperfection in the 
CMB:the temperature anisotropy at the level of one part in 10~^ on angular scale 
~ 10° and roughly agree with the scale-invariant (n = 1) initial spectrum [Q. The 
temperature anisotropy provides not only the information about the physical con- 
ditions at redshifts z ~ 10'^ but also provides an important probe of the spatial 
geometry on the present horizon scales. On superhorizon scales which are much 
larger than the Jeans scale, the interaction between radiation and matter is negligi- 
ble. Therefore the photon-baryon fluid can be seen as a single fluid. The anisotropy 
is determined by only gravitational effects. There are two kinds of effects which 
account for the generation of the CMB anisotropy |]7T|, p^ . One is the so-called or- 
dinary Sachs- Wolfe (OSW) effect. When the photons climbs out the gravitational 
potential at the last scattering, the temperature of the photons are redshifted. Be- 
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Figure 4.1: COBE-DMR Anisotropy Sky Map 

Temperature anisotropy detected by the COME-DMR experiment is at the 
level of AT/T ~ 10"^ The plotted sky map of temperature fluctuations 
AT (in the Aitoff projection) has been obtained by averaging data at fre- 
quencies 31 GHz, 53GHz, and 90GHz The monopole and the dipole 
components have been removed. 



cause the intrinsic temperature fluctuations at the last scattering is specifled by the 
gravitational potential if the initial condition is given, the temperature anisotropy 
is associated with the spatial fluctuation of the gravitational potential at the last 
scattering. Another one is the so-called integrated Sachs- Wolfe (ISW) effect. If the 
gravitational potential decays while the photons are climbing out the potential, the 
expected redshift would be reduced. The decay occurs when the fluctuation scale 
enters the sound horizon while the radiation is still dominant and when the curva- 
ture term or the cosmological constant dominates the total energy. 

It is of crucial importance to test whether the CMB anisotropy in spatially mul- 
tiply connected models is consistent with the COBE data. The observable effect of 
the non-trivial spatial topology becomes much prominent for models with smaller 
volume since apparent fluctuations beyond the size of the fundamental domain are 
strongly suppressed. In other words the periodical boundary conditions restrict the 
wave numbers k to have discrete values that are equal to or larger than the that of 
the flrst non-zero mode ki (the flrst eigenvalue). As we have seen, ki is related with 
the diameter, the shortest length of periodic geodesies and volume of the space. 

If the CMB anisotropy is completely determined by the OSW effect then the 
"mode-cutoff" in the eigenmodes implies the suppression on the large-angle power 
spectrum. Fluctuations on large angular scales which corresponds to fluctuations on 
large physical scales beyond the actual size of the space are strongly suppressed. In 
fact, for closed flat toroidal models without the cosmological constant, the COBE 
data constrain the size of the fundamental domain L > 0.8R^, where is the co- 
moving radius of the last scattering surface and L is the side of the cube []TU|, 0, |T2|, 
131 |74| . It should be emphasized that the constraint itself is not enough stringent 
to "rule out" closed flat toroidal models since it corresponds to < 8 where A^ 
is the number of copies of the fundamental domain within the observable region at 
present. 

In contrast, for low matter density models, the constraint could be considerably 
milder than the locally isotropic and homogeneous flat (Einstein-de-Sitter) models 
since a bulk of large-angle CMB fluctuations can be produced by the ISW effect 
which is the gravitational blueshift effect of the free streaming photons by the de- 
cay of the gravitational potential |1^, |75]. As the gravitational potential decays in 



either A-dominant epoch or curvature dominant epoch, the free streaming photons 
with large wavelength (the light travel time across the wavelength is greater than or 
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comparable to the decay time) that chmbed a potential well at the last scattering 
experience blueshifts due to the contraction of the comoving space along the trajec- 
tories of the photons. Because the angular sizes of the fluctuations produced at late 
time are large, the suppression of the fluctuations on scale larger than the topolog- 
ical identification scale does not lead to a significant suppression of the large-angle 
power if the ISW effect is dominant. In other words, the relevant fluctuations which 
have been produced at late time 2; = 1 ~ 3 are "small" enough compared with the 
size of the fundamental domain that they are not affected by the "finiteness" in the 
size of the spatial hypersurface. In fact, recent works |TB|, |7S|, |7S| have shown 
that the large-angle power (2 < / < 20) is completely consistent with the COBE 
DMR data for closed hyperbolic (CH) models whose spatial geometry is described 
by a small CH orbifold, the Weeks and the Thurston CH manifolds with volume 
0.36, 0.94 and 0.98 in unit of the cube of the curvature radius, respectively. Note 
that the Weeks manifold is the smallest and the Thurston manifolds is the second 
smallest in the known CH manifolds. 

These results are clearly at odds with the previous constraints [ PU| , [HT[ | on CH 



models based on pixel-pixel correlation statistics. It has been claimed that the sta- 
tistical analysis using only the power spectrum is not sufficient since it can describe 
only isotropic (statistically spherically symmetric) correlations. This is true inas- 
much one considers fluctuations observed at a particular point. Because any CH 
manifolds are globally anisotropic, expected fluctuations would be also statistically 
globally anisotropic at a particular point. In order to constrain CH models, it is 
necessary to compare the expected fluctuation patterns at all the place with the data 
since CH manifolds are globally inhomogeneous. However, it should be emphasized 
that the stringent constraints obtained in the previous analyses [^, |81| are only for 
CH models at a particular observation point Q where the injectivity radius is locally 
maximum for 24 particular orientations. Q is rather a special point at which some 
of the mode functions have a symmetric structure. It is often the case that the base 
point Q becomes a fixed point of symmetries of the Dirichlet fundamental domain 
or the manifold. As we have seen, the mode functions for closed hyperbolic spaces 
are well described by a set of Gaussian numbers which corresponds to a choice of the 
position and orientation. Therefore, we expect that the likelihood should be highly 
dependent on the choice of the place and orientation of the observer which had been 
ignored in the previous analysis. 

In this chapter we first derive the evolution of the linear perturbations expressed 
in the gauge-invariant formalism. Next, we solve the Boltzmann equation on the 
perturbed background for computing the CMB anisotropy. Then we investigate the 
angular power spectra for estimating the suppression on the large angular scales. 
Finally we carry out the Bayesian analysis using the COBE DMR 4-year data incor- 
porating the effects of inhomogeneity and anisotropy of the background geometry 
for testing closed flat and hyperbolic models with low matter density. 
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4.1 Linear perturbation theory 

The observable physical quantity should not be dependent on choice of gauge. How- 
ever, it is also true that a different time slice gives a different value of fluctuation. 
In other words, observers with different dynamics would observe different universe. 
If the physical state of the observer is fixed then no freedom is allowed. For instance 
the synchronous gauge is a natural choice for the present free falling observer. If 
the observer were not free faUing but at the rest frame of the total matter then the 
total matter gauge is a useful choice. Since the choice of the physical state of the 
imaginary observer in the past is arbitrary, physically, none of gauge choices has 
priority over others. Mathematically speaking, the gauge freedom arises due to the 
freedom of correspondence between a perturbed manifold M and a reference (locally 
homogeneous and isotropic) manifold M. 

However, it is useful to choose a gauge in which we can "easily" understand the 
evolution of perturbation. Newtonian gauge is useful in analyzing propagation of 
photons while the description of matter becomes simpler in the total matter gauge. 
All we have to do is to choose a completely fixed gauge in which physical interpre- 
tation and evolution of the fluctuations are simplest. 

The line element of the perturbed FRWL metric takes the form 

ds'^ = -0^(77) [(1 + 2A)dr]^ + Bjdr]dx^ - {-fij + 2HLiij + 2HTij)dx'dx^] , (4.1) 
where 

^ j _ f -K'^{dx^ + sinh'^xidO'^ + si'n?9d(f)'^) K<0 , . 

A,Bj,HL,a.nd HTijSnc functions of space and time that describe a perturbation and 
Hxij denotes a three-dimensional traceless tensor. It is convenient to decompose 
the metric perturbation variables in the following way. First, a 3- vector Bj can be 
decomposed as the sum of the flrst derivative of a longitudinal 3-scalar Bl and a 
transverse 3- vector Bxj- 

vector Bj = VjBl + Btj , VjBt^ = 0, (4.3) 

where Vj is deflned as the three dimensional covariant derivative with metric 7ij. 
Similarly, a 3-tcnsor Hij can be decomposed as the sum of a longitudinal 3-scalar H^, 
the second derivative of a traccfrcc 3-scalar Ht, the first derivative of a transverse 
3-vector Hlj and a transverse tracefree 3-tensor Httij as 

tensor Hij = jijHL + ^iVj - jij —J Ht + V^Htj + V^Hti + HrTij ■ (4.4) 



For locally isotropic and homogeneous background geometry, Einstein equations are 
divided into a decoupled set of equations each one of which consists of only one type 
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of perturbation such as scalar, vector or tensor type perturbation^ Consequently, we 
can argue separately for each type of components. 
The infinitesimal gauge transformation 

7] fj = rj + T{x), —>■ X-' = + L\x), (4-5) 

transforms the metric as 

= g^u + 9''u{5x^'),a-g^^{5x^),u-g^.,x5x\ (4.6) 

which gives the following transformation law, 

scalar type A = A — T' T, 

a 

Bl = Bl + L'l~T, 

6 a 

Ht = Ht-L, (4.7) 

vector type B^j = Btj + L'j.- , 

Hrj = Htj-Ltj, (4.8) 

tensor type HxTjk = HrTjk, (4.9) 
where ' denotes d/drj. 

The perturbed energy momentum tensor is written as 
T\ = {P+~p)u°u, + P5l 

Tj = {P + p)u% + P6'j + 6PT'j, (4.10) 



where 



P = P + 6Pl, p = p + 6p, = u^ ^ Su", 

3 

= (a-\ 0,0,0), = (-a, 0,0,0), ^5P^\ = 0. (4.11) 

i=l 



^If the perturbed spatial geometry is locally homogeneous but anisotropic, we have coupled 
equations between different types of perturbations so that the decomposition is no longer convenient 
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P is the averaged pressure, p is the averaged energy density, 6Pl is the longitudinal 
pressure perturbation and j the anisotropic pressure of the fluid. Note that 
we have chosen coordinates for which the averaged velocity of the fluid vanishes. 
Let us introduce perturbation variables for the density and the stress, 

TTL = TTt j = — d = — , (4.12) 

and for the three-dimensional velocity, 

V^^J-^^^-f- — ■ (4-13) 
dt dr] u" ^ ^ 

Since the three dimensional velocity is defined on the space with metric 7^^, the 
index of is raised or lowered with 7'-' , 7jj, respectively. 

In terms of these variables, the energy momentum tensor is written as 

T\ = n = (p+p)(y,-s,), 

T\ = -{p + p)V\ T^^^p(S\ + 7rLS^k + 7^Jk)- (4-14) 

Like the metric perturbation, the perturbation variables for matter are decomposed 
into different types of components. The transformation law is described as 

scalar type 5 — 5 + 3(1 + w)—T, 

a 



Al + w) a', 



TTL = TTz, + 3Cs T, 

w a 

TTt — TTt, (4-15) 
vector type Vrj — Vrj + L'rpj, 

Tttj = TTTj, (4.16) 

tensor type TfrTjk = T^TTjk, (4-17) 
where w = P/p, and = dP/dp. 

For simplicity, we will consider only scalar type perturbations. It should be 
noted that the degree of apparent freedom of the metric perturbation is four and 
the degree of gauge freedom is two, so that the degree of actual freedom of the 
metric perturbation is two for scalar type perturbations. 
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In linear theory, each normal mode evolves independently. Therefore, mode 
expansion is a very useful method in analyzing the evolution of perturbed quanti- 
ties. From now on we follow the notations by Kodama and Sasaki |j7^. A harmonic 
(=a mode function of the Laplacian) is defined as a function which satisfies the 
Helmholtz equation 

iV^ + e)uk = 0, (4.18) 

The line element for the scalar type perturbation can be written in terms of har- 
monics as 



ds = —a 



[1 + 2Au)dr]'^ + BiUj dr] dx^ - 



{{l + 2HLu)jij + 2HTUij)dx'dx^ , (4.19) 



where Ui and defined as 



V 



2 



Uj 



~k-'V,u, Ui, = k'\ViVj - —iij)u. (4.20) 



Note that we have omitted the labels for the wave number k in Eq.( 4.1£ ) since there 
is no mode coupling. The energy momentum tensor T^^, is given by 

T\ = -p(l + H, T", = {p + P){Vl-Bl)u,, 

T\ = -(p + P)Vlu^, Th = P{5'k + 'nLu5^k + T^Tuh). (4.21) 

We have to make a choice of gauges in which the physical interpretation and 
evolution are simplest. The Newtonian gauge, defined as Bi = Ht = giving a 
time slicing in which the expansion is isotropic, is convenient gauge in analyzing 
the propagation of photons. On the other hand, the evolution of matter and metric 
becomes simplest in the total matter gauge, defined as Vr = -Bl, Ht = 0, giving a 
time slicing in which the total matter is at rest. We choose a hybrid choice of repre- 
sentation for fluctuations. For the metric fluctuation, we choose the perturbations 
defined in the Newtonian gauge, namely the Newtonian curvature $ and potential 
\Ef, which are written in terms of quantities in an arbitrary gauge as 



$ = n-k-^[^yg, (4.22) 



^ = A-k^^l^^^g-k-^a'g, (4.23) 

where 

n = HL + \ht, ag = k-'H^ - Bl. (4.24) 
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Note that in the Newtonian gauge, $ = Hl, and \E' = A. 
For matter perturbations, we choose the four variables as 

V ^ V,-^, (4.25) 

A = S + 3{l + w)-{VL-BL)k-\ (4.26) 

r = TTL-^S, (4.27) 
w 

n = ttt, (4.28) 

where V is the velocity of the fluid in the gauge with Ht = 0, A is the density 
perturbation in the gauge with Bl = Vl where the total matter is at rest, and 11 is 
the anisotropic pressure of the fluid. F can be considered as the entropy perturbation 
because it can be written as 

SP- cl5p 4 1 - 3w 

F = ^— — = s 

P 3 l + w ' 



5{s^/n) 



(4.29) 



where is the entropy density of photons, and n is the number density of non- 
relativistic particles. In other words, F is proportional to the fluctuation in the ratio 
of photon number density to the matter density since the entropy density of the 
photon is proportional to the photon number density. 

Now, we are ready to write the perturbed Einstein equations. First, the evolution 
of background is described as 

G\ (^)' + ^ = ^Z"'. ("-SO 

G\ ^ -K = 87^Gf>a^ (4.31) 



a V a 

where G^jy denotes the Einstein tensor. From the perturbed Einstein tensor 60^,^, 
we obtain the generalized Poisson equations, 

5G"o, 5G°i {e-3K)<^ = AnGa^pA, (4.32) 

5G\ A;2(^ + $) = -STrGa^Pn. (4.33) 

Prom the energy conservation T"^.^ = 0, we obtain the continuity equation 

A' - 3w-A = -(1 - ^)(1 + w)kV - 2(1 - ^)-wU, (4.34) 
and from the momentum conservation T'^^.^ — 0, we have the Euler equation, 

V' + -V^ -^kA + k<lf + T^kT - 1(1 - ^)-^kU. (4.35) 
a 1 + w 1 + w 3^ k^ + w ^ ^ 
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We derive the evolution equation for the Newtonian curvature $ assuming that 
the anisotropic pressure is neghgible 11 = 0. In this case, the Newtonian potential 
coincides with the Newtonian curvature except for the sign $ = — \E'. Let us consider 
the perturbed Einstein equation for dG\, 

- (^^)$ - $' = 47rGa2(P + p)A;~V, (4.36) 

Note that this can be also derived from the background equations and the continuity 
and the generalized Poisson equations. From the derivative of ( |4.36|) we have 

AirGa^iP + p)k-\V' + HV) = -(7^$)' - $" + 7^^$ + $'7^ 

+ {n<i> + <^'){P' + p')iP + p)-\ (4.37) 

where Ti = a' / a. On the other hand, the Euler and the Poisson equations yield 

kwT kcl (k'^-3K)^ ,^ , 

V' + n'V = - + ^ . ^ 9 - 4.38 

1 + w 1 + w ATcGa^p ^ ' 

and the background equations give 

47rGa2(P + p) = + + (4.39) 

and 

(J = -m[P^p). (4.40) 

From ([4.37|) , ([4.38|) , (|4.39| ) and (|4.4CI| ,) we finally obtain the second order ordinary 
differential equation for the Newtonian curvature 

$" + m{\ + c^)$' + clk^^ + (2H' + (1 + "icWn^ - K))^ = -AirGa'^PT. (4.41) 

Let us consider the adiabatic case (F = 0) and assume that the energy consists of 
two components: radiation and non-relativistic matter for which the sound speed 
is written in terms of energy density as = 4pr/(12pr- + 9pm)- First we consider 
the early universe when the radiation is dominant (P = p/3) For 77 << 1, the effect 
of the curvature is negligible. The background equations yield the evolution of the 
scale factor as a = a^r/ where is a constant. Then ( [4.41 ) is written as 

4 

$" + + — $ = 0, (4.42) 
1] 3 

which has an analytical solution 

$(77) = ri~'^{{uri cos^uT]) — sin(ci;?7))Ci + {ur] sin^urj) — cos{ujri))C2}, (4.43) 



4.2. SACHS -WOLFE EFFECT 



55 



where co = k/y/S and Ci and C2 are constants which depend on k. Around i] = 0, 
( |4.43| ) can be expanded as 



$(7/) = r]-^{C2 + ^a;V - - ^^V + (^(M)')}- (4.44) 

zoo 

Thus for the long- wavelength modes(ti;?7 << 1), the non-decaying mode (C2 = 0) 
has constant amplitude and the time derivative at the initial time vanishes, 

$ ~ const, $'(0) = 0. (4.45) 

During the radiation-matter equality time the equation of state change from w = 1/3 
to w = 0. Then the amplitude of the long-wavelength mode changes by a factor 
of 9/10. Using the explicit form for the scale factor for hyperbolic models in the 
matter dominant epoch, 

^" (cosh(y^r/) - 1), (4.46) 



2(1 -n 



the non-decaying mode of $ in the matter dominant epoch for hyperbolic models is 
given by 

^ . , ^, , 5(sinh^ r? — 3?7 sinh ?7 -|- 4 cosh ?7 — 4) 



where fj = \J—Kri and r/j„j denotes the initial time normalized by the present cur- 
vature. In the curvature or A dominant epoch the amplitude decays as 1/a. 

Using the Newtonian curvature $, the total density fluctuation can be readily 
obtained from the Poisson equation ( [4.32| ). In the radiation dominant and the mat- 



ter dominant epoch A evolves as but the amplitude freezes in the curvature or 
A dominant epoch. As for the velocity of the total matter, the continuity equation 
( |4.34|) gives a solution 



V = a-\7])C + a-\r]) j dr] a(r/) (^^-^kA - A;$) , (4.48) 

where C is a constant. The decaying mode corresponds to a~^{ri)C which can 
be omitted for considering the non-decaying mode. In the radiation and matter 
dominant epochs, the amplitude of V increases as a and as a}/"^ respectively but it 
decays in the curvature or A dominant epoch. 

4.2 Sachs -Wolfe effect 

The fundamental equation of radiative transport is derived from the Boltzmann 
equation which is written as 

|=C[/], (4.49) 



56 



CHAPTER 4. CMB ANISOTROPY 



where / = f{'K,p,t) is the distribution function and C [/] denotes the colhsion 
terms. We assume that the photons propagate freely after the decouphng {i.e. no 
reionization is considered), and hence we can neglect the collision terms (this is 
a good approximation for superhorizon perturbations even before the decoupling 
time) , 

df^dldx^ dld^^ (4.50) 
dt ^x^' dt dpt" dt ^ ^ 

which is also known as the Liouville equation. 

Although the phase space density of photons is conserved along their trajectories, 
each of the configuration space density and the momentum space density changes 
with time and space because the gravitational force affects the photon trajectories. 
As we shall see later, the redshift of photons can be interpreted as the change in the 
momentum space density which is affected by the background metric, gradients in 
the perturbative potential, and dilation effects caused by the curvature perturbation. 

First of all, we need to fix the gauge in order to solve the geodesic equation 
coupled with the Liouville equation on the perturbed FRWL space-time. We choose 
the Newtonian gauge in which the physical interpretation of the dynamics of photons 
is simplest. The line element takes the form, 

ds'^ = -(1 + 2*(x, t))dt'^ + a^(l + 2$(x, t))-fijdx'dx^ . (4.51) 

The geodesic equation for the photon is 

+ r^:,-7^-7^ = 0, P'^^-T^, (4.52) 



dX^ dX d\ ' ^ rfA ' 

where A is an afiine parameter. Since the photon momentum satisfies the relation, 

dx'' 

p° ~ dt' 

the geodesic equation can also be written as 



(4.53) 



dp" ^.u{}^dgap_dg^\p^ , . 

^ \2dx- dx^ ) v° ' ^ ' 



dt V2 dx'' dxf^ ) p 

Now, the Liouville equation is expressed as 



df df dx' , dfdp df df 



dt dx' dt dp dt (?7* dt 

where 



p2 = (p«)2^ y = a(l + $)-. (4.56) 

P 

Note that 7' is a three-dimensional vector defined on the space with metric 7^-. 
Since static curvature effects are unimportant in determining the redshift, it will be 
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assumed that the background metric is flat without loss of generahty {i.e. jij = Sij). 
The geodesic equation ([4 .541) with scalar type perturbations gives 



— — = \ (\& - 1) 2 — r — . (4.57) 

p dt dt dt a dt dx'^ a 

With the definition of p, we obtain the photon energy equation 

Idp _ Ual 9$ d^f\ 
p dt \dt a dt dx^ a J' 



The first term in the right-hand side of Eq.( [4.58|) is interpreted as the redshift by 



the isotropic expansion of the universe. The second term represents the dilation 
effect caused by the fluctuation of the spatial curvature. The third term can be 
considered as the gravitational redshift caused by the gradient in the potential. It 
is known that the CMB spectrum is well described by a blackbody spectrum, which 
obeys the distribution 

f{p) = {27r)-%e^/^ - (4.59) 

Then the fractional shift in frequency from the gravitational effect does not depend 
on frequency. In other words, the blackbody distribution will not change its form. 



f'ip'. 



1 1 



(27r)3 e(p'-^p)/^) - 1' 
11 

(27r)3ef'/T'_i' ^^-6°^ 



where 

T' = T + 6T = T{l + 6p/p). (4.61) 

Because the energy density of photons satisfies oc in the blackbody radiation, 
the temperature fluctuation 6 is defined as 

e = ^ = (4.62) 
T 4p, ^ ^ 

Note that p^ is a global averaged quantity (background value) which is averaged 
over the ensemble and 7 while p^ = p^ + is a quantity which is locally averaged 
over the ensemble, 



P7 



2 J f{p) Anp" dp, 

^,j-^^dp, T = T{r^,^,^). (4.63) 
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Substituting the total time derivative of p in the geodesic equation ( 4.54 ) into 
Eg. ( ^.581 ), and integrating over p, we obtain in the zero-th order, 



at 



-4Hp. 



11 



{H = a/a) 



(4.64) 



which represents the free-streaming of the photons Py oc a ^, and in the first order, 



A 2 dQ o 



The latter equation can be written as 

d 



9$ 9^ , 



(4.65) 



(9^ (9$ 

dr]^" ~^ '^^ dr] dr] ' 



(4.66) 



where dt = adrj. Integrating over the time, the temperature fluctuation is finally 
written as 



e(?7o,x„,7j + ^'(r/o,x,) = e(r/„x„7j + ^(r7„x,) + y i^— - —jdr], (4.67) 

where rj^: is the decoupling time, and rjo is the present time. The contribution of the 
intrinsic temperature fluctuation and potential redshift [0+^^] [r]^, x^,, 7^,] is called the 
ordinary Sachs- Wolfe (OSW) effect. The contribution of the remaining integrated 
term is called integrated Sachs- Wolfe (ISW) effect |7T|. For adiabatic perturbations, 
when the curvature perturbation begins to decay before the present time, the ISW 
effect dominates over the OSW effect. This is the case if the density parameter Qo 
is much smaller than the unity . 

Now, we estimate the temperature fluctuations on superhorizon scales at the de- 
coupling time. Using Eq.( 4.62| ), the temperature fluctuation can be written in terms 
of fc-component of the Newtonian density fluctuation of radiation (gauge invariant 
form)[| 

Ok = ^. (4.68) 



7 
4 



This can also be written in terms of quantities in total matter gauge as 

^ A, a'Vr 

- -j 

4 a k 



(4.69) 



where is the photon density fluctuation in the total matter gauge. Let us assume 
n = and the effect of the radiation at the last scattering is negligible. Substituting 



^More rigorously, 0k must be understood as the monopole {I = 0) component (Ok)o of the 
density fluctuations. On superhorizon scales, the contribution of higher multipoles of density 
fluctuations at the decoupling time is negligible in estimating the present anisotropics [[73[ . 
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the explicit form of the adiabatic perturbation ( 4.47|) for hyperbohc models into the 
Euler equation ( [4. 351 ), the non-decaying mode of the velocity of the total matter is 
given by 



V{v) 



5k^{ri^) 2r] + 7] cosh rj — 3 sinh rj 



(cosh 7] — 1] 



(4.70) 



where rj^, is the (normalized) last scattering time. On the other hand, from the 
continuity equation ( |4.34| ) the energy density A is explicitly written as 



1 3K 

fiv) = ^2\°^^*"'^2/ ^'^^^ 2 ~^ ^^^^ 



5 
3 



(4.71) 



From (1^, (FfO|) and (|TTD, we have 



ek(^)= (^-0(77^))$(77*), r/«l. 



(4.72) 



Up to the first order in rj, we obtain 



0k(r/*) ^ -$(r?*) 



(4.73) 



Because r]^ ~ 10~^ for the typical range of cosmological parameters, ( |4.73| ) gives a 
good approximation. Substituting Eq.( [4.73|) into Eq.( [4.67|) , we obtain the explicit 
formula for the temperature fluctuation. 



1 n° (9$ 

Qk{7]o, Xo, 7o) + ^(r/o, Xo) = --$(r7,x,) - 2 / —dr]. (4.74) 

3 J^, or] 



It should be noted that one should employ the full form of equation ( 4.67 ) instead of 
(|4.74|) for much rigorous calculation including the effect of the radiation contribution 
at the last scattering and the anisotropic pressure. However, for superhorizon scale 



temperature perturbation, ( [4.74|) gives a fairly good approximation provided that 
the matter dominates the energy at the last scattering time. 
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4.3 Angular power spectrum 

In order to analyze the angular dependence of the temperature anisotropy in the 
sky it is convenient to expand the anisotropy in terms of spherical harmonics Yim as 

AT 

— (n) = ^ai„l^„(n), (4.75) 

Im 

where n is the unit vector along the line of sight. Then the temperature 2-point 
correlation is given by 

AT AT \ 

^n) — (n'))= M,^>)YMY;^,{n'), (4.76) 



T ' ' T 



LmJ'm' 



where <> denotes the ensemble average taken over the initial conditions and the ob- 
servation positions and orientations. If the temperature fluctuations are described as 
a Gaussian random field, then no mode correlation occurs. Therefore, the statistics 
can be characterized by the angular power spectrum Cim defined as 

{O'lmOqim.i) = Cim5ll'5mm' ^ Clm = (k/m^)- (4.77) 

If the Gaussian random field is isotropic then the angular power spectrum has no 
dependency on m and can be written as Ci. From observation Ci is best estimated 
by simply averaging over m. 

The fact that we have only a finite number of 2/ + 1 samples for estimating Ci is called 
cosmic variance ^ which sets a limit on how precisely we can measure the 
angular power spectrum. Q obeys a distribution with 2/ + 1 degree of freedom. 
The fractional uncertainty in the estimate of Ci is given by 

^ ,/;^. (4.78) 

Ci V2/ + 1 ^ ^ 

For example, there is a 63% uncertainty in the quadrapole C2, 31% in Cio, and 7.1% 
in C200- 



4.3.1 Closed flat models 

Let us now consider cosmological models where the spatial geometry is represented 
as a flat 3-torus which is obtained by gluing the opposite faces of a cube by three 
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Figure 4.2: Suppression in Large-angle Power for Toroidal Models 

A large-angle suppression in AT; = ^Jl{l + 1)Ci/{2it) (all the plotted val- 
ues arc normalized by AT20 with infinite volume) occurs for the "stan- 
dard" toroidal model ($7^7 ^^a) = (1-0,0) at I < Icut ~ lixR^jL — 1 while 
such a prominent suppression is not observed for the toroidal model with 
(p^rai J7a)=(0. 1,0.9). The transfer function r;(A;) describes how each fc-mode 
contributes to the power for a particular angular scale I (lower figures). The 
corresponding first eigenvalues k\ are plotted as horizontal lines. The unit 
of k is B.Q. 



translations. Then the wave numbers of the square-integrable eigenmodes of the 
Laplacian are restricted to the discrete values hi — 27rni/L, {% — 1, 2,3) where n^'s 
run over all integers. Because the background geometry is globally anisotropic, the 
temperature anisotropy is no longer statistically isotropic. Therefore, the effect of 
the non-diagonal elements I ^ V oy m ^ m' m the correlation is non-negligible in 
general. The 5*0 (3) invariant quantity in the 2-point correlation consists of only 
diagonal elements as 

1 ^ 

m=—l 

where <> denotes an ensemble average taken over the initial conditions. For brevity, 
we also call Ci as the angular power spectrum. Q's describe the 2-point correlations 
for an isotropic ensemble since the non-diagonal elements in the 2-point correlation 
vanish if they are averaged over 50(3) transformations. Note that Q's are not 
independent each other since the fluctuations form a non-Gaussian random fleld 
for an isotropic ensemble. Thus C;'s still gives us information about the 2-point 
correlations although they are non-Gaussian and have a larger cosmic variance. 

Let us choose the Euclidean coordinates such that the z-axis is perpendicular to 
one of the face and the origin of the coordinates is located at the center of cube. 
Then the angular power spectrum is written as 



2 



87r-^P,(/,-)i^, 



kl 



1 r° 

-^t{ri*)3i{k{rio-ri.)H2 drj -^ji{k{rjo-v)), (4.81) 
3 Jr?, dV 



where Vi^{k) is the initial power spectrum for $ , A: = ^/kf + + A';| and ji denotes 
the spherical Bessel function and $t = $/$(0). From now on we assume the scale- 
invariant Harrison-Zel'dovich spectrum (V^{k)=k'^~^ — const.) as the initial power 
spectrum as predicted by the standard inflationary scenarios. Let us estimate the 
angular scale below which the power spectrum is suppressed owing to the mode- 
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cutoff kcut = 271/ L. First of all we consider the transfer function Ti{k) 

21 + 1^ f^2^,^^^.^dk 



Q = J Tf{k)r{k)-, (4.82) 

which describe how the spatial information is contained in the angular power. For 
a given k, the contribution to the angular power on smallest angular scales comes 
from the OSW effect where the transfer function can be written in terms of the 
spherical Bessel function as Tf^^ (k) = ji^krio). Because ji(x)'s have the first peaks 
at X ~ 1 + /, the angular cutoff Icut is determined by Icut = ^nrjo/L — 1. For instance, 
hut ~ 5 for L = r/o = 2Hq^. Thus the large-angle suppression scale is determined 
by the largest fluctuation scale at the last scattering. On smaller angular scales 
I > hut, the second peak in the power corresponds to the fluctuation scale of the 
second eigenmode at the last scattering. This behavior is analogous to the acoustic 
oscillation where the oscillation scale is determined by the sound horizon at the last 
scattering. On angular scales larger than hut one observes another oscillation feature 
in the power for models with a smaller cell where kcutVo is sufficiently large which is 
apparently determined by the behavior of the first eigenmode. From the asymptotic 
form 

~ -sin(x-/7r/2), x » ^ ^\ (4.83) 
X 2 

one notices that the angular scale of the oscillation for asymptotic values 

kcutVo + is A/ = 2. For intermediate values of kcutVo, A/ take much larger 

values (see figure |47^) . 

So far, we have studied the effect of the non-trivial topology on the OSW con- 
tribution only which is sufficient for constraining the topology of the "standard" 
toroidal model with Qtot = 1- However, for low matter density models, one cannot 
ignore the (late) ISW contribution which is generated by the decay of the grav- 
itational potential at the A-dominant epoch. The crucial point is that they are 
produced at later time 1 -|- 2; ~ (^a/^o)^''^ well after the last scattering. Although 
it is impossible to generate fluctuations beyond the size of the cell (in 3-dimensional 
sense), that does not necessarily mean that any fluctuations on large angular scales 
(in 2-dimensional sense) cannot be produced. Suppose that a fluctuation is pro- 
duced at a point nearer to us, then the corresponding angular scale becomes large 
if the background geometry is flat or hyperbolic. Therefore, one can expect that 
the suppression on large angular scales owing to the mode-cutoff is not stringent for 



low matter density models. As shown in figure the angular powers for a model 
with (l^A; ^m) = (0.9, 0.1) are almost flat. In contrast to the "standard" model, the 
transfer function for low matter density models distributes in a broad range of k that 
implies the additional late production of the fluctuations which contribute to the 
angular power. Surprisingly, in low matter density models with small volume, the 
slight excess power due to the ISW effect is cancelled out by the moderate suppres- 
sion owing to the mode-cutoff which leads to a flat spectrum. However, as observed 
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in the "standard" toroidal model, the power spectra have prominent oscillating fea- 
tures. The oscillation scale for / < Icut is determined by the first eigenmode. The 
peaks in the angular power correspond to the first SW ridge and the first and the 
second and other ISW ridges. Is such an oscillating feature already ruled out by the 
current observation? We will see the results of our Bayesian analyses for testing the 
goodness-of-fit to the COBE data in the next section. 



4.3.2 Closed hyperbolic models 

Assuming that the initial fluctuations obey the Gaussian statistic, and neglecting 
the tensor-type perturbations, the angular power spectrum Ci for closed hyperbolic 
models can be written in terms of the expansion coefficients C,uim as 



(4.84, 



where 



1 n° d^t 

3 ir,. drj 



(4.85) 



Here, z/ = \Jk? — 1, V^i^u) is the initial power spectrum, and i]^ and rjo are the 
last scattering and the present conformal time, respectively. X^^i denotes the ra- 
dial eigenfunctions of the universal covering space and ^^im denotes the expansion 
coefficients, $t = ^/^(O) and Vol(M) is the volume of the space. The (extended) 
Harrison-Zel'dovich spectrum corresponds to Vq>{v) = v^~^ = Const, which we shall 
use as the initial condition. Using Weyl's asymptotic form one can easily show that 



in the limit z/ ^ oo, ( [4.85|) coincides with that for the infinite counterpart provided 
that < l^jy^mP >oc Because we have assumed that the normalization factor 

in < l^jy^mP > does not depend on the volume of the space, the volume factor is 
necessary for the manifold-invariant form which cancels out the volume factor term 
in Weyl's formula. However, it might not be necessary if < l^^^mP > is proportional 
to l/Vol(M). As we have seen in chapter 3, the ratio of the normalization factor in 
< l^i/imP > is 1-03 for the Weeks and the Thurston manifolds which is close to the 
inverse ratio in the volume 1.04. 

In contrast to the previously studied closed fiat toroidal models closed hyperbolic 
spaces are globally inhomogeneous. Therefore, the angular power spectrum depends 
on the place in the space. It seems that evaluation of the full spatial dependence is 
intractable. Fortunately, we can understand the spatial dependence in terms of the 
pseudo-Gaussian random property of the eigenmodes. As we have seen in chapter 
II, each set of expansion coefficients ^^im at a particular place with particular ori- 
entation of the coordinates can be well approximated as a "realization" of random 
Gaussian numbers. Hence the uncertainty in the power is well described by the 
statistical property of the random Gaussian numbers. The uncertainty in the power 
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Figure 4.3: Geometric Variance 
The geometric variance in tlie angular power spectra AT; = 

^J l{l + 1)Ci/{2tt) (normalized by Ar;(^Q^), l{max) = 100,70,50 for 0^ = 
0.1,0.2,0.3 — 0.6, respectively) for the Thurston models (volume=0.98). 
The ensemble averaged values (boxes) are compared with those at a place 
where the injectivity radius is locally maximal (diamonds) and those for the 
infinite hyperbolic models (stars) for different matter density at present. 
Cosmological constant is not included. Two-sigma geometric variance is 
shown in vertical lines which has been obtained by 500 random realiza- 
tions for the expansion coefficients. The initial fluctuations are assumed 
to be adiabatic and obey Gaussian distribution specifled by the extended 
Harisson-Zel'dovich power spectrum. Since we neglected the 'early' ISW 
effect, on intermediate scales (/ > 20) the plotted powers are slightly sup- 
pressed than the expected values. 

owing to the globally inhomogeneous geometry can be called the geometric variance. 
Because we do not have any information about the position in the spatial hypersur- 
face, any particular places do not have priority over others. However, there are some 
places at which one observes symmetries of the background geometry. For instance, 
in the Thurston manifold having a D2 symmetry group there are three axes that 
are fixed by a tt rotation. If we were in a position near these axes by chance then we 
would observe Z2 symmetry in the sky where the pseudo-Gaussian behavior breaks 
down since there are non-trivial correlations in the expansion coefficients. In gen- 
eral, the chance that we actually observe any symmetries (elements of the isometry 
group of the manifold or orbifold or of their finite sheeted covers) is expected to be 
very low since a fixed point by an element of the isometric group is either a part 
of a line (for instance, an axis of a rotation) or an isolated point (for instance, a 
center of an antipodal map). On the other hand, the symmetric structures can be 
a crucial sign which implies the non-trivial topology. On large angular scales, the 
contribution to the power only comes from the eigenmodes with long wavelength. 
Therefore, the degree of freedom which is equivalent to the number of the terms in 
the summation that gives the power is small. In figure (|4.3|) , one can notice a large 
uncertainty in the power on large angular scales. However, the geometric variance 
is always smaller than the "initial" variance owing to the uncertainty in the initial 
conditions since the degree of freedom {k,m) is larger than 21 + 1. Therefore, the 
net cosmic variance (initial variance + geometric variance) cannot be significantly 
greater than the values for the infinite counterparts. On small angular scales, the 
contribution to the power comes from a number of eigenmodes which causes an in- 
crease in the degree of freedom and the geometric variance becomes small. In that 
case, one cannot distinguish the closed model with the infinite counterpart. When 
the power coincides with that for the infinite counterpart, the geometric variance 
becomes negligible. Thus we have two kinds of uncertainty: one is the initial variance 
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Figure 4.4: Comparison with COBE data 

Plots of the angular power AT; = y''^(^~+^O^V(^^ for the Thurston models 
with different matter density parameters in comparison with the COBE- 
DMR data analyzed by Gorski (diamonds) and Tegmark (boxes). The 
light-gray band denotes one-sigma "initial" variance around the power av- 
eraged over initial conditions for a particular choice of the observation 
point (where the injectivity radius is locally maximal) and the darkgray 
band represents the net variance of one-sigma "initial" variance and one- 
sigma "geometric" variance around the power averaged over both initial 
conditions and positions. 

Figure 4.5: Suppression in Large-angle Power for CH models 

Suppression in large-angle power AT; = ^J l{l + 1)Ci/{2tt) for 5 closed 
hyperbolic models with different volume (name,volume)=A:(m003(3,- 
1) (Weeks), 0.94), B:(m010(-l,3),1.9), C:(m082(-2,3),2.9), D:(m288(- 
5,1), 3. 9) E:(s873(-4,1), 4.9) in comparison with the infinite counterpart 
(denoted as O). All the plotted values are normalized by AT20 for the 
infinite hyperbolic model. First eigenvalues ki for the five models are 
plotted as horizontal lines (lower figures). The unit of k is equal to the 
inverse of the present curvature radius. 

arising from the uncertainty in the initial conditions; the other one is the geometric 
variance arising from the uncertainty in the position of the observer owing to the 
global inhomogeneity in the geometry. The latter effect is negligible on small scales 
compared with the size of the fundamental domain. 

Compared to flat A models, the ISW contribution is much significant for hyper- 
bolic (open) models since the curvature dominant epoch comes 1 -|- z ~ (1 — Qq)/Qq 
earlier than the A dominant epoch 1 + z (^a/^o)^^'^- Therefore, the suppression 
in the large-angle power owing to the mode-cutoff is less stringent compared with 
closed flat A models. Similar to A models, for a given angular scale the transfer 
function for hyperbolic low matter density models distributes in a broad range of k 
owing to the additional production of the fluctuations in the late epoch. As shown 
in ( |4.5|) for Qq = 0.2 the contribution from the ISW effect dominates the large-angle 
power on large angular scales / < 20. Therefore, the mode-cutoff owing to the finite 
size of the fundamental domain does not lead to a significant suppression even for 
small CH models like the Weeks or the Thurston models. The angular cutoff scale 
Icut is determined by the angular scale of the fluctuations of the flrst eigenmode at 
the last scattering. The cutoff scale is approximately given by the volume of the 
space V and the smallest non-zero wavenumber ki 

k 

lent = -J Vw-H(sinh(2(i?* + rave)) - sinh(2(i?, - rave)) - 4r„„e)) (4.86) 
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where v = 7r(smh(2rat,e) — 2rave) and i?* is the comoving radius of the last scattering 
surface. For the Weeks manifold [v = 0.94 and ki = 5.27), Icut = 26 for Qq = 0.2 
and Icut = 7 for Qq = 0.6. On angular scales I > 2lcut, the power for a closed 
model coincides with that for the infinite counterpart. Compared with closed flat 
models the "topological oscillation" on scales Icut < ^ < '^hcut is not prominent. The 
oscillation feature in the power is purely determined by the first eigenmode. The 
peaks in the power corresponds to the OSW ridge and the ISW ridges. 

As shown in figure ([4.4|) , the quadrapole in the COBE data is very low and there 
is a peak at / ~ 4. A Thurston model with = 0.6 has the first peak in this scale. 
Therefore, the fit to the observed power is much better than do any FRWL models. 
For the infinite hyperbolic models, the ISW contribution leads to the excess power 
on large angular scales which is not compatible with the behavior of the observed 
power on very large angular scales. Interestingly, the angular suppression owing to 
the mode-cutoff effectively reduces the excess power as observed in closed fiat A 
models. Thus we expect that the COBE constraint for CH models are less stringent 
than the "standard" flm = 1 closed flat toroidal models which are claimed to have 
been "ruled out". In the next section, we will statistically test the likelihoods of 
closed models with non-trivial topology using the COBE-DMR data. 
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4.4 Bayesian analysis 

In general, the covariance in the temperature at pixel i and pixel j in the sky map 
can be written as 

=< TiTj >=J2< (^irnai'm' > WiWi>Yim{n^)Yi'm'{nj)+ < N,N^ > (4.87) 
I 

where aim is an expansion coefficient with respect to a spherical harmonic Yim, <> 
denotes an ensemble average taken over all initial conditions, positions and orienta- 
tions of the observer, Tj represents the temperature at pixel i, Wf is the experimental 
window function that includes the effect of beam-smoothing and finite pixel size, rij 
denotes the unit vector towards the center of pixel i and < NiNj > represents the 
noise covariance between pixel i and pixel j. If the temperature fluctuations form 
a homogeneous and isotropic random Gaussian field then the covariance matrix can 
be written in terms of the power spectrum Ci as 

Mij = i- ^(2/ + l)W^CiPiini ■ nj)+ < NiNj > (4.88) 

where Pi is the Legendre function. Then the probability distribution function of the 
pixel temperature T for the Gaussian field is 

f(f\Ci) = ^—j. exp(-f^-M-\Ci)-f\ (4.89) 

^ ' ^ (27r)^/2deti/2^(G) V2 ^ v ; 

where N is the number of pixels. Bayes's theorem states that the probability distri- 
bution function of a set of parameters A given the data T is 

f{A\f) (X f{f\A)f{A). (4.90) 

If we assume a uniform prior distribution, i.e. taking f{A) to be constant, the 
probability distribution function of a power spectrum Ci is then 

A(Ci\f) oc expf -f^-M-^Q) -f). (4.91) 

^ ' ^ det^/'M(Q) V2 ^ ^) J ^ ^ 

In the following analysis, we use the inverse-noise-variance-weighted average map 
of the 53A,53B,90A and 90B COBE-DMR channels. To remove the emission from 
the galactic plane, we use the extended galactic cut (in galactic coordinates) [p^] . 
After the galactic cut, best-fit monopole and dipole are removed using the least- 
square method. To achieve efficient analysis in computation, we compress the data 
at "resolution 6" (2.6°)^ pixels into one at "resolution 5" (5.2°)^ pixels for which 
there are 1536 pixels in the celestial sphere and 924 pixels surviving the extended 
galactic cut. The window function is given by Wi = GiFi where Fi are the Legendre 
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coefficients for the DMR beam pattern |85| and Gi are the Legendre coefficients 
for a circular top-hat function with area equal to the pixel area which account for 
the pixel smoothing effect (which is necessary for "resolution 5" pixels since the 
COBE-DMR beam FWHM is comparable to the pixel size) To account for the 
fact that we do not have useful information about monopole and dipole anisotropy, 
the likelihood must be integrated over Co and Ci in principle. However, in practice 
we set Co = Ci = lOOmK^ which renders the likelihood insensitive to monopole 
and dipole moments of several mK. We also assume that the noise in the pixels is 
uncorrelated from pixel to pixel which is found to be a good approximation . 



4.4.1 Closed flat models 

Before applying the method to the flat toroidal models one must be aware that 
a fiat 3-torus we are considering is globally anisotropic although it is globally ho- 
mogeneous. In contrast to the standard infinite models, the fluctuations form an 
anisotropic Gaussian field for a fixed orientation if the initial fluctuation is Gaussian. 
In other words, for a given /, a set of ai^s {—l<m<l) are not 2/ + 1 independent 
random numbers. In order to see this, we write a plane wave in terms of eigenmodes 
in the spherical coordinates, 

e^^-^ = hum ji{kx)YiUn), him = ^7t{{}' Yi*^{k), (4.92) 

Im 

where k denotes the unit vector in the direction of k. Then we have 

aim = J2^(^) bklmFkl (4.93) 

k 

where F^i is given by ([4.81|) . Because aim is linear in $(fc), it is Gaussian. However, 
they are not independent since bkim^s are proportional to spherical harmonics. The 
statistical isotropy is recovered iff the size L of the cube becomes infinite where k 
takes the whole values for a given k. If one marginalizes the likelihood with respect 
to the SO (3) transformation, or equivalently, the orientation of the observer, the 
distribution function of a;^ becomes non-Gaussian since a;^ is written in terms of 
a sum of a product of a Gaussian variable times a variable that is determined by a 
spherical harmonic Yjm(k) where k is a random variable which obeys a uniform dis- 
tribution. Thus in principle, in order to compare a whole set of fluctuation patterns 
over the isotropic ensemble one should use the likelihood function which is different 
from the Gaussian one. 

Nevertheless, we ffist carry out a Bayesian analysis using the Gaussian likelihood 
function ( ^.91| ) which depends only on C;'s (which is SO (3) invariant) in order to 
estimate the effect of the mode-cutoff imposed by the periodic boundary conditions. 

From figure [4.6| one can see that the constraint on the size of the cell is less 
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Figure 4.6: Relative Likelihoods of Toroidal Models I 

The likelihoods of the toroidal models with (Qrm^A) = (1-0,0) (diamond) 
and (0.1,0.9) (box) relative to the infinite models with the same density pa- 
rameters are plotted. The likelihoods are marginalized over the quadrapole 
normalization Q = (5C2/(47r))^/^. An isotropic Gaussian approximation 
has been used for computing the likelihoods. The COBE-DMR data is 
compressed to 924 pixels at "resolution 5" . 



stringent for the low matter density model as expected from the shape of the 
power spectrum. The effect of the suppression of the large-angle power is not 
prominent unless the cell size is sufficiently smaller than the observable region 
(L < O.SHq^ = 0.08 X 2i?^, for = 0.1). However, the likelihood function varies 
rapidly as the size increases since the power is jagged in /. Unfortunately, the peaks 
in the power at / ~ 6 and / ~ 11 which correspond to the first and the second ISW 
ridge give a bad fit to the COBE data for L ~ 1.7Ho = 0.17 x 2R^,. However, the 
parameter range which gives a good fit to the data is wider for the low matter den- 
sity model. It should be emphasized that for both fiat models there is a parameter 
region in which the fit is much better than the infinite counterpart. For example, the 
peaks at / ~ 4 and / ~ 9 in the angular power of the low matter density model with 
L ~ 2.8Hq^ give a much better fit to the COBE data than the infinite counterpart. 
In fact the quadrapole component in the COBE data is very low and the angular 
power is peaked at / ~ 4. For infinite fiat A-models with a scale-invariant initial 
spectrum such features can be a problem since the ISW contribution gives an excess 
power on large angular scales. 

Next, we carry out a full Bayesian analysis in which all the elements of < 
aimdi'm' > are included. Because of the limit in the CPU power we further com- 
press the data at "resolution 5" pixels to "resolution 3" (20.4°)^ pixels in galactic 
coordinates for which there are 60 pixels surviving the extended galactic cut. Al- 
though the information on smaller angular scales / > 10 is lost, they still provide 
us a sufficient information for discriminating the effect of the non-trivial topology 
since it is manifest only on large angular scales. The computation has been done 
in a similar manner as previous analysis except for the covariance matrix for which 
we use (|4.87| ) instead of ( [4.88|) . The likelihoods A are computed for a total of 2000 
random orientations for each model. The approximated likelihoods which depend 
on only the power spectra are also computed for comparison. As shown in figure 
[4.7| , the discrepancy between the likelihood marginalized over the orientation of the 
observer and the approximated likelihood is prominent for models with small volume 
in which the effect of the non-trivial topology is significant. Let us estimate the size 
of the cell for which the effect becomes insignificant. For the "standard" model with 
Qm = 1-0, the volume of the observable region is 47ri?^/3~34 which gives the critical 
scale Lc^3.2Hq^ for which the volume of the cell is comparable to the volume of the 
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Figure 4.7: Relative Likelihoods of Toroidal Models II 

The hkelihoods of the toroidal models relative to the infinite models with 
the same density parameters marginalized over 2000 orientations (star) and 
that using only the power spectrum (diamond). The inverse- noise- variance- 
weighted average map of the COBE-DMR data is compressed to 60 pixels 
(resolution 3). All the likelihoods are also marginalized over the quadrapole 
normalization Q = (5C2/(47r))-'^/^. 

observable region at present. For the low matter density model a significant amount 
of large-angle fluctuations are produced at the A-dominant epoch z^l. Therefore 
one should compare the volume of the sphere with radius r]{z = 0)—r]{z = l) ~ 0.9iJ(7^ 
to the volume of the cell instead of the observable region which gives Lc^1.5Hq^ . 
We can see from flgure that these estimates well agree with the numerical result. 

We can give two explanations for the discrepancy although they are related each 
other. One is the non-Gaussianity in the fluctuations and another one is the corre- 
lations between aim's owing to the global anisotropy in the background geometry. 

Suppose a random variable Z = XY in which X obeys a distribution function 
E{X) which is even and Y obeys a distribution function F{Y). Then the distribution 
function G of Z is given by 

GiZ) = I ^i^FiY)dY, (4.94) 

which is apparently even. Because the fluctuations are written in terms of a sum 
of products of Gaussian variable $(fc) (with zero average) times a non-Gaussian 
variable Yjm(k), the skewness in the distribution function of aim marginalized over 
the orientation is zero although the kurtosis is non-zero. 

The correlations in aim's are the consequence of the gap between the degree of 
freedom of (k) and (/, m). The degeneracy number in a fc-mode (=the number of the 
direction k) is much less than the number of relevant "quantum numbers" {I, m)| 
if the scale (=27r//c) is comparable to the length of the side L. Taking an ensemble 
average over the initial condition we have 

< aimai'm' >= jjj^V^{k)Yimik)Yi,m'ii^)FkiFki', (4.95) 

where Fki is given by (|4.81| ). The sum does not vanish in general when (/,m) ^ 
(/', m') which is the consequence of the global anisotropy in the background geome- 
try. However, if one takes an average of ( [4.95|) over k, one flnds that all off-diagonal 
elements vanish because of the orthogonality of the spherical harmonics. Similarly 



^Here we assumed that for each degenerated mode specified by k, $(fc) are independent Gaussian 
variables. 
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one can consider 4-point correlations < ai^rn-i_o,i^m20'i3m30,i^m4 >■ In this case, all the 
off-diagonal elements {li,mi) ^ {lj,mj) do not necessarily vanish even if one takes 
an average over k. Thus the non-Gaussianity for the flat toroidal models contrasts 
sharply with the one for the compact hyperbolic models in which the pseudo-random 
Gaussian property of the expansion coefficients bkim (which are obtained by expand- 
ing the eigenmode in terms of eigenmodes in the universal covering space) EO, ET 



renders off-diagonal elements always vanish if an average is taken over the position 
of the observer (although the kurtosis is non-zero). The difference can be attributed 
to the property of eigenmodes. In a less rigorous manner the property of the eigen- 
modes which are projected onto a sphere with large radius can be stated as follows: 
Projected eigenmodes of compact hyperbolic spaces are "chaotic" whereas those of 
compact flat spaces are "regular". 

As for the constraint on the size of the cell, we set a slightly severe condition 
for the low matter density model since the likelihood of the infinite counterpart is 
~ 10^^ of that of the infinite "standard" model due to a slight boost on the large 
angular scales caused by the ISW effect. Together with the previous analysis us- 
ing the data on the "resolution 5" pixels, the conditions on the relative likelihood 
logio(A/Aoo) > -2 and logio(A/Aoo) > -1 yield L > O.AOR^ and L > 0.22i?, for the 
"standard" toroidal model {flm, ^a) = (1-0, 0) and the low matter density toroidal 
model (Qmy^A) = (0.1,0.9), respectively. Here Aoo denotes the likelihood of the 
infinite counterpart with the same density parameters. The maximum number N of 
images of the cell within the observable region at present is 8 and 49 for the former 
and the latter, respectively. Note that the constraint on the "standard" toroidal 



model is consistent with the previous result|10, 11, 12 



4.4.2 Closed hyperbolic models 

For CH models, the temperature fluctuations for a homogeneous and isotropic en- 
semble can be written as a sum of modes which consist of a product of :a Gaussian 
random number $^(0) originated from quantum fluctuations; a Gaussian pseudo- 
random number ^,^1^ originated from geometry. Thus the amplitude of the temper- 
ature fluctuations in (/,m) space is written as 

aim = J2M0)UmF,i (4.96) 

which are non-Gaussian. In contrast to flat models, correlations between 2/ + 1 vari- 
ables for a given I are very small since the case for which one observes any symmetric 
structure in the temperature fluctuations is rare. 

First of all, we consider the effect of the suppression on large angular power. 
To do so, we fix the initial conditions as $^(0) = Const. jiyiy^ + 1)) so that they 
are not random number but satisfying the standard assumption that the power 
spectrum is described by the extended Harrison-Zel'dovich spectrum. Then the 
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Figure 4.8: Relative Likelihoods of Smallest CH Models 

Plots of the ratio of approximated hkelihoods for two smaUest CH mod- 
els (Weeks and Thurston) to a hkelihood for the Einstein-de Sitter model 
(r^TO = 1.0) with scale invariant initial spectrum (n = 1). All the likeli- 
hoods are marginalized over the normalization of the power. Here we have 
assumed "hxed initial fluctuations" <1>^(0) = Const. + 1)) for CH 
models. 



Figure 4.9: Relative Likelihoods of Open and Thurston Model 
Plots of the ratio of approximated likelihoods for the inhnite hyperbolic 
"open" models (n = 1) to one for the Einstein-de Sitter model (n = 1) 
(left) and likelihoods for the Thurston models relative to one for infinite 
hyperbolic models (right). All the likelihoods are marginalized over the 
normalization of the power. We have assumed "fixed initial fiuctuations" 
$^(0) = Const. /{v{v'^ + 1)) for the Thurston models. The slight improve- 
ment in the likelihood of infinite hyperbolic models with Vtrn < 0.1 is caused 
by the absence of the supercurvature modes. For the Thurston models with 

= 0.1 0.3 the likelihoods are significantly improved. 



fluctuations approximately form a Gaussian random field owing to the pseudo- 
Gaussianity of the eigenmodes and the property of the fiuctuations is fully specified 
by the corresponding angular power Ci. As in the previous sections, we use the 
inverse-noise-variance-weighted average map of the 53A,53B,90A and 90B COBE- 
DMR channels at 924 pixels ((5.2°)^) surviving the extended galactic cut. As is 
expected from the shape of the angular powers, for a wide range of parameters 
{Qm + > 0.1) the likelihoods for the smallest (for manifolds) CH models are 
better than one for the Einstein-de Sitter models with the scale-invariant spectrum 
(n = 1) where AT/ = (/(/ + l)C//2'7r)"^/^ is almost constant in /. One can see that the 
fit to the COBE data is much better than the Einstein-de Sitter models for three 
parameter regions:l. Qm = 0.5 ~ 0.7 with small for which the angular power is 
peaked at Z ~ 4 which corresponds to the first ISW ridge of the transfer function 
;2. Qrn = 0.85 ~ 0.9 where the angular scale I ~ 4 corresponds to the OSW ridge;3. 
flm ~ 0.2 where the slope of the power on large angular scales / < 10 fits well with 
the data. 

For low matter density models, the ISW effect leads to an excess large-angle 
power. Therefore the likelihoods of low matter density models Q^n = 0.1 ~ 0.3 with 
infinite volume are relatively small because of the low quadrapole moment in the 
COBE data. However, for small CH models, the excess power owing to the ISW 
effect is mitigated by suppression owing to the mode-cutoff. Therefore, likelihoods 
for small CH models with low matter density are significantly improved compared 
with the infinite counterparts (figure [4.9|) . 

As we have seen in chapter 3, the property of eigenmodes is associated with the 
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Figure 4.10: Relative Likelihoods of 27 CH Models 

Plots of the ratio of approximated likeUhoods for 27 CH models with volume 
(0.94—6.15) to a likelihood for the Einstein-de Sitter model (f^m = 1-0) with 
scale invariant initial spectrum (n = 1). All the likelihoods are marginalized 
over the normalization of the power. Here we have assumed "fixed initial 
fluctuations" ^1(0) = Const./ {u{u'^ + 1)) for CH models. 



volume and the "anisotropy" which can be characterized by the length of the shortest 
periodic geodesic Imin or diameter d. For instance, as can be seen in figure ( [4.8| ) the 
likelihoods of the Weeks models (volume=0.94, /mm = 0.585 and (i = 0.82) are similar 
to those of the Thurston models with similar geometrical property (volume= 0.98, 
^mm = 0.578 and (i = 0.81). If Imin is comparable to the "average" radius dave defined 
as the radius of a sphere with the same volume of the CH space then the statistical 
property of CMB anisotropy is described by Weyl's formula. 

For CH models with large volume, the likelihoods will converge to those of 
the infinite counterparts although the convergence rate depends on the cosmolog- 
ical parameters. One can see in figure (|4lO|) that the conspicuous difference for 



Vim = 0.1 — 0.3 still persists for volume~ 6 whereas such difference is not observed 
for nearly fiat cases (f2t = f2m + fiA = 0.9). The difference depends on the number N 
of the copies of the fundamental domain in the observable region at present. Sup- 
pose that = then the comoving radius of the last scattering surface in unit of 
the present curvature radius Rcurv is approximately given by 

Rlss = Rcurv cosh-i(2/fi^ - 1). (4.97) 

The comoving volume of the ball inside the last scattering surface in hyperbolic 
space is 

V = nRl^,,{smh{2RLSs/ Rcurv) - 2RlSs/ Rcurv) (4.98) 

For example, v ~ 490/2^^^^ for (fi^, ^a) = (0.2, 0) whereas v ~ Rl,urv (^m, ^a) = 
(0.9,0). Apparently, in nearly fiat cases, the volume of the spatial hypersurface 
should be smaller than R^^urv detecting the effect of the non-trivial topology. 
However, inclusion of the cosmo logical constant leads to an increase in Rlss because 
of a slow increase in the cosmic expansion rate in the past. Thus it allows a large 
N , for instance, is 8.7 for a Weeks model with f^A = 0.7 and Vim = 0.2 whereas 
A^ = 1.2 for one with r2A = and fim= 0.9. 

Next, we consider the effect of the non-diagonal elements which we have ne- 
glected so far. The likelihood for homogeneous and isotropic ensemble is obtained 
by marginalizing the likelihoods all over the positions Xo^^ and the orientations a of 
the observer, 

A = j ^ dy.obsdaK{y.obs,Ci). (4.99) 
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Plots of the ratio of "rigorous" likelihoods incorporating the effect of in- 
homogeneity and anisotropy of the background geometry for two closed 
multiply connected models relative to one for the infinite counterparts in 
ascending order(full curve) for 2000 random realizations of orientation: 
closed flat toroidal model with = 0.1, = 0.9, L = 2.0Hq^ (left) 
and 1500 random realizations of position and 10 realizations of orientation 
;a Thurston model with = 0.4 (right). The dashed lines denote the en- 
semble averaged values. The number of the copies of fundamental domain 
inside the last scattering surface is 64.7 for the closed flat toroidal model 
and 72.3 for the Thurston model. 



We assume a constant distribution for the volume element ^(hi.ohs and da which is 
the volume element of a Lie group 5*0(3) with Haar measure. Assuming that the 
initial fluctuations are Gaussian, the likelihood A(xof,s,a)) is given by (|4.87| ) and 
( [4.91|) where < aimdi'm' > is written in terms of the expansion coefficients ^^im and 



the initial fluctuation as 

7^ (^) 

< O'lmO'l'm' >0C . ^. C,ulmC,ul'm'FulFiyi', (4.100) 

//(z/"' + 1) 

where <> denotes an ensemble average taken over the initial conditions, Vqt^u) is 
the initial power spectrum and F^i describes contribution from the OSW effect and 
the ISW effect, respectively. iuVm'^^ are functions of Xq^^ and a. As a typical ex- 
ample we choose the Thurston manifold m004(— 5, 1) which is the second smallest 
known CH manifold. In order to compute likelihoods, we use a compressed data 
at "resolution 3" (20.4°)^ pixels in galactic coordinates for which there are 60 pixels 
surviving the extended galactic cut as in closed flat cases. 

In comparison with the flat toroidal model with low matter density (^a = 
0.9, VLrn = 0.1, = 65) with similar value in N the likelihoods of a Thurston model 
(Qa = 0,f2m = 0.4, A^ = 72) distribute in a wide range of values(flgure [4.4.2|) . 



For the toroidal model, the likelihoods are larger than the average for 11 percent 
out of the total of 2000 orientations which are randomly chosen. On the other 
hand, for the Thurston model, only 0.09 percent of the total of 1500 positions and 
10 orientations are larger than the average. We have also computed likelihoods 
of the Thurston model for 5100 positions and 40 orientations. Then the percent- 
age has reduced to 0.02 from 0.09. The log of ratio of the likelihood marginalized 
over 5100 positions and 40 orientations to the one of the inflnite hyperbolic model 
(Vim = 0.4) is Logio(A/A(oo)) = —1.4). The maximum value for the likelihood is 
Logio(^('^ax)/A(oo)) = 3.6). Thus the flt to the COBE data at a particular place 
with a particular orientation is usually not good for the Thurston model but it does 
not constrain the model. For a particular set of position and orientation likelihoods 
are much better than the standard inflnite counterpart. 



The likelihood analyses in [BO, 81] are based on correlations for a particular 



choice of position Q where the injectivity radius is locally maximal with 24 ori- 
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entations. It should be emphasized that the choice of Q as an observing point is 
very special one. For instance it is the center of tetrahedral and Z2 symmetry of 
the Dirichlet domain of the Thurston manifold although they are not belonging to 



the symmetry of the manifold In mathematical literature it is standard to 

choose Q as the base point which belongs to a "thick" part of the manifold since 
one can expect many symmetries. However, considering such a special point as the 
place of the observer cannot be verified since it is inconsistent with the Copernican 
principle. Because any CH models are globally inhomogeneous, one should compare 
fluctuation patterns expected at every place on the space. On the other hand, the 
previously studied closed flat models in which the spatial geometry is obtained by 
gluing the opposite faces by three translations are globally homogeneous^ for which 
one does not need to compute fluctuations at different positions. Any CH models 
cannot be constrained until all the possible fluctuations are compared to the data. 

The result is not surprising if one knows the pseudo-random behavior of eigen- 
modes on CH spaces. For each choice of the position and orientation of the observer, 
a set of expansion coefficients ^^im of eigenmodes is uniquely determined (except 
for the phase factor), which corresponds to a "realization" of independent random 
Gaussian numbers. By taking an average over the position and the orientation, the 
non-diagonal terms proportional to < iimivm' ^ I' ,m ^ m' vanish. In other 
words, a set of anisotropic patterns all over the place in the CH space comprises an 
almost isotropic random field. Consider two realizations A and B of such an isotropic 
random field. The chance you would get an almost similar fluctuation pattern for A 
and B would be very low but we do have such an occasion. Similarly, the likelihood 
at a particular position with a certain orientation is usually very low but there are 
cases where the likelihoods are considerably high. Thus we conclude that the COBE 
constraints on small CH models are less stringent compared with the "standard" flat 
toroidal model with = 1.0) as long as the Gaussian pseudo-randomness of the 
eigenmodes holds. 



*As described in chapter 1, other closed flat spaces are globally inhomogeneous. 
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Chapter 5 

Search for Topology 



Clever persons cannot fall in love. Love is blind. To become a scientist, one 
must love Mother Nature. She would open her heart to only those who truly 
love her. 

(Torahiko Terada, 1878-1935) 

Suppose that the spatial hypersurface of the universe is muhiply connected on scales 
smaller than the comoving distance to the last scattering surface, then, in principle 
we would be able to detect the signature of the non-trivial topology of the spatial 
geometry which would appear as the periodicity or non-Gaussianity in the CMB sky 
map. 

The periodic structure in the CMB sky map would certainly be the direct evi- 
dence of the non-trivial topology. Assuming that the temperature fluctuations are 
produced entirely at the last scattering surface(determined by the OSW effect only), 
the periodical structure would appear as a pair of circles in the sky on which the 
fluctuation patterns are identified [|]. In practice the effect of the periodical struc- 
ture in the CMB would appear as an anomalously large correlations on a particular 
set of directions in the sky [Q. Suppose a sphere with thickness comparable or less 



than the typical correlation length of the relevant fluctuations in the universal cov- 
ering space. Then one can expect a number of copied "thick" spheres that intersect 
the original one at a particular place. In other words, one would expect a number 
of copies of a particular realization of fluctuations in different directions in the sky. 
Note that the set of the points with large correlation does not necessarily belong to 
the "matched circles" because of the thickness of the spheres. 

Another feature is the Non-Gaussianity in the CMB temperature fluctuations 
which is the consequence of the break of the global homogeneity or isotropy. For 
a particular choice of the position and the orientation of the observer, the temper- 
ature fluctuations form an anisotropic Gaussian random field. If one marginalizes 
the likelihood over the positions and the orientations, then the fluctuations are re- 



garded as an isotropic non- Gaussian random field] 89]. Because we do not have any 



information of the position and the orientation of the observer we should give the 
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same chance for each set of position and orientation. In contrast to cosmic string 
models, the non-Gaussianity only appears on large-angular scales assuming that the 
spatial hypersurface is multiply connected on scale not significantly smaller than the 
present horizon scale. Furthermore it does not depend on the type of the promor- 
dial perturbation(adiabatic or isocurvature). In some non-standard (with multiple 
scalar fields) inflationary models, non-Gaussianity in the initial fluctuations are also 



expected but with isocurvature modes in most cases|9y]. It seems that one needs 
unnatural conditions for the potential(not smooth) and the infiaton for generating 
adiabatic non-Gaussian fiuctuations|^. The detection of adiabatic non-Gaussianity 
will certainly be the positive sign of the non-trivial topology of the universe. 

As we shall see, in order to search for the periodic patterns, one needs a sky 
map with high signal-to-noise ratio. Therefore, absence of such periodical structure 
in the COBE-DMR data in which the noise is comparable to the signal (on angu- 
lar scales ~ 10°) does not imply the trivial topology (=simply connectivity) of the 
spatial geometry. In near future, satellite missions MAP and Planck will provide us 
CMB sky maps with high resolution and high signal-to-noise ratio compared with 
COBE. Using such maps we may detect the imprint of the non-trivial topology if 
we live in a "small universe" . 



5.1 Circles in the sky 

Suppose a "bubble" B in the form of a 2-sphere with small radius in a multiply 
connected constantly curved space M. As one inflates B in M, B will eventually 
meet itself when the radius is equal to a half length of a closed piece-wise geodesic 
curve (which cannot be contracted to a point by continuous deformation). Let us 
further assume that B is able to cross itself without any deformation. Then the 
self-intersection becomes a circle C when inflated further. For an observer sitting 
at the center of B, the shortest geodesic distance from the center O to a point on 
C is equivalent to the radius. Therefore, the length of a geodesic segment si that 
connects O and a point p on C is equal to that of another geodesic segment s2 that 
also connects O and a point p on C. In other words, we can see a pair of circles 
in two different directions which are actually two images of a circle C at the same 
time. Now let us consider the last scattering surface as B. If the injectivity radius 
is less than the radius of the last scattering surface, then one can find a pair (or 
pairs) of circles around which the temperature fluctuation patterns are identical. 
In the universal covering space, one can consider the copied last scattering surface 
g{B) by the element g of the discrete isometry group. Then the matched circles 
correspond to the intersection of B and g{B) and that of B and g~^{B){see figure 
p.l| ). It should be noted that the distance between the observer and the copied 
observer should be less than the diameter of the last scattering surface for detecting 
the matched circles. 
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Figure 5.1: Matched Circles 



For multiply connected models identical fluctuation patterns around a pair 
of circles might be observed. The circles correspond to the intersections of 
the last scattering surface and its clones. 

In order to statistically test the significance of the periodical structure, we need to 
develop a statistical tool. Let us choose the directions ni and n2 which correspond to 
the centers of circles Ci and C2 with angular radius a in the sky. We now introduce 
a statistic for the correlation of the temperature fluctuations AT along the circles 
which is known as a Pearson's (product-moment) correlation constant r 

,2... ^ (E,Ar(»,,)A7-(»..+,))^ 

where riji, nj2, ■ ■ ■ , rij^, njN+i — njijj — l, 2 denotes the directions pointing to pixels 
around Ci and C2 (numbering the pixels clockwise for Ci but counterclockwise for 
C2). From the Cauchy- Schwartz inequality, it is apparent that r ranges in [—1,1]. 
r — 1 corresponds to the perfect correlation while r = — 1 corresponds to the perfect 
anti-correlation. In spatially orientable models, for a pair of matched circles we 
expect a nearly perfect correlation r ~ 1 ii the noise is negligible. If we include 
spatially unorientable spaces, we should also test for correlation in which pixels are 
numbered in the same orientation (clockwise-clockwise) for both circles. 

In actual experiments, we should treat the noise that might smear the clear 
detection of the periodicity. Let us assume that the temperature fluctuations at each 
pixel i for a circle is written as a sum of independent Gaussian random numbers for 
signal s and noise 77 as 

ATim) ^Si + rii (5.2) 

which obey the distribution 



p{s,r]i,r]2) = , ,3)^ ^ exp 



2V 



(5.3) 



where a1 and denote the variance of the signal and the noise, respectively. 
Marginalizing over the signal, the distribution of the temperature fluctuations for a 
pair of matched circles C\ and C2 is given by the Gaussian distribution 



p(ATi, AT2) = ^ exp 



1 / ATI + - 2pAri ATa 



2 1 - p2 V a 



2 



(5.4) 
where 

_ cov(A7\.Ar2) ol e 
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and (7^ = 0"^ + 0".^, ^ = as/drj- Let Nc be the total number of pixels in each circle. 
Then the distribution function of the correlation coefficient r for a pair of matched 



circle is written as [92 



p(r) = n-\N^ - 1)(1 - r^)(^-3)/2(i _ ^2^nj. I d(3 







(cosh P — pr)^'=~^ 



For sufficiently large A''^ 



1 1 + r 

z=-ln^— (5.7) 

2 1 — r 



is known to obey the Gaussian distribution with average 



<.>=il„i±^ (5.8) 



and variance = l/{Nc — 2). Therefore, the distribution (|5.6|) is peaked at 



^2 

'"max = P = -j^ _|_ ^2 ' (^-9) 

where ^ denotes the signal-to-noise ratio. Thus the large signal-to-noise ratio implies 
that p{r) is peaked at near r = 1 which makes it easy to detect the signature of the 
non-trivial topology. As for the distribution of r for unmatched circles we have to 
consider the spatial correlations of the temperature fluctuations. However, for pairs 
of circles which are sufficiently separated each other, the spatial correlations can be 
negligible. Approximating the temperature fluctuation at each pixel as a random 
Gaussian number (p = 0), the distribution for unmatched circles can be obtained 
from ( ^TB| ) 

(l-r2)(^-3)/2r(Ar,/2) , , 

In this case, the variance of r is approximately given by 

^^'^"^ - Ar.r(Ar./2) ^ AT, " ^^'^^^ 

The distribution of r for unmatched circles is peaked at zero and the variance is 
inversely proportional to the number of pixels Nc on the circle. Thus an increase in 
Nc will also make it easier to detect the imprint of the non-trivial topology. 

Let us now estimate the effect of the noise in actual experiments. First we con- 
sider the detectability by using the COBE-DMR data for which the signal-to-ratio 
is approximately ^ ~ 2 on angular scales ~ 10°. The number of pixels on a 
circle with angular radius a is given by Nc{a) = 2-7rsinQ;/A^. Ignoring the galactic 
cut, for the COBE-DMR data with resolution = 10°, N^a) ~ 36 sin a and the 
total number of pixels in the sky is Nt ~ 4.1 x 10^. If one considers the galactic 
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Figure 5.2: Detectability in COBE Data 

Distribution n for unmatched circle pairs (dashed curve) and a matched 
circle pair (full curve) with angular radius a = 70° as measured by the 
COBE-DMR (usinff pixels (10°)2). n{ro)dr gives the number of pairs of 
circles with r = tq 



cut, it reduces to Nt ~ 3.0 x 10^. The total number of the circle with angular 
radius a is Nt and the total number of pairs of circles with angular radius a is 
Npt{a) ~ Nc{a)Nt{Nt — l)/2. As shown in figure 5^, one can see that the detection 
at most would be marginal for the COBE-DMR data. It should be noted that the 
above estimate is done under idealistic circumstances. In reality, we should also 
consider physical effects such as the ISW effect and the Doppler effects which might 
smear the signal of the periodicity. For low matter density models, the ISW effect 
significantly reduces the signal-to-noise ratio (if we regard the ISW contribution as 
the "noise") on large angular scale ~ 10°. 

Next, we estimate the detectability by the future satellite missions MAP (launched 
in 2001) and P/ancA; (launched in 2007). According to the technical information at 
||93|| , the MAP satellite will measure the temperature fluctuations covering the full 
sky at five frequencies (22GHz, 30GHz, 40GHz, 60GHz and 90GHz) with angu- 
lar resolution (in FWHM of the Gaussian central beam) (0.93°, 0.68°, 0.53°, 0.35°, < 
0.23°), respectively. The sensitivity at a pixel with angular scale (0.3°)^ is ~ 35fiK 
for all the frequencies. If the Galactic emission is negligible at high Galactic latitudes 
and frequency above 40GHz then the three highest channels will give the sensitiv- 
ity ~ 20 fiK since the noise of the combined maps can be approximately given by 
r] = (Yl Vr"^)'^^"^ where r]i denotes the noise for map i. On the angular scale 1°, cor- 
responding to the acoustic peak, the expected sensitivity will be ~ 6fiK. Because 
the observed amplitude of the acoustic peak at / ~ 200 is 60fiK ~ SOfiK |95 



we have ^ = 10 ~ 13. Planck will provide us sky maps with angular resolution 
5 to 10 arcminutes[|96|. Certainly the signal-to-noise ratio will be much improved 



compared with the data supplied by MAP. Thus the prospect for the detectability 
is far brighter for the future satellite missions. 

However, one should be cautious about the physical effects that will smear the 
clear periodical structure. If one uses the filtered sky map on 1° scale, the 'early' 
ISW effect owing to the decay of the gravitational potential at the radiation dom- 
inant epoch equality epoch and the Doppler effect owing to the dipole component 
of the temperature fluctuations at the last scattering cannot be completely negli- 
gible although they are sub dominant. Below the Jeans scale the photon pressure 
resists the gravitational compression of the photon-baryon fluid leading to driven 
acoustic oscillations. When the scale of the fluctuations enters the sound horizon 
the gravitational potential starts to decay because of the resistance of the photon 
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Figure 5.3: Detectability in MAP Data 

Distribution n as measured by the MAP (using pixels (1°)^) for unmatched 
circle pairs with the total number of pixels on the circle Np = 30, 60, 90 
(dashed curve) and a matched circle pair with Np = 90 (full curve) and 
Np = 60(dashed-dotted curve). ^ denotes the signal-to-noise ratio. The 
pair of numbers in the figure in the right represent Np and ^. n{rQ)dr gives 
the number of pairs of circles with r = rg 



pressure to the gravitational force leading to an enhancement in the amplitude for 
the adiabatic modes. At the start of the oscillation, the amplitude of the monopole 
and dipole increase with R = 3pb/4:p^ (where pb and p^ denote the baryon density 
and the photon density, respectively) owing to the reduction in pressure |73[|. As the 
gravitational potential starts to decay, the amplitude of the monopole or the dipole 
decreases. However, the amplitude of the dipole decreases rapidly compared to the 
monopole by a factor of ^/Sdrg/drj = (l+i?)~^/^ where rj denotes the conformal time. 
Let us estimate the difference in the amplitude. The radiation density is given by 
the present temperature T normalized by 2.77^ and the Hubble parameter h as 
2.38 X IQ-^h-^TlM. Hence for flh = 0.06, h = 0.75 and 2 (recombination) = 1300, 
the factor (1 -|- R)^"^ is ~ 0.8. However, one should also consider the effect of the 
zero point shift of the oscillation for the monopole owing to the gravity since the zero 
point shift increases the amplitude of the observed temperature fluctuations. Note 
that the shift is absent for the dipole since it is approximately given by the time 
derivative of the monopole. Assuming that R is constant in time, the zero point 



shift for the temperature fluctuation 60 + \1/ is — i?\l/[98|. In reality, one should 



also consider the effect of the evolution of R. If one assumes that the zero point 
shift is equal to the amplitude of the temperature fluctuation plus the gravitational 
potential at the last scattering(6o + ^)('7*) ~ ^('7*)/3, then the relative amplitude 
of the monopole plus the gravitational potential to the dipole (for idiabatic case) is 
approximately equal to ~ 2. Thus the Doppler effect is not completely negligible. 

As shown in figure |5.3| , we need approximately 60 pixels on the circle for ^ = 4 
whereas 90 pixels for ^ = 2 in order to detect one pair of matched circles by using 
the expected sky map with resolution (1°)^ supplied by MAP. The decrease in the 
signal-to-noise ratio implies that it is necessary to have a pair of matched circles 
with large angular radius. Thus if the observed peak in the angular power is actu- 
ally relevant to the acoustic oscillations, then it is possible to detect the non-trivial 
topology provided that the injectivity radius at the observation point is smaller than 
the radius of the last scattering and the angular radius of the circles are sufficiently 
large. 
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5.2 Non-Gaussian signature 



Let us write the temperature fluctuations in the sky in terms of real spherical har- 
monics Qim and real expansion coefficients bim as 

AT 

^ = ^W;m. (5.12) 

Im 

If one assumes the spatial homogeneity and isotropy on the background geometry 
then the temperature fluctuations in the sky form an isotropic Gaussian random 
field, in other words, the distribution of the expansion coefficients bi^ is given by 

fikm) = exp(-6L/2Q). (5.13) 

However, in locally isotropic and homogeneous models which are spatially multi- 
ply connected, Gaussian temperature fiuctuations for a particular choice of position 
and orientation of teh observer are no longer isotropic since the spatial hypersurface 
is not globally isotropic. In other words, the variance of fe^mS may depend on m 
or him^ may not be independent each other for a given /. Marginalizing the distri- 
bution over the position and the orientation the temperature fiuctuations becomes 
non-Gaussian. Assuming that the initial fiuctuations are Gaussian, then the distri- 
bution has vanishing skcwncss but non-vanishing kurtosis since aim is written as a 
sum of a product of a Gaussian variable times an expansion coefficient of an eigen- 
mode which can be considered as a "random" variable. Non-Gaussianity owing to 
the break of the global isotropy and homogeneity in the background geometry can 
be a prominent signature of the spatial non-trivial topology. 

In order to test for Gaussianity, it is necessary to study the ensemble average 
of higher-order correlations (A^-point correlations N>2). Assuming vanishing aver- 
age, a Gaussian distribution can be defined as one for which the odd correlations all 
vanish but the even correlations all factorize into products of 2-point correlations. 

As we have seen, a temperature fiuctuation in locally homogeneous and isotropic 
background geometry is written as a sum of a product of two factors each one of 
which corresponds to the initial condition and the geometric property, respectively. 
For CH models, the both factors can be described by random Gaussian numbers. Let 
us calculate the distribution function F{Z, az) of a product of two independent ran- 
dom numbers X and Y that obey the Gaussian (normal) distributions A"(X; 0, ux) 
and A^(y;0,(Ty), respectively, 

N{X-ti,a) = —^e-^^-'^^'/'^' . (5. 14) 

v27r(T 



Then F{Z — XY, az) is readily given by 

z/Y. n rr^Wfy. n. ^-v^ 



~ dY 
F{Z,az) = 2 / N{Z/Y,0,ax)N{Y,0,aY)^ 
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Figure 5.4: Non-Gaussian Distribution 

On the left, the distribution function F{Z, 1) for a product of two ran- 
dom Gaussian numbers is plotted in solid curves. On the right, the dis- 
tribution function G{Z, 1) (la = 1) of a sum of two random variables that 
obey l/\/2). The dashed curves represent the Gaussian distribution 
N{Z- 0,1). 



(5.15) 



where Kf){z) is the modified Bessel function. The average of Z is zero and the 
standard deviation satisfies az = As is well known, Kq{z) is the Green 

function of the diffusion equation with sources distributed along an infinite line. 
Although Kq{z) is diverged at 2; = its integration over (—00,00) is convergent. 
From the asymptotic expansion of the modified Bessel function 
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one obtains in the lowest order approximation 

1 



FiZ,a) 



-\Z\/a 



27Ta Z 



Z » 1. 



(5.16) 



(5.17) 



Thus F{Z, a) is slowly decreased in Z than the Gaussian distribution function with 
the same variance in the large limit while it has a prominent peak around zero. The 
skewness is zero but the kurtosis is positive 



/oo 
Z^F{Z,l)dZ -3 = 6. 
-00 



(5.18) 



For fiuctuations on large angular scales only the eigenmodes with large wavelength 
(= 27r//c) can contribute to the sum. Owing to the finiteness of the space, the 
number of eigenmodes which dominantly contribute to the sum is finite. Therefore, 
the fiuctuations become distinctively non-Gaussian. On small angular scales the 
number of eigenmodes that contribute to the sum becomes large and the distribution 
function converges to the Gaussian one as the central limit theorem implies. One 
can see from figure ^ (right) that the distribution function G{W, 1) of W = Zl-\- Z2 
where both Zl and Z2 obey F{Z, \/2) is much similar to the Gaussian distribution 
N{Z, 0, 1) than F{W, 1) although divergence at Z = still persists. 



5.2.1 Cosmic variance 

In order to test for Gaussianity a simple approach is to measure the variance of the 
power spectrum Ci which contains the information of 4-point correlations. If the 
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Figure 5.5: Deviation from Distribution 

Plots of the distributions of = {21 + l)Ci/Ci for the Weeks models in 
comparison with the distributions with 2Z + 1 degree of freedom (solid 
curves) . The distributions are calculated using the eigenvalues obtained by 
the DBEM and the Gaussian approximation for the expansion coefficients. 
The Monte-Carlo simulations are carried out based on 200 realizations of 
the initial Gaussian fluctuations $;/(0), and 200 realizations of the base 
points. 



Figure 5.6: Cumulative Distribution 
Plots of 1 — P{Z) {P{Z) is the cumulative distribution function) which 
gives the probability of observing X > Z. The solid curves correspond to 
1 -P(X^) for the Weeks model 1^^ = 0.2, / = 5 (left) and l = lb (right). The 
dashed curves correspond to 1 — P{x^) of the Gaussian model. 



expansion coefficients him of the temperature fluctuation in the sky are Gaussian, 
X^ = (2/ + \)Ci/Ci must obey the distribution with 2/+1 degrees of freedom. In 
the case of CH models with small volume, the distribution of him for a homogeneous 
and isotropic ensemble has vanishing skewness but positive kurtosis. The non- 
Gaussian features in the distribution of appears as the slight shift of the peak 
of the distribution to the center (zero) and the slow convergence to zero for large 
amplitudes (flgure |^ ) . 

Roughly speaking the cosmic variance consists of the "initial" variance owing 
to the uncertainty in the initial conditions and the "geometric" variance owing to 
the uncertainty in the location and the orientation of the observer. The "initial" 
variance is equal to the "standard" cosmic variance ACi/Ci = ^2/21 + 1 while the 
additional "geometric" variance is smaller than the "initial" variance. The effect of 
the "geometric" variance is prominent for low density models since the relative size 
of the spatial hypersurface to the present observable region becomes small (flgure 
|5.7|) . Enhancement in the large-angle initial power ("blue-shifted power") also leads 
to a large "geometric" variance since the effective degree of freedom of the relevant 
fluctuations is decreased. 



Figure 5.7: Cosmic Variance 

Excess cosmic variances owing to the global inhomogeneity in the back- 
ground geometry. The ratios of the standard deviation in the angular 
power spectrum AC; for the Weeks models relative to that for the Gaus- 
sian model AC/ (Gauss) are plotted. In order to compute AQ Monte-Carlo 
simulations based on 200 realizations of the initial perturbation $i/(0) and 
200 realizations of the position of the observer have been carried out. 
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5.2.2 Topological measurements 

Topological measures:total area of the excursion regions, total length and the genus 
of the isotemperature contours have been used for testing Gaussianity of the temper- 



ature fluctuations in the COBE DMR dataP9|, |100|| . In integral geometry, they are 
known as the Minkowski functionals which characterize the morphological property 
of the fluctuations that are additive and invariant under translations and rotation. 
Let us first summarize the known results for Gaussian fields (see [|101| , |102|] ). 

The genus G of the excursion set for a random temperature field on a connected 
and simply-connected 2-surface can be loosely defined as 

G = number of isolated high-temperature connected regions (5.19) 
— number of isolated low-temperature connected regions. 

For instance, for a certain threshold, a hot spot will contribute +1 and a cold 
spot will contribute —1 to the genus. If a hot spot contains a cold spot, the total 
contribution to the genus is zero. 

The genus can be represented as the integration of the local properties of the 
field which can be easily simulated by a computer. From the Gauss-Bonnet theorem, 
the genus of a closed curve G being the boundary of a simply-connected region Qc 
which consists of arcs with exterior angles ai,a2, ...aN can be written in terms 
of the geodesic curvature and the Gaussian curvature K as 



G 

271 



/ Kgds + ai+ KdA 



(5.20) 



For a random field on the 2-dimensional Euclidean space E"^ where the N arcs are all 
geodesic segments (straight line segments), K and Kg vanish. Therefore, the genus 
is written as 



1 ^ 

Gi?2 = ^5Z«i (5.21) 



=1 



The above formula is applicable to the locally fiat spaces such as E^ x and 
which have E"^ as the universal covering space since K and Kg also vanish in these 
spaces. In these multiply connected spaces, the naive definition Eq.( |5.20|) is not 



correct for excursion regions surrounded by a loop which cannot be contracted to a 
point. 

In order to compute the genus for a random field on a sphere 5*^ with radius 
equal to 1, it is convenient to use a map ip-.S'^ — {pi\ — {p2}^ x (0, vr) defined as 

: (sin6'cos0, sin6'sin0,cos6') ^ (0, 61), < < 27r, < ^ < vr, (5.22) 



where pi and p2 denote the north pole and the south pole, respectively. Because 
X (0, vr) can be considered as locally flat spaces (0, 9) with metric ds^ = dO"^ + d(f)'^ 
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which have boundaries ^ = 0,7r, the genus for excursion regions that do not contain 
the poles surrounded by straight segments in the locally flat (0, 9) space is given 
by Eq. (|5.21|) . It should be noted that the straight segments do not necessarily 
correspond to the geodesic segments in S"^. If a pole is inside an excursion region 
and the pole temperature is above the threshold then the genus is increased by one. 
If the pole temperature is below the threshold, it does not need any correction. Thus 
the genus for the excursions is 

i 

where is the exterior angles at the intersection of two straight segments in the 
(0, 9) space and Np is the number of poles above the threshold. 

Now consider an isotropic and homogeneous Gaussian random temperature field 
on a sphere S"^ with radius 1. Let {x,y) be the local Cartesian coordinates on 
5"^ and let the temperature correlation function be C(r) =< (AT/T)o(AT/T),. > 
with r^ = x^ + y"^ and Cq = C{0) = a"^, where a is the standard deviation and 6*2 = 
— {d?C/dr'^)r=Q- Then the expectation value of the genus for a threshold AT /T = ucr 
is given as ||102 

<Gs^>= .[^%e-^'/' + erfc f ^ Y (5.24) 



where erfc(x) is the complementary error function. The first term in Eq. (|5.24| ) is 
equal to the averaged contribution for the excursions which do not contain the poles 
while the second term in Eq. (|5.24|) is the expectation value of Np. 

The mean contour length per unit area for an isotropic homogeneous Gaussian 



random field is flOl, 102 



and the mean fractional area of excursion regions for the field is the cumulative 
probability of a threshold level, 

<a>=^erfc(-^^, (5.26) 



1 



1 fC,\~' 



which gives the second term in Eq. ( |5.24D 



The CMB anisotropy maps for the Weeks and the Thurston (adiabatic) models 
are produced by using eigenmodes k < 13 and angular components 2 < Z < 20 for 
f^o = 0.2 and 0.4. The contribution of higher modes are approximately 7 percent 
and 10 percent for Qq = 0.2 and 0.4, respectively. The initial power spectrum is 
assumed to be the extended Harrison-Zel'dovich one. The beam-smoothing effect is 
not included. For comparison, sky maps for the Einstein-de-Sitter model with the 
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Figure 5.8: Simulated CMB Sky Map 

Contour maps of the CMB (not smoothed by the DMR beam) for 

the Thurston model = 0.4 and a flat (Einstein-de-Sitter) Harrison- 

Zel'dovich model Ci oc + 1)) in which all multipoles I > 20 are re- 
moved. 



Harrison-Zel'dovich spectrum Ci oc !/(/(/ + 1)) are also simulated. 

In order to compute the genus and the contour length for each model, 10000 
CMB sky maps on a 400x200 grid in the {4>,6) space are produced. The contours 
are approximated by oriented straight segments. The genus comes from the sum of 
the exterior angles at the vertices of the contours and the number of poles at which 
the temperature is above the threshold. The total contour length is approximated 
by the sum of all the straight segments. Typical realizations of the sky map are 
shown in figure [57^ . 

One can see in figure ^.9| that the mean genera and the mean total contours for 
the two CH models are well approximated by the theoretical values for the Gaussian 
models. This is a natural result since the distribution of the expansion coefficients 

is very similar to the Gaussian distribution in the modest threshold levels. On 
the other hand, at high and low threshold levels, the variances of the total contour 
lengths and the genera are much larger than that for the Gaussian models, which 
can be attributed to the positive skewness in the distribution function of 6;^- The 
excess variances for the Weeks model with i7o = 0.2, 0.4 compared with the Gaussian 
flat Harrison-Zel'dovich model are observed at the absolute threshold level approx- 
imately \i>\ > 1.4 for genus and \i>\ > 0.7 for total contour length, respectively. The 
reason why the non-Gaussian signature appeared in the variance of the topological 
measures rather than the means is explained as follows: Suppose that the tempera- 
ture map consists of hexagonal pixels whose distance between the center of any two 
adjacent pixel is A. From the additivity of Minkowski functionals, one can show 
that they can be decomposed into the components for each pixel and the intersec- 
tions of adjacent pixels which allows the explicit calculation of means and variances 
of them ||103|| . For instance, the mean of the total length per unit area for the map 
is given by 

<s>=^{Pi{u)-P2{u;A)) (5.27) 



where 



poo poo 

Pn{^, Tij) = dXi... dXn Pn{Xl, . . . , X„; Vij) (5.28) 

J J u 



that describes the probability of AT/T > z/ at points (xl, . . . ,x„) with distance r^j 
for each pair (xj,Xj) and Pni^^-, ■ ■ ■ i^n'i'^ij) is the n-point distribution function in 
real space. For "smooth" random fields, the limit of (|5.27|) as A ^ gives the mean 



value. In fact, for isotropic random Gaussian field, one can show that the limit of 
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( |5.27 ) gives ( 5.25 ). The variance of s is written as 

4 \2 



Var(s) = ]im(^) [P2(/^; Ha) - 2P3(^^; A, ri2) + P4(^^; A, ns)] 



Sl,S2 



■{^y[Piiu)-P2iu;A)r, 



(5.29) 



where is the total number of pixels and Si denotes a pixel and = \s2 — si\. 
Psi^u; A, ri2) represents the distribution of AT /T>u for 3 points configured as equi- 
lateral triangles with sides ri2 and A. Similarly, for genus one can show that the 
mean and the variance are written in terms of n-point distributions n < 3 and n < 6, 
respectively. Thus the variance has information of much higher-order correlations. 
It is concluded that the non-Gaussian signature is relevant to the higher-order cor- 
relations (n>3) in real space. 

As is well known, the COBE DMR data excludes grossly non-Gaussian models 
100|| . However, one should take account of a fact that the signals in the 10" 

This 



99 



104| 



smoothed COBE DMR 4-year sky maps are comparable to the noises 
makes it hard to detect the slight non-Gaussian signals in the background fluctua- 
tions. Because the mean behavior of the topological measures in CH models is well 
approximated by those of Gaussian models, the COBE constraints on CH models 
cannot be so stringent. In fact, some recent works using different statistical tools 
have shown that the COBE DMR 4-year sky maps are non- Gaussian ||, [Tol, Hol 



although some authors cast doubts upon the cosmological origin of the observed 
non-Gaussian signals ||107| , |108|| . If any non-Gaussian signals in the CMB on large 
angular scales are detected by the future satellite missions, then the non-trivial 
topology of the spatial geometry surely gives a natural explanation of the origin. 



Figure 5.9: Minkowski Functionals 
The means of genus and total length of isotemperature contours averaged 
over 100 realizations of the initial fluctuations and 100 realizations of the 
base points and itlo" run-to-run variations at 27 threshold levels for the 
Weeks model with VL^ = 0.2. The dashed curves denote the mean values 
for a corresponding Gaussian model. The solid curves denote the mean 
values for a Gaussian model that are best-fitted to the average values. 
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Figure 5.10: Non-Gaussian Signature in Variance 

The variances of genus and total length of isotemperature contours for 
the Weeks models (dashed and dashed-dotted curves) for a homogeneous 
and isotropic ensemble in comparison with those for the infinite flat model 
Cljn = 1-0 (dashed-dotted-dotted curve) with Gaussian initial fluctuations. 
At large threshold level u, the deviation from the Gaussian model is con- 
spicuous. 
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Chapter 6 
Summary 



In this thesis, we have investigated the CMB anisotropy in closed flat and closed 
hyperbolic models with multiply connected spatial hypersurface and obtained con- 
straints on these models using the COBE-DMR data and discussed various observ- 
able signatures of non-trivial topology in the CMB which can be tested by using 
more accurate data which will be supplied by the future satellite missions. 

In chapter 2 mathematical property of closed 3-manifolds with constant curva- 
ture K of three types, namely flat(iir = 0), spherica^/T > 0) and hyperbolic(ii' < 0) 
geometry has been reviewed. 

In chapter 3 we have formulated two types of numerical methods, namely the di- 
rect boundary element method (DBEM) and the periodic orbit sum method (POSM) 
for computing mode functions of the Laplace-Beltrami operator for closed hyper- 
bolic(CH) manifolds and have investigated the statistical property of mode func- 
tions and the relation between the low-lying eigenvalues and various diffeomorphism- 
invariant geometric quantities. The numerical accuracy of the DBEM is much better 
than that of the POSM but it needs mesh generation for each model. On the other 
hand the POSM is much suited for analyzing a large samples of spaces (manifolds 
and orbifolds) which needs only the fundamental group of the space. 

In section 3.1 we have seen that the low-lying eigenmodes that are continued onto 
the universal covering space of the two smallest CH manifolds are well described by 
a set of random Gaussian numbers. Finding such property for the low-lying modes 
is rather surprising since they are relevant to purely quantum nature. The descrip- 
tion breaks only for cases where the symmetries of the mode functions which are 
equivalent to the symmetries of the manifold and perhaps of the finite-sheeted cover 
of the manifold arc apparent. The origin of the pscudo random Gaussianity of the 
mode functions has not yet been well understood. The "random" property of a large 
sphere in the universal covering space may give a clue for finding the origin. We 
have conjectured that the pseudo random Gaussian behavior is universal for all CH 
spaces which can be used for simulating the CMB anisotropy in {l,m) spaces. 

In section 3.2, we have numerically analyzed the length spectra and the eigen- 
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value spectra of CH manifolds with small volume. First, the asymptotic behavior in 
the classical staircase is found to be consistent with the known analytical formula 
which does not depend on the topology or symmetry of the manifold although the 
arithmeticity of the manifolds makes a difference for behavior in the multiplicity 
number. Next, we have applied the trace formula to these manifolds and obtained 
consistency with those by the DBEM. No supercurvature mode are found for a to- 
tal of 308 CH 3- manifolds (volume < 3). We need a further investigation to check 
whether small manifolds that are sufficiently similar to cusped manifolds can sup- 
port supercurvature modes. Thirdly, the first eigenvalues are compared to diameter, 
volume and the shortest length of the periodic orbits. The numerical results imply 
the existence of much shaper bounds for the first eigenvalues in terms of diameter. 
Some fitting formulae have been introduced and their validity has been checked. CH 
3-manifolds can be roughly divided into two categories: "slightly anisotropic" and 
"almost anisotropic" ones. The former has not any very short periodic orbits while 
the latter has. For example, manifolds which are very similar to the original cusped 
manifold are belonging to the latter category. It is found that the deviation of the 
spectrum from the Weyl asymptotic formula for these manifolds is conspicuous even 
for manifolds with small volume. Finally, the global "anisotropy" in the spatial ge- 
ometry has been measured by C-function and the spectral distance for 3 examples of 
CH 3-manifolds. It is found that the angular power spectra in the CMB are greatly 
affected by the globally anisotropic structure in the spatial geometry. 

In chapter 4, we have explored the current constraints on the closed fiat and 
hyperbolic models using the COBE-DMR data. In section 4.1 and 4.2, we have 
derived the Sachs- Wolfe formula which describes the CMB anisotropy based on the 
linear perturbation theory in a gauge-invariant way. In section 4.3, we have investi- 
gated the effect of the non-trivial topology on the angular power spectra. Assuming 
adiabatic initial perturbation with scale-invariant spectrum (n = 1), a prominent 
suppression occurs for a "standard" closed ffat toroidal model with = 1-0 since 
fluctuations beyond the size of the cube at the last scattering are strongly sup- 
pressed. However, for low matter density models such a strong suppression does not 
occur since one cannot ignore the contribution from the (late) ISW effect caused by 
the decay of the gravitational potential at the A or curvature dominant epoch. The 
crucial point is that they are produced at late time well after the last scattering. If a 
fluctuation on scale smaller than the fundamental domain is produced at a point suf- 
ficiently nearer to us, then the corresponding angular scale becomes large provided 
that the background geometry is fiat or hyperbolic. Therefore, the large angular 
fiuctuations produced late time "survive" even if the size of the cell is finite. On 
the other hand, slight suppression in large-angle temperature correlations in such 
models explains rather naturally the observed anomalously low quadrupole which 
is incompatible with the prediction of the "standard" FRWL models. In section 
4.4, we have carried out the Bayesian analyses using the COBE-DMR 4year data. 
In the case of fiat topology, the discrete eigenmodes have "regular" features which 
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leads to significant correlations in a/m's which describes the fluctuations in the (/, m) 
space. Even if marginalized over the orientation, the correlations do not completely 
disappear. The analysis using only the angular power spectrum is not enough since 
the background geometry is globally inhomogeneous or anisotropic. However, even 
including the non-diagonal elements in the correlations the constraint is still less 
stringent for low matter density models. This is because the physical size of the 
relevant fluctuations which are produced at late time are much smaller than the size 
of the fundamental domain. In the case of hyperbolic topology in which the discrete 
eigenmodes are "chaotic" the correlations in aim's for a given / almost disappear 
if one takes an average over the position of the observer. As we have seen, the 
likelihood of CH models is very sensitive on the position of the observer. Although, 
for a large number of choices the fits to the COBE data are bad, in some places 
the fits are much better than those of the FRWL models. Because all CH manifolds 
are globally inhomogeneous p7| , in order to constrain the models, one must compare 
the expected sky map at every place and every orientation to the data. Even after 
marginalized over the position and orientation of the observer, the likelihoods are 
still comparable to that of the infinite counterparts. It is natural that the previous 
analyses on CH models |HT]] have given stringent constraints since only fluctu- 
ations at a point where the injectivity radius is maximal with only 24 orientations 
were compared to the COBE data. Nature wouldn't choose such a special point 
where the symmetry of the fluctuations is maximal as our place in the universe. 

In the last chapter, we have investigated the observable signatures of the spatial 
non-trivial topology, namely the periodical structure and the non-Gaussianity in the 
CMB. In section 5.1 we have re-examined the method of finding matched circles. We 
have crudely estimated the detectability of the periodical patterns in the CMB sky 
map assuming that the signal and noise obeys the Gaussian statistics as in but 
using slightly different statistics. The detectability depends on the signal-to-noise 
ratio and the angular resolution of the sky map as well as the number of circles and 
their angular radius. We have also considered the effect from the "background noise" 
owing to the Doppler effect on the scale of acoustic oscillations which might smear 
the perfect identical patterns. The chance of detecting the non-trivial topology us- 
ing the COBE data has found to be very low but we can expect a good chance of 
the first discovery of the imprint of the non-trivial topology using CMB maps with 
much better signal-to-noise ratio and the angular resolution which will be supplied 
by the future satellite missions, namely MAP and Planck. In section 5.2 we ex- 
amined the non-Gaussianity of the CMB temperature fluctuations. Inhomogeneous 
and anisotropic Gaussian fluctuations for a particular choice of position and orien- 
tation are regarded as non-Gaussian fluctuations for a homogeneous and isotropic 
ensemble. For CH models, the distribution of the expansion coefficients bim of the 
2-dimensional temperature fluctuations has positive kurtosis but vanishing skewness 
provided that the initial fluctuations are Gaussian. The non-Gaussianity leads to 
a slightly larger cosmic variance where the possibility of having a rare fluctuation 
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is high. The non-Gaussianity is much prominent if one considers the statistics of 
Minkowski functionals, namely the total area, the total length and the genus of 
isotemperature contours. The non-Gaussian signature would appear as large vari- 
ances of statistics at high and low threshold levels on large angular fluctuations. If 
such a signature should be found in the forthcoming GMB sky maps, then it would 
certainly be a strong sign of the non-trivial topology. Analyzing such statistics for 
other closed multiply connected models with low matter density, namely, flat and 
spherical models with non-trivial topology will be also interesting issues. Even if we 
failed to detect the identical patterns, the imprint of the non-trivial topology -the 
effect of the "finiteness"- could be still observable by measuring such statistics. 

The determination of the global topology of the universe is one of the key issues 
of the modern observational cosmology. The observation of the cosmic microwave 
background, the distribution of clusters and QSOs and other astronomical objects 
in the distant places provides us precise information about the global topology of 
the universe as well as other cosmological parameters, such as matter contents, pri- 
mordial spectrum, Hubble parameters and so on. In the first decade of the 21st 
century, we will be able to answer the fundamental and the simple question whether 
the universe is finite or "sufficiently big" . In the former case, it surely gives a great 
leap in our comprehension of the universe. Even in the latter case, our quest for 
searching the global topology will continue until we obtain a well established unified 
theory which integrates the microscopic and the macroscopic worlds. 
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APPENDIX A. BOUNDARY INTEGRAL EQUATION 



Appendix A 

Boundary integral equation 



Here we derive the boundary integral equation ( |3.12| ) in section 3.1. For simplicity, 
we prove the formula in 3-spaces. 

First, we start with Eq.( |3.11| ) with dimesnsion M = 3. Although the integrand in 
Eq. (|3.11|) is divergent at x = y G dfl, the integration can be regularized as follows. 



Let us draw a sphere with center y G dil with small radius e and let be the outer 
spherical boundary and a and (3 be the internal solid angle and external solid angle 
as shown in figure 



u{y) + / G^(x, y)V,u ^dS' - [ (V,Gs(x, y))^ ^dS' = 0. (A.l) 
The singular terms in Eq.( |A.l| ) can be separated from non-singular terms as 

lim / GE{^,y)V^u^dS' = I GE{^,y)V,u^dS' 
•^^0 Jar+r^ Jar 

+ lim / GE{^,y)V^u^dS\ 

lim [ {V,GEi^,y))u^dS' = [ {V,GE{^,y))u^dS' 
Jar+Ve Jar 

+ lim / {V,GE{^,y))u^dS\ 
Jr, 

If e is sufficiently small, the region enclosed by F^ can be approximated as an Eu- 
clidean subspace. In this region, the asymptotic form of the free Green's function 
Ge takes the form 

^ / N exp(ikd) 1 ik ,^ . . x 

where d is the Euclidean distance between x and y. Taking the spherical coordinates 
(e, 9, 0) with center y, the singular terms in Eq. ( |A.2| ) are estimated as 

lim [ GE{^,y)VM^)V9dS' = lim- [ -^^^e^dQ = 0, 
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Figure A.l: Boundary Integral 



lim / iV,GEi^,y))u{^)V9dS' 



lim 

e — >■ 
I 



1 9e(x) 
Ane"^ dn 



(A.3) 



where dQ denotes the infinitesimal solid angle element. Taking the limit e 
Eq.( |A.l| ), we have the boundary integral equation in 3-spaces, 



in 



—a{y)u{y) 



GE{^,y)V.u^dS'- {V,GE{^,y))u^dS' = 0, (A.4) 
an Jan 



where a(y) denotes the internal solid angle at y. If the boundary is smooth at y, 
a(y) is equal to 2tt which gives the coefficients 1/2 in Eq. (p.l2|) . Similarly, one can 
prove the formula for M = 2 and M > 3. 
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APPENDIX B. ALGORITHM OF COMPUTING LENGTH SPECTRA 



Appendix B 

Algorithm of computing length 
spectra 

For a given Dirichlet domain D, SnapPea computes 

1. Neighboring copies of D (tiles) recursively and stores the corresponding ele- 
ments g of the discrete isometry group T. If it has already been found out, it 
is discarded. The computation proceeds until for all the neighborhoods of gD, 
d{x,hx) > cosh~^ (cosh i? cosh Z/2) satisfies where hD's are neighborhoods of gD{R 
is the spine radius and h is an element of F). Since there is no g other than identity 
where all h's satisfy d{x,hx) > d{x,gx), this algorithm will not miss any tiles gD 
where d{x,gx) < 2cosh~^(coshi?cosh//2). 

2. A hst of geodesies for all g's — {g} where l:the real part of length is not zero and 
less than I; 2:the distance from x to the geodesic is at most R. 

3. The conjugacy class g' = hgh~^ or its inverse for each g where h is an ele- 
ment of {g}. If an identical complex length is found, the complex length which 
corresponds to g' is omitted from the list of geodesies. If the complex length of g is 
conjugate to that of g~^, the geodesic is topologically a mirrored interval, otherwise 
it is a circle. 

4. Multiplicity of geodesies. If a pair of geodesies with two complex lengths be- 
ing identical within an error range is found, the multiplicity number is increased by 
one. 
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Appendix C 

Estimate of volume in terms of 
diameter 



Suppose a CH manifold M which resembles the original cusped manifold Mf. with 
one cusp. Let us divide into two parts, the neighborhood of a cusped point K^, 
and the complementary part Kcq. Similarly, one can divide M into K and Kq where 
Kq and K corresponds to a "thin" part. Since the neighborhood of a cusp 
is represented as a "chimney" (but having infinite length) in the upper half space 
coordinates { ), K can be well approximated by an elongated box defined by 

(— Aa;/2 < xi < Aa;/2,— Aa;/2 < X2 < /S.x/2,xzq < x^ < x^i). Then the physical 
length (^diameter) d of X in the direction is given \yy d — ln(x3i/x3o). On the 
other hand, the volume v of K satisfies 



{Axf ( 1 



r)> (C.l) 



which gives the ratio of the volume v to the volume of Vao of -f^ci vjvaa = 1 — 
(3:^30)^/ (2^31)^ = 1 ~ exp(— 2(i). If we approximate the diameter 0? of M as d = 0? + do 
where do is the diameter of Kq^ then we finally have the ratio of the volume v — v^^v 
of M to the volume Vc — vq^ Voo of Mc, 

V exp{-2{d - dp)) 

— = 1 — — , d = vo (C.2) 

Vc + 1 

where Vq denotes the volume of Kq. 



Figure C.l: Estimate of Volume of "thin" Part 
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